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^ ■ Abstract 
cn 

^ ! We compute the axial anomaly of a Lifshitz fermion theory with anisotropic scaling 
lO ■ z = 3 which is minimally coupled to geometry in 3 + 1 space-time dimensions. We 
2 ■ find that the result is identical to the relativistic case using path integral methods. An 
independent verification is provided by showing with spectral methods that the rj- 
invariant of the Dirac and Lifshitz fermion operators in three dimensions are equal. 
, Thus, by the integrated form of the anomaly, the index of the Dirac operator still ac- 
'. counts for the possible breakdown of chiral symmetry in non-relativistic theories of 



gravity. We apply this framework to the recently constructed gravitational instanton 
backgrounds of Horava-Lifshitz theory and find that the index is non-zero provided 
that the space-time foliation admits leaves with harmonic spinors. Using Hitchin's 
construction of harmonic spinors on Berger spheres, we obtain explicit results for the 
index of the fermion operator on all such gravitational instanton backgrounds with 
SU{2) X U{1) isometry. In contrast to the instantons of Einstein gravity, chiral sym- 
metry breaking becomes possible in the unimodular phase of Hof ava-Lifshitz theory 
arising at A = 1/3 provided that the volume of space is bounded from below by the 
ratio of the Ricci to Cotton tensor couplings raised to the third power. Some other 
aspects of the anomalies in non-relativistic quantum field theories are also discussed. 
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1 Introduction 



Since the discovery of the axial anomaly and its profound success in explaining the 
two-photon decay mode of the meson 111 13, there has been a lot of activity in the 
study of anomalous symmetries in quantum field theory and their geometrical and 
topological origin (see, for instance, [3] for a recent brief overview of the subject and 
its history) as well as in model building The axial anomaly is accounted pertur- 
batively by the one-loop fermion contribution (triangle diagram) of the two-photon 
coupling to an axial current, as in spinor electrodynamics, but it also has an alternative 
description in terms of the non-invariance of the Euclidean path integral measure un- 
der chiral rotations of the fermions coupled to a background gauge field [5] (see also 
the textbook [6]). This framework was soon afterwards generalized to fermions cou- 
pled to a background metric by computing the gravitational corrections to the axial 
anomaly associated to the two-graviton decay mode of the meson (Zl El [3, which 
also enjoys a path integral derivation EEI as for the case of gauge fields (see also ||T0| 
for the structure of gravitational anomalies in general dimensions). These results refer 
to the so called local form of the anomaly 

Another important development is related to the integrated form of axial anomaly 
which can be non-zero when the background gauge field is non-Abelian and topo- 
logically non-trivial. In this case, the four-dimensional Dirac operator exhibits nor- 
malizable zero modes with unequal number of positive and negative chirality solu- 
tions (denoted by n± respectively) in the background of an tnstanton configuration. 

^It is beyond the scope of the present work to give a comprehensive account of the many important 
contributions made over the years by several groups of people, as we will only focus on the axial 
anomaly in its simplest form. Several expository contributions on the subject can be foimd, for instance, 
in 14J. 
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Thus, there is chiral symmetry breaking induced by instantons that may in turn lead 
to baryon and lepton number violation in Yang-Mills theories | fTT| . This novel pos- 
sibility is a physical manifestation of the Atiyah-Singer index theorem [[T2]| (but see 
also the textbooks IITSl [T4|| ) that equates the index of the four-dimensional Dirac op- 
erator, namely n+ — n_, to the topological charge k of the background gauge field on 
the compactified Euclidean space-time S^. (In the actual instanton solutions one has 
n_|_ = /c > and n_ = 0, whereas for the anti-instantons the situation is reversed, i.e., 
U- = —k > and n_|_ = 0.) 

Likewise, the integrated form of the gravitational correction to the axial anomaly 
is a topological invariant given by (one eighth of) the Hirzebruch signature of space- 
time, which can be non-zero on topologically non-trivial gravitational backgrounds 
and may lead to similar results for chiral symmetry breaking by the index theorem 
of Dirac operator. This is certainly true for the Dirac operator on K3, which is the 
unique compact self-dual gravitational instanton without boundaries Ifl5l [T6l : in this 
case there are two covariantly constant negative helicity modes and no positive helic- 
ity modes so that n+ — ti- = —2. However, gravitational instantons with boundaries 
do not admit any normalizable zero modes of the Dirac operator, provided that the 
Ricci scalar curvature of space-time is non-negative (as it is by Einstein's field equa- 
tions when the cosmological constant is non-negative), and by a general result based 
on Lichnerowicz's theorem fTTj ] the index of Dirac operator vanishes on such back- 
grounds |[T8|. An alternative derivation of the same result is based on the Atiyah- 
Patodi-Singer index theorem |fT9l |20l [T3ll taking into proper account the boundary 
terms of space-time from the asymptotic regions ||2ni22ll23ll24ll25ll26ll27ll28B . Thus, 
although the signature of space-time is non-zero, the index of the Dirac operator turns 
out to be zero for all gravitational instanton backgrounds with boundaries, in which 
case there is no chiral symmetry breaking induced by quantum tunneling effects. For 
a systematic exposition of all these matters we refer the interested reader to the report 
||29| . More explanations will also be given in the text. 

In this paper we examine the chiral anomaly for non-relativistic fermions of Lif- 
shitz type that are coupled to background fields in four space-time dimensions. Our 
investigation was prompted by the current activity on non-relativistic theories of grav- 
ity with anisotropic scaling, serving as alternatives to the ultraviolet completion of 
Einstein theory of gravitation with higher order spatial derivative terms [|30i|3T|. The 
resulting theory sacrifices relativistic invariance for perturbative renormalizability at 
short distances and despite criticism related to the presence of an unphysical mode 
and the difficulty to flow to ordinary gravity at large distances (see, for instance, 
Il32l |33ll and references therein, as well as EH]), it provides an interesting framework 
that should be given the benefit of the doubt to play a useful role in quantum gravity 
(see also | |35| for a more recent proposal that defies some of the criticism). Coupling 
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the theory to Lifshitz fermions with the same anisotropic scaling in space and time, as 
for gravity, is a step that we take in this work and which, apparently, has not been con- 
sidered in the literature so far. The computations we present in the following provide 
the gravitational contribution to the axial anomaly at a Lifshitz point and, as such, 
they generalize rather naturally the results on the axial anomaly for Lifshitz fermions 
coupled to gauge field backgrounds [36|. In either case, the local form of the anomaly 
turns out to be identical to the relativistic case. These results are surprising at first 
sight, but clearly they should be related to the topological character of the anomaly 
that is proportional to the characteristic classes Tr(f A f ) and Tr{R A R) for gauge 
and metric fields, respectively. They are also quite general since they treat the gauge 
and metric fields as background without reference to their field equations. As such, 
they provide the local form of the axial anomaly at a Lifshitz point and they are in- 
dependent of the particular form of the bosonic action. This can make them useful 
for applications in other fields of current research that incorporate ideas of quantum 
criticality. 

The integrated form of the axial anomaly at a Lifshitz point is also identical to the 
relativistic case, and, therefore, it equals the index of the relativistic Dirac operator in 
a given background. This will be shown explicitly by proving an index theorem for 
the Dirac-Lifshitz operator and by computing the //-invariant of the operator on the 
three-dimensional leaves of space-time foliation. These results can be regarded as an 
alternative proof that the anomalies are the same in both theories. The index is also 
inert to the back-reaction of fermions on the gauge and metric fields. Thus, it makes 
good sense to compute the index on the background of an instanton solution, as in 
the relativistic case. Here, we pass directly to (3 + 1) -dimensional Hofava-Lifshitz 
gravit}J§ coupled to Lifshitz fermions with anisotropic scaling z = 3 using its simplest 
form (said to satisfy the detailed balance condition [301 |3T1) that exhibits instanton 
solutions in the Euclidean regime, p8|. We will find explicit examples that are capa- 
ble to produce a non-vanishing index and, hence, lead to chiral symmetry breaking 
by gravitational instanton effects. These instantons are very different in nature from 
the gravitational instanton solutions of ordinary Einstein theory in that they are not 
self-dual spaces; they resemble more the instanton solutions of point particle systems 
that interpolate continuously between different vacua as the Euclidean time runs from 
— oo to +00 and they have the topology IR x E3 (e.g., E3 ~ in the typical examples 
that we will consider later). These configurations are also chiral in that they exist 
for one choice of orientation on Z3 and not for the other, but as it turns out by explicit 
computation they should be "sufficiently chiral" to allow for non-vanishing index (the 

^The other interesting case of a Lifshitz type gauge theory (in the spirit proposed in |37l ) coupled to 
Lifshitz fermions and the integrated form of the axial anomaly in the background of gauge fields will 
not be treated at all, for a good reason that will be discussed later, even though it looks more elementary. 
Hopefully, we will turn to this and related issues elsewhere. 
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precise meaning of these words will be given in due course). This possibility will be 
made available in the phase of Horava-Lifshitz theory that exhibits spatial unimod- 
ular invariance in sharp contrast to ordinary Einstein gravity that never allows for 
chiral symmetry breaking induced by instantons. The results can be regarded as an 
application of our instanton solutions Il38i to non-relativistic theories of gravitation. 

The material of this paper is organized as follows: In section 2, we present our re- 
sults for the axial anomaly of a Lifshitz fermion theory which is minimally coupled to 
gauge and metric backgrounds and show that in either case they are the same with the 
relativistic theory of Dirac fermions, hereby extending the results of Dhar et.al. ||36l 
to the gravitational form of the anomaly. In section 3, we present a short overview 
of Horava-Lifshitz gravity in 3 + 1 dimensions and describe its instantons in terms 
of eternal solutions of a certain geometric flow equation called Ricci-Cotton flow, fol- 
lowing [38|. We also summarize all instanton solutions with SU{2) isometry that will 
be used as working examples in the applications. In section 4, we compute the index 
of the Dirac-Lifshitz operator on these gravitational instanton backgrounds by taking 
proper account of the boundary terms. This index counts the difference of positive 
and negative chirality zero modes of a Lifshitz fermion and it turns out to be non- 
zero on certain instanton backgrounds in the unimodular phase of the theory, which 
arises when the parameter A of superspace assumes a particular value (to be explained 
later) that is special to the deep ultra-violet regime of the theory. In this context, we 
also obtain a critical value for the couplings, equivalently for the volume of space, be- 
yond which zero modes become possible. The index is always zero in all other phases 
of the theory. The index can become non-zero for geometrical rather than topolog- 
ical reasons (this is a special feature associated to the presence of harmonic spinors 
in three dimensions according to Hitchin |i39|, but see also |j40l), which explains the 
lower bound on the volume of space that will be obtained later. Although the results 
are derived using the class of S 11(2) symmetric solutions, they serve as prototype for 
the more general situation. In section 5, we compare our results with the index of the 
Dirac operator on gravitational instanton backgrounds of ordinary Einstein theory fo- 
cusing, in particular, to instantons with 51/(2) isometry for which direct comparison 
can be easily made. It is reconfirmed that the index is zero in that case. In section 
6, we present our conclusions and discuss some of the implications of anomalies in 
quantum field theories of Lifshitz type. Directions for future work are also outlined. 

Finally, there are four appendices with some technical material. Appendix A refers 
to the algebra of Dirac matrices in flat Euclidean space IR^ and contains some useful 
mathematical identities that involve the trace of products of gamma-matrices. They 
will all be needed for the evaluation of the axial anomalies. Appendix B contains some 
geometrical apparatus that is needed for the coupling of fermions to geometry and the 
evaluation of the anomaly based on the notion of geodesic interval and the associated 
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Synge-DeWitt tensors [1411 142| (but see also l|43ll ). Appendix C summarizes the Eu- 
clidean path integral method used to obtain the local form of the axial anomaly in the 
background of gauge and metric fields in relativistic quantum field theory, following 
ISHl]]. It sets up the notation and framework for carrying out the same computation 
in the non-relativistic theory of Lifshitz fermions. Appendix D summarizes some 
mathematical details related to harmonic spinors on equipped with homogeneous 
metrics, following ||39l |40l , which are needed for computing the index of the fermion 
operator on our gravitational instanton backgrounds. Explicit results are included for 
the case of Berger spheres that support gravitational instantons with SU{2) x Lr(l) 
group of isometries. 



2 Axial anomalies of Lifshitz fermions 

In this section, we consider the non-relativistic theory of Lifshitz fermions and dis- 
cuss some general properties including the symmetry of axial rotations in the massless 
case. We also set up the formalism and compute the quantum anomaly to the diver- 
gence of the corresponding axial current in close analogy with the relativistic Dirac 
theory. The generalization to fermions ip with mass m can be easily made by adding 
the classical term mip^^ip to the divergence of the axial current, as usual, but this will 
not be discussed further. 



2.1 Lifshitz fermions and axial symmetry 

The theory of Lifshitz fermions in 3 + 1 space-time dimensions is defined by the fol- 
lowing action 

S = j dtd^XsJ\deiG\ipiY^}i^ (2.1) 

that resembles that of a Dirac fermion ip. The conjugate field is defined, as usual, 
by i/7 = 1/7^7*^. The action is generally taken in curved space-time with local coor- 
dinates {t,x^) (assuming the presence of a privileged time direction) and metric G^v 
with signature — h ++, whereas 7^' are the standard Dirac matrices satisfying the 
anti-commutation relations 

[7^7^+ = 2G?'^ (2.2) 

The difference, however, comes in the definition of the operator D^,, which is of mixed 
order, 

-1 

Do = , D, = - (d,(-D2) + (-D2)D,) , (2.3) 
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where D^^ are the usual covariant derivatives expressing the minimal coupling of the 
fermion field ip to the background geometry and/ or background gauge fields and 
D2 = DiDK Clearly, this is a non-relativistic theor}|§ with anisotropic scaling z = 3, 
which is seen by assigning scaling dimensions [L] = — 1, [T] = —3 and [ip] = 3/2. 
Then, the analogue of the Dirac equation and its conjugate equation read as follows, 

j^lD^tp = 0, {ID^,tp)j^ = 0. (2.4) 



The action (|2.1[) is naturally defined on space-times that are topologically of the 
form R X E3 using the privileged time direction and the corresponding foliation by 
three-dimensional spatial slices L3. Then, the space-time metric G^/ takes the general 
form 

ds^ = -N^dt'^ + gij l^dx' + NUt^ (^dxi + N'dt^ (2.5) 

using lapse and shift functions N{t,x) and Ni{t,x), respectively, and gij{t,x), as in 
the ADM (Arnowitt-Deser-Misner) formulation of Einstein gravity (see, for instance, 
1141). The metric on E3 is gij providing the induced metric on the leaves of the fo- 
liation, i.e., gij = Gij, whereas G^^ = g'J - N^N' with = N'/N. We also have 



detG I = Ny^detg. Here, we confine ourselves to the so called projectable case, 
meaning that the lapse function N associated to the freedom of time reparametriza- 
tion is restricted to be a function of t alone, whereas the shift functions N, associated 
to diffeomorphisms of E3 can depend on all space-time coordinates. This is natu- 
ral in all Lifshitz theories that do not exhibit general coordinate invariance but only 
a restricted space-time symmetry associated to foliation preserving diffeomorphisms 
t = i, x^ = x^{t,x). Furthermore, in order to simplify certain aspects of the presenta- 
tion, and, also, in view of the applications we will make in subsequent sections, we 
adopt the choice of lapse and shift functions 

N{t) = l, Ni{t,x)=0 (2.6) 

unless stated otherwise. Setting N{t) = 1 amounts to using proper time in the pro- 
jectable case and Ni{t,x) = amounts to having G'^ = g'^ so that the spatial space- 
time indices can be lowered and raised using the three-dimensional metric on E3 with- 
out worrying about off-diagonal space-time terms. 

Next, we examine the symmetries of the classical action (|2.1|) . First of all there 
is a non-relativistic variant of the vector current }^{t,x) = ip'y^ip associated to the 
invariance of the action under xp — ^ exp{ia)ip and xf> — ^ exp{—ia.)if>, but its details 



■'Note that the theory of Lifshitz fermions we are considering here is rather different from that of 
non-relativistic fermions based on Galilean-covariant field theory (45l . The latter can be obtained 
from a relativistic massless fermion theory in one dimension higher by performing a Scherk-Schwarz 
reduction along a null direction and they exhibit different behavior. 
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will not be relevant to the present work. In any case, the vector symmetry is also 
present quantum mechanically, since there is no anomaly to obstruct its conservation 



law. Also, since the mass of the fermions is zero, the action (|2.1[) is invariant under 
chiral rotations, 

(2.7) 



using 75 = (z/24)e^i/KA7''7^7'^7'^ (it takes the standard form 75 = i^y^^^^'^^^ in terms 
of tangent space-time indices associated to a local Lorentz frame), as for the relativis- 
tic case. It follows that there is an associated axial current /g , which is classically 
conserved by virtue of the equations of motion. 



where its time and space components are chosen to be 



(2.8) 



and 



Ut,x) 



]l{t,x) = {p^\^xp 



[Djxp)jjj5iD'xp) + (D'xp)j}j5{Djxp) + xpYjsi-D^xp) + 



-D^xp)Yj5xp - xpyjj5{D'Djxp) - (D'D^f )7%t/^ 



(2.9) 



(2.10) 



These are the components of a real current that should be compared to the corre- 
sponding components of the axial current }^(t,x) = tpj^Jstp of the relativistic Dirac 
theory, which are, of course, much simpler. The situation is more analogous to the 
current conservation law associated to the time dependent Schrodinger equation in 
non-relativistic quantum mechanics. 

Before we proceed further, an important remark is in order regarding the unique- 
ness of the Dirac-Lifshitz operator. There is a factor ordering ambiguity in the con- 
struction of its spatial components because the covariant derivatives that we raise 
to the third power do not commute - their commutator is proportional to the field 
strength of the background field (it is the Riemann curvature in the case of a gravita- 
tional field). Another choice that could be made equally well is provided by 



= -DjDiDj , 



(2.11) 



in which case the spatial components of the associated axial current should be re- 
placed by 



r,\t,x) = {Djxp)yj5{Djxp) - xpj}j5{DjD'xp) - {DjD^xp)jjj5^ (2.12) 
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to ensure current conservation with respect to the new (prime) Dirac-Lifshitz equa- 
tion. Note, however, that this ambiguity is irrelevant for the purposes of the present 
work, since D, and differ from each other by a total derivative term of the back- 
ground field strength, which can not affect the local form of the anomaly, since its 
functional form is protected by topology. Throughout this paper we will adopt the 
particular choice (|2.3|) , even though the alternative choice (|2.11[) looks simpler as it 



admits a simpler current. This will not affect much the intermediate steps of the cal- 
culation of the axial anomal}]^. Finally, we note for completeness that other factor 
orderings also provide viable choices for the Dirac-Lifshitz operator (for instance, one 
may choose the Weyl-ordered third order operator = (D; + D-)/2), but they all 
have the same degeneracy. 

We can perform path integral quantization of the Lifshitz fermions coupled to 
background gauge or metric fields and examine the fate of axial current conservation. 
We will derive the corresponding axial anomaly following the method of Fujikawa 
EUS by taking the theory in the Euclidean domain. First, using the standard Noether 
procedure, the fermionic action transforms as follows under the chiral rotation (|2.7|) 
with parameter a{t,x), 

S^S + j dtd^x^J\deiG\ a{t,x)V^}^{t,x) . (2.13) 

In the quantum theory, this transformation is compensated by the Jacobian of the path 
integral measure for the fermions ('Dip)('Dip) and their combination yields the quan- 
tum correction to the axial current conservation law. 

In the Euclidean regime, where we will work from now on, the signature of space- 
time is assumed to be + + ++ (opposite to the convention often adopted 

in the literature). For this, we perform the Wick rotation t it and 70 270 so 
that (7^')^ = for all space-time indices. In the Euclidean regime, the Dirac op- 
erator i^^D^, as well as the Dirac-Lifshitz operator 27'' D^, we are considering here 
are both Hermitian, whereas 75 becomes 75 = — 7o7i7273 arid it anti-commutes, as 
usual, with the Dirac matrices, i.e., [75, 7?^]+ = 0, and (75)^ = 75. The Euclidean 
space 75 follows from its Lorentzian counterpart by letting 75 — ^ 275; this should be 
properly accounted when comparing (|2.13|) with the transformation of the fermionic 
measure that is most conveniently described in the Euclidean regime. The algebra of 
Dirac matrices in R"^ together with their trace identities that will be relevant for the 
computation of the axial anomalies are summarized in Appendix A. 



^Actually, in the paper |36| that only discusses the gauge field contribution to the axial anomaly of 
Lifshitz fermions, the issue of factor ordering is not addressed at all (not even the Hermiticity properties 
of the Dirac-Lifshitz operator) because this turns out to be irrelevant for the intermediate steps of the 
calculation, and, of course, for the final result. 
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The path integral over the fermions with Euclidean action is defined as 

Z = [ (r>i/;)(r>i/7)e-5E (2.14) 



setting Planck's constant equal to 1. For now and later use, we consider the complete 
set of eigen-functions cpn of the Hermitian operator i^^Y)^, 

{ij^n^)(pn{t,x) = \n(pn{t,x) , (2.15) 

which are normalized as 

j dtd^xV detG (p'l^{t,x)(pm{t,x) = Snm ■ (2.16) 

Then, a general fermion configuration is decomposed into eigen-states as ip{t,x) = 
Iln'^n(Pn{t>x) and xf>{t,x) = Ew (i, using independent elements fl„ andfc„ of the 
Grassmann algebra, whereas the fermionic path integral measure is formally written 
as {'Dxf)){Vxp) = Ylm bm Hn ^n- Under the chiral rotation (|2.7|) , the coefficients fl„ and 
bn can be easily seen to transform as Un — > Em Cnm^m and fc„ — > Em Cnmbm, where 

Cnm = J dtd^xVdetG cpl {t, x)e'<^'''^'^^ cpm {t, x) . (2.17) 

As a result, the fermionic path integral measure picks up a Jacobian factor and trans- 
forms as 

{V\p){V\p) {Vip)(Vip) exp ^-2i J dtd^xVdetG ci(t,x)Y^(pl{t,x)j5(pn{t,x)^ . 

(2.18) 

Consequently, the primitive local form of the axial anomaly follows by combining the 
result of transformations (|2.13|) and (|2.18)) in the Euclidean regime and it reads 



V^jf (i,x) = 2^<p+(i,x)75(p„(i,x) . (2.19) 



Note that all steps outlined above are identical to the path integral formulation of 
the axial anomaly for Dirac fermions followed by simply replacing by E)^. How- 
ever, the actual evaluation is more tricky and it requires computing the formal sum 
(trace of 75) shown in (|2.19|) over the complete set of eigen-functions of the interact- 
ing theory. In all cases, this sum is ill-defined, as it formally diverges, and appropriate 
regularization is required to extract the local form of the anomaly. The standard pro- 
cedure is to regularize the large eigen-values, i.e., | A„ | < (accounting also for their 
scaling dimension), and then obtain finite result for the right-hand side of equation 
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(|2.19|) by computing 



A{t,x) = lim 

M^-oo 



(2.20) 



Note the appearance of the anisotropy scaling parameter z in the cut-off M, which is 
z = 1 for Dirac fermions and z = 3 for the Lifshitz fermions we are considering here. 

The final expression for the anomaly depends upon the background gauge and/ or 
metric field via the corresponding field strength and it should be gauge invariant. It 
should also be topological density so that A (f, x) can be locally written as total deriva- 
tive term to account for the anomaly in the divergence of the axial current (|2.19|) and. 



as such, its form is very restricted (up to an overall factor). The actual terms that 
contribute to the evaluation of the anomaly depend crucially on the operator and the 
associated eigen-values and functions A„ and cpn that enter into (|2.20|) ; in a diagram- 



matic approach to the same problem different loop diagrams contribute to the answer, 
depending on the available couplings to the background fields. Nevertheless, the out- 
come is proportional to Tr(F A F) for the gauge and Tr(R A R) for the metric field 
contribution to the axial anomaly, as required on topological grounds that severely 
constrain the form of anomalous terms. It can also be seen without much effort that 
the other topological density in four dimensions, Tr(R A *R), can not possibly con- 
tribute to the answer for it has an excess of anti-symmetric Levi-Civita tensors. 

In the following, we present the explicit computation of the anomaly terms (|2.20|| 



for the case of Lifshitz fermions coupled to background gauge and metric fields and 
show that they are the same, including the overall numerical factors, as for the Dirac 
fermions that are summarized in Appendices B and C. This particular result is not ob- 
vious from the beginning and has important consequences to the integrated form of 
the anomaly. Here, we choose to work with Fujikawa's path integral method for eval- 
uating the anomaly and leave to the interested reader the diagrammatic interpretation 
of the individual terms that are contributing to the final answer. 



2,2 Gauge field contribution to the anomaly 

We will first compute the axial anomaly of Lifshitz fermions coupled to an Abelian or 
non-Abelian gauge field, thus reconfirming the results reported in 1361 . Our presen- 
tation includes all intermediate steps of the calculation and parallels the derivation of 
the axial anomaly for Dirac fermions (see, in particular. Appendix C.l for compari- 
son). According to equation (|2.20|) , we have to compute the regularized sum 



A(t,x) 



lim 

M->oc 



(2.21) 
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using the coupled Dirac-Lifshitz operator ij^D^^ in flat space-time which is mini- 
mally coupled to the external gauge field with = — iA^. The explicit gauge 
field dependence will be extracted using the plane wave basis of solutions of the free 
Dirac-Lifshitz operator. Thus, we use the alternative expression 

A(f,x) = lim Tr / 75 ,-'fc.-"e-('7''iD,)VM''g,7:,x^' (2.22) 

M->oo J {2ny 

and perform the computation acting with the sixth order operator {iy^Dn)-^ on the 
plane waves. The trace is taken over anything available on the right-hand side. 



By definition of the Dirac-Lifshitz operator (|2.3|) , we have the following relation. 



{ij^^f = -D^D?' - l[7^ Y] [ID^, D,] . (2.23) 

The commutator term is not the field strength of the gauge field, as usual, but an 
operator valued quantity with "electric" and "magnetic" components 

[Do, D,] = (fo,D2 + D%, + FoiD^D, + D.D^Fqi + D^Fq/D^ + D^Fq^D,) (2.24) 

and 

Dfc] = (f;)c(D2)2 + {DYFjk + FikD^D^Dj + DjD^&Fik+ 
FjiD^D^D^ + D^DkD^Fji + D^^D^D^ + DjD^FikD^+ 
D%D^D,, + DkD^FjiD^ + DjFjkD^D^ + D^&Fj,,Dj+ 
D^FjiD^D^ + D^D^FjiDk + D^FjiD^D^ + D^DkFjjD^ + 
D^FikD^Dj + DjD%kD^ + 2D%D^^ , (2.25) 

respectively. Here, the derivatives act as operators on anything that appears on the 
right. Thus, many more terms will participate in the evaluation of the axial anomaly 
compared to the relativistic case. 

The action of the interacting Dirac-Lifshitz operator on the plane waves amounts 
to replacing by + ik^ everywhere, so that we may formally write 

when the plane wave factor exp{ikj^x^) passes to the far left and eventually cancels the 
other factor exp(— zfc^x^') appearing in expression (|2.22|) . Here, we set for notational 
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convenience 



Do + iKo = Do + iko , + zK, = ((D, + iki){D + ikf + (D + zl)2(D, + zfc,)) . 

(2.27) 

Obviously this fc-shift produces many terms in the exponential, but the most relevant 
ones are easily selected by rescaling /cq to M'^/cq and A:, to M/c, (due to z = 3 anisotropy) 
and expanding the result around exp[— — {k]k^'f'\ in power series of 1/M. Setting 
= kik\ we have, in particular, 

+ -^b', 7*1 (fitfc* + iFjlhkV) +■■■}, (2.28) 

where • • • denote all other sub leading operator terms that arise from (|2.23|) by expand- 
ing the shifted covariant derivatives; it can be easily seen that these do not contribute 
to the final result. Taking into account the trace identities of products of gamma- 
matrices and, in particular, Tr^'y^'y^'y^'yj'y^) = —Ae^'^'^ for Euclidean space gamma- 
matrices, it turns out that only the quadratic term in the power series expansion of the 
exponential gives a non-vanishing contribution, whereas all other terms either have 
zero trace or they vanish as M — ?■ oo. Thus, we arrive at the expression 

A(f, x) = -2e'''' I ^ e-'^l-'^'TriiFo^k^ + lF,ik,k'){Fj^k^ + ^Fjjc^k^k^)] (2.29) 

that provides a finite contribution to the anomaly when the regulator is finally re- 
moved. 

Next, we perform the Gaussian integration over fco, picking up a factor of y/n, and 
introduce a unit three-momentum vector with components ki (i.e., ki = kki) so that 

A{t,x) = - e^'j' J e-^'k' Tr (fo^f^fc + ^hrFjmhk''' 

+2FoiFjM^ + SFoiFjJ.k^k'^^ . (2.30) 

To complete the calculation we introduce spherical coordinates in the space of three- 
dimensional momenta of length k and note that 

/ e-^'k' = An r dk kh-^' = — Tdx xh-'' = -n'/^ (2.31) 
J Jo 3 Jo 3 



and 



/ d^k e-^k^ kkj = — S,j r dk k^e~^^ = -n^/^j . (2.32) 
J 3 J 9 
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The last identity differs from the first in the angular integration and it can be read- 
ily verified using the orthogonal components of the unit vector ki = sin^sin^, ^2 = 
sin^cos(/) and = cosO in momentum space. Clearly, the last term in expression 
(|2.30|) vanishes, since it is contracted with the totally anti-symmetric symbol e^^^^ (one 
should simply rename the indices i and k), whereas the remaining terms in the trace 
are all proportional to fo/f/fc after performing the integration over k. 

Collecting all terms together, we obtain the final result for the axial anomaly of 
Dirac-Lifshitz fermions in the background of gauge fields (they can be Abelian or 
non-Abelian), which equals 2A{t, x), 

his = - i^^^'^'Tr (Fo.F;ic) = - ^Tr(F A F) . (2.33) 

It is identical to the relativistic case, as advertised before, even though it was obtained 
by combining more terms with a different weight. 

2.3 Metric field contribution to the anomaly 

Next, we obtain the gravitational contribution to the axial anomaly of Lif shitz fermions, 
which is new and constitutes one of our main results. We refer the reader to Appen- 
dices B and C.2 for the notation and comparison with the relativistic case. Throughout 
this subsection, Latin letters a, b, c, d from the beginning of the alphabet are used to de- 
note tangent space-time indices, whereas capital Latin letters I, }, K, L from the middle 
of the alphabet are reserved for the tangent space indices (not to be confused with the 
space indices that are denoted by small case Latin letters i,j,k,l). We also assume 
that the space-time is of the form R x E3 and it comes equipped with a metric G^i/ 
satisfying the projectability condition, as explained before. 

Starting from the general expression (|2.20|) , the primitive form of the anomaly is 
conveniently expressed as 

= lim lim Tr / 75 e-(^7"D„)VM^g,fc,v/V(.y) (2.34) 

using the Dirac-Lifshitz operator i'y^T)^ with = 9^ + (1/8) [70, 7j,]a;^"^ that is 
minimally coupled to geometry via the spin connection. Here, the Dirac gamma- 
matrices are expressed as 7^ = e^ji^a using tangent space-time indices and [7^, 7^,] + = 
25 (fi,. The trace is taken on anything that is available on the right-hand side of the 
equation. We also use the analogue of plane wave basis in curved space that naturally 
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involves the notion of geodesic interval cr{x, x'), as it was originally defined and used 
in the literature | |4Tll42ll43| , in order to extract the background field dependence of the 
anomaly. 

As before, using the definition (|2.3|) of the Dirac-Lifshitz operator, we have the 
following relation, 

ih^'D^)^ = - ^[7^ Y] [D^, D,] . (2.35) 

The commutator terms are operator valued involving the components of the Riemann 
curvature and they turn out to be 

[Do, A] = -1[7„ 7b] [r"\,D^ + D^R"\, + R"\iD'd,+ 

D.D^R'^l + D,R'\jD^ + D^R'\jDi^ (2.36) 

and 

[D^, Dfc] = 1[7„ 7,] [R'^jkiD^f + {D^ fr'jk + R"'ikD'D^Dj+ 

Dp^D^R^'^lj, + R'^jiD'Dj^D^ + D^D,,D^R'^ji + D^R'^ij,D^Dj+ 
Dp^R^^i^D^ + D^R'^jiD^D^ + D^D^R^^jiD^ + DjR'^ikD^D^+ 
D^D^R"^lkDj + D^R^^jiD^D'^ + D^D^R'^^jiDk + DkR"^jiD^D^+ 
D^DkR"^jiD^ + D^R^^j^D^Dj + Dp^R'^i^D^ + 2D^R"^jkD^^ ,(2.37) 

thus providing the metric field analogue of the electric and magnetic components of 
the field strength (|2.24[) and (|2.25|) , respectively. Here, = D/D^ and the derivatives 



act as operators on anything that appears on their right. Putting these together, we 
obtain a relation for (z7^'D^,)^ that generalizes (|B.9|) to the case of Dirac-Lifshitz op- 
erator, but the result is quite lengthy; there is no simple analogue of Lichnerowicz's 
formula for higher order fermion operators. Only those terms that can contribute to 
the axial anomaly will be selected later. 

Acting with the Dirac-Lifshitz operator on exp(zfc^; V^V(x, x')) amounts to replac- 
ing by Dj^i + iApi everywhere with 

(x, x') = fc^ V^, VV(x, x') . (2.38) 

Indeed, one has 
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This relation and its higher derivative generalizations will be used later in explicit 
calculations. For example, a closely related identity that involves the action of two 
covariant derivatives on curved space waves is 

e"'!'^'^^''''^ Q[7«, lb] R'\v + ^V,A, - iVyA^^ . (2.40) 

All terms should be accounted properly before taking the coincidence limit. For now 
it only suffices to note that since Vf'(r(x, x') vanishes in the limit x — > x', we have 

(-^) ' ^ (-^^^^^^) . <^.«) 

after passing exp{ik^iV^cr(x, x')) to the far left and setting for notational convenience 

Do + zKo = Do + zAo, D, + zK, = -- (^{D, + iA,){D + iAf + (D + iAf(D^ + iA^)^ . 

(2.42) 

Then, the primitive form of the anomaly takes the rather simple looking form 

AU ^ V r T r f {hn^v + i'^,)f \ 

A(t,x) = lim lim Tr / - — -r 75 exp — ^ —7 (2.43) 

^ ^ M^oox^x' J (27r)4 -t^ \^ m6 J ^ ^ 

that needs to be evaluated carefully before removing the point split and the cutoff M 
from the integral. 

Note that the operators (|2.42[) resemble (|2.27[) introduced for the computation of the 



axial anomaly in the background of a gauge field, justifying the use of the same nota- 
tion. There is an importance difference, however, apart from the fact that the coupling 
is now taken with respect to the background geometry. The quantity A^(x, x') equals 
kpi in the limit x — ^ x' provided that there are no more derivatives acting on it. The op- 
erators (|2.42|) naturally involve terms with up to two derivatives of A^^ (x, x') that need 
to be extracted before taking the limit x — > x'. There can also be more derivatives act- 
ing on A J, (x, x') when the power series expansion of the exponential is employed for 
the computation of the anomaly. Thus, the Synge-DeWitt tensors | j4Tl 142] . which are 
the multiple covariant derivatives of the geodesic interval cr(x,x') in the coincidence 
limit X — > x', are expected to play essential role in the calculation (see Appendix B.2 
for the mathematical details). As it turns out, Synge-DeWitt tensors with up to four 
derivatives of cr can and will contribute to the calculation of the gravitational anomaly 
of z = 3 Lifshitz fermions, whereas for ordinary Dirac fermions the corresponding 
tensors involve up to two derivatives of a, as in the quantity A^(x, x'). Thus, in the 
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latter case, it is legitimate to replace by k^ from the very beginning and proceed 
with the calculation, as described in Appendix C.2, without worrying much about the 
intricacies of the point split method on curved spaces, which, otherwise, can lead to 
an error. In a nut-shell, one may say that the use of the geodesic interval for evaluat- 
ing the anomaly is rather "cosmetic" for relativistic fermions, serving only the rigorous 
derivation of the final result, whereas for non-relativistic fermions it has very essential 
role. We will say more about this in due course. 

The actual computation of the axial anomaly proceeds in several steps that are 
sketched in the following. First of all , we rescale Aq to M^Aq and Aj to MAf. This is 
not necessarily equivalent to rescaling ko to M^ko and fcf to M^fcf, which is naturally 
implied by the anisotropic scaling of the Lifshitz theory. Note in this respect that if we 
had rescaled the time and space components of both A^ and kfi as just described, the 
rescaled quantities would have been related to each other as 

Ao{x,x') = fcoVoVV(x,x') + -^fcfVoVV(x,x') , (2.44) 
Ai{x,x') = M%ViV°a{x,x')+kjViVja{x,x') . (2.45) 

Although Vo W and V°£r vanish in the limit x — )■ x' (irrespective of the choice of 
shift functions Nf in the ADM decomposition of the four-dimensional metric), their 
multiple covariant derivatives do not vanish in general. The scaling of A^ is the same 
as kji provided that there are no derivatives acting on it prior to the coincidence limit. 
On the other hand, the scaling of the multiple derivatives of A^ do not follow the scal- 
ing of A:^ because there are additional terms with anomalous scaling which are given 
by the appropriate components of the Synge-DeWitt tensors. Such terms can and will 
become relevant in the calculations. Using the notation for the Synge-DeWitt ten- 
sors given in Appendix B.2, i.e., [V^ • • • Vjttr], we have, in particular, the following 
anomalous scaling relations, 

Jim/V^VvAo(x,x')) = fco[V^VvVoVV] + -^fci[V;,V,VoVV], (2.46) 
lim(V„VA-(^.^')) = M%[V.VvV,-VV]+L[V„VvV,-VV] . (2.47) 

Here, we only give the result for terms that involve two derivatives of A^ (x, x') . Terms 
with only one derivative vanish identically in the coincidence limit since the third or- 
der Synge-DeWitt tensors are identically zero, i.e., [V^ VjcC^] = 0. The fourth order 
Synge-DeWitt tensors [V^ V^V^VaC^] are proportional to the Riemann curvature ten- 
sor and their form will be employed later to treat these and other relevant terms. When 
there are more than two derivatives acting on A^ the resulting tensors are non-zero, 
in general, but such terms will not be encountered in the present work. 
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For now, we adopt the rescaling of Aq and A, as a good book keeping device to 
start organizing the power series expansion of the Dirac-Lifshitz operator and the ex- 
ponential of its square and worry later about the possible relevance of the anomalous 
scaling terms when passing to the rescaled momenta Rq and kj. In terms of the rescaled 
quantities, the square of the A-shifted Dirac-Lifshitz operator arising in (|2.43|) , 

- (zy (D^ + iK^)f = (D^ + zK^)(Df' + iK^ + ^[7^ 7l [E*^ + ^'K^, IDi. + i^v] , 

(2.48) 

involves several terms that are conveniently organized in powers of 1 /M as outlined 
below step by step. 

First, we expand the Lifshitz analogue of the A-shifted Bochner Laplacian, setting 

A2 = A,A', 



^(D, + ^K,)(]D^' + m) = -AoAO - + ^3 



(D^ + zK^O(lD^' + = -AqAO - A^ + [2AODo + VqA^' 

+ -^[2A4a"D„ + A^VnA" + 4A2a"A'«V„A^] 

+ -^[SA^D^ + 12A2a"A"^D„D„, + 12A"D„ + B] 

+ irrelevant terms , (2.49) 

where the coefficient functions shown in the third line turn out to be 

A" = a2a^(V"A'" + V^A") + A2A"VmA'" + 2A"A'«A''V«A, , (2.50) 



B = A^A" {5V^An + 5V„ V^A'" + 2V^ V„A'«) + 8A" A'«A''V„ V«A, + 
a2[2(V„A-)(V„,A«) +3(V„A«)(V«A-) +4(V„A^)(V«A-)] + 
A«A'«[12(V„A^)(V,A'-) + 7(V„A,)(V«A'-) + 7(V,A„)(V''A«) + 
10(V„A,)(V''A„,)] . (2.51) 



The first line of equation (|2.49|) contains all terms associated to (Dq + zAo)(D'^ + zA*^) 
apart from DqD*^ which is of order 1/M^ and it is irrelevant. The sub leading terms 
of order 1/ or higher that arise from (D, + zKf)(D' + z'K') are also irrelevant and 
they are omitted from the expansion. It can be seen that they do not contribute to the 
anomaly based on the criterion that wiU be given below. 

Next, we expand the operator [7?^, 7^][D^, + iK^, T)y + zKy] in powers of 1/M 
using the expressions (|2.36|) and (|2.37[) with shifted to + zA^. We obtain, in 
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particular, the following result for the electric commutator after rescaling Aq and A,, 



2M6 



I 



■[7°, 7'][>, 7b] (2R''Va,d'+2R''%,A'd,+2R''%,A,D'+ 



16M5 

R''%,(V.A«) + R'\,{V,^') + i^%(V^A,; 
Y^[7^ 7'][7«, 7b] (k''%.d2 + R%d'd, + R%D,D' 

7'] fiCi + - -^C3 - -^C4 + ^Csl + • • • , (2.52) 



2^' ' ' ' M2 ^ 2M3 2NP 

where • • • include the terms that involve one or two derivatives of the Riemann cur- 
vature tensor, which, in fact, are irrelevant for the axial anomaly. All other terms are 
relevant for the computation and for this reason we give the explicit form of the coef- 
ficient functions Cj, 

Ci = A^Vf Ao + 2A/A" V„Ao , (2.53) 

C2 = 2(A,V"Ao + A"V,Ao)D„ + 2(V„Ao)A"D, + A«(V,V„Ao + V„ViAo) + 

A,V% + (V,Ao)(V„A«) + (V„Ao)(V,A« + V'^A^ , (2.54) 

C3 = 2(V,Ao)d2 + 2(V"Ao)(D,D„ + D„D,)+2(V«V,-Ao + V,V«Ao)D„ + 

2(v2Ao)D, + V^V^Ao + V^Ao + 2A2 VqA, + 4A,A" VoA„ , (2.55) 

C4 = 2(A«VoA, + A,VoA'')D„ + 2(VoA«)A„D, + A,V„VoA« + (VoA0(V„A«) + 
A"(V,VoA„ + V^VoA,) + (VoA")(V,A„ + V„A,) , (2.56) 

C5 = 2(VoA,-)D2 + 2(VoA«)(D,D„ + D„D0+2(V,-VoA" + V"VoA0D„ + 

2(V„ VoA")D,- + V^VoA, + V,V„ VqA" . (2.57) 



It should be emphasized that the results shown above follow from equation (|2.36|) 
by simply replacing ^ + zA^ as well as 



^ba, 7b] ^"\v ^ g[7«, 7b] ^"'liv + ^(V^A, - V,A^) (2.58) 
based on (|2.40|) . If there were no more derivatives acting on it, the coincidence limit 
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of VuAi/ would have been zero, rendering the excess terms in (|2.58|) obsolete. Note, 



however, that the multiple derivatives of no not vanish in the coincidence limit, 
in general, and they should be accounted properly before letting x — ^ x' in all terms 
that originate from the operators DiR^^qi, D^R^^qi, etc. The same remarks apply to the 
terms that arise from equation (|2.37|) and they are discussed next. 

The magnetic commutator term admits the following expansion in powers of 1 / M, 
using the rescaled variables Aq and A;, 

+ 7'] (a^(V^A, - V,A^) +4A,A'a2(V,.A, - V;A^)) 

+ i^bK 7']{(V,Afc - VfcA^)[2A«A2D„ + V„(A«a2)] + (V^-A; - V/A^) x 

[4A"AfcA'D„ + IA^A^D' + lA^A^Dj, + 2V„(A"AfcA') + Vfc(A'A2) + 
V'(AfcA2)] + A"A2v„(VyAfc - VkAj) + 2A"AfcAV„(V;-A, - V/A^-) + 

AfcA2v'(V^A; - ViAj) + A^A^VkiVjAi - ViAj) } 
+ irrelevant terms, (2.59) 

where the irrelevant terms are of order 1 / or higher, as all terms of the same order 
that arise in the power series expansion of (D, + zKj) (D' + iK})/ M^. 

Next, we use the Campbell-Baker-Hausdorff formula for any two operators X and 
y, which is conveniently written here as 



^^"''exp (^x + y + l[x, y] + l[x, [X, y]] - l[y, [x, y]] 



^[y, [X, [X, y]]] + higher commutator terms j (2.60) 



and choose, in particular. 



1 



X = AoAO + A^ y = --^(D,, + zK,,)(Df' + zK?^) . (2.61) 

For practical reasons and in view of the applications that will be discussed in subse- 
quent sections, we consider four-dimensional metrics with vanishing shift functions, 
Ni{t, x) = 0, in the proper time gauge, N{t) = 1, so that we can simply set AqA*^ = Aq 
here and in the following. In any case, (|2.60|) allows to pull out a factor exp(— Aq — A6) 
on the left, which is derivative free, so that we can safely assert irrespective of the 
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remaining operator terms that 



lim e ^0 



A2-A6 



= e 



(2.62) 



after rescaling Rq to M^ko and kj to M/cf and setting kik^ = kp-. Adopting this rescaling 
everywhere, we find that A{t,x) takes the form 



A{t,x) = limM^Tr/-— ^e-'^^'^Ss limexp|x + y + ^[X, y] 





[X, Y]] - ^[Y, [X, Y]] - ^[Y, [X, [X, Y]]] + • • • } (2.63) 



which is the basis for all subsequent manipulations. Note that the successive commu- 
tators of X and Y appearing inside the curly bracket will generate many more terms 
beyond those appearing in X to all orders of 1 / M. We simply have to select the rel- 
evant ones, but they are too many to display them here. These complications do not 
arise in the relativistic case, since the coincidence limit can be taken safely at an early 
stage of the calculation. Even if we had used the Campbell-Baker-Hausdorff formula 
there, the additional terms associated to successive commutators would have been 
irrelevant at the end. 

The rest of the computation proceeds by series expansion of the exponential oper- 
ator in powers of 1 / M, as in a Dyson expansion. The relevant terms for the anomaly 
are selected by the following criteria: 

• they should be of order 1/M^, and 

• they should involve quadratic curvature terms or sufficient number of powers 
of the derivative operator that can yield quadratic curvature terms at the end of 
the calculation, and 

• they should involve sufficient even number of gamma matrices (at least four) so 
that their trace together with 75 can be non-zero, and, finally, 

• they should contain even powers of ko so that their integral with exp(— Aiq) can 
also be non-zero (recall that ko runs from —00 to +00). 

Note at this end that all other terms of order 1 / that do not involve four derivatives 
will not be capable to produce a topological density in four space-time dimensions; 
it can also be verified by direct computation that all such terms cancel against each 
other. Terms of the higher order obviously give zero in the limit M — )■ 00, whereas 
terms of lower order ought to cancel for, otherwise, the anomaly would be infinity; it 
can also be verified by direct calculations such terms indeed cancel order by order. 
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The first three criteria are obvious and they are similar to the ones used for the eval- 
uation of the gravitational contribution to the axial anomaly of relativistic fermions, 
apart from the fact that the relevant terms there are of order 1 / rather than 1 / by 
the difference in time scaling. Then, because of this difference, and the higher order 
structure of the Dirac-Lifshitz operator, the terms contributing to the anomaly turn 
out to be very different from those arising in the relativistic case. For example, as will 
be seen later, terms that involve products of up to eight gamma matrices with 75 con- 
tribute to the final answer, whereas in the relativistic case only the trace of 75 with four 
gamma matrices comes into play. Most importantly, higher order derivatives of the 
geodesic interval are generated (either in X, as summarized by adding together the 
individual contributions (|2.49|l , (|2.52|) and (|2.59|) , or in the successive commutators of 



X with Y or in the subsequent expansion of the exponential operator shown in (|2.63|) ) 
and they all contribute to the axial anomaly of Lifshitz fermions in various ways. 

The fourth criterion looks rather superfluous at first sight, but, in fact, it helps to 
select some additional terms that could have been easily missed otherwise. Recall at 
this point the mismatch in the scaling of Aq and A; as compared to Rq and fc, when 
multiple derivatives act on A^(x, x') (see, in this respect, equation (|2.47|l ). Certain 
terms that are seemingly of order 1 / M^, like the cross product of A^Dq/M^ appearing 



in the first line of (I2.49D and R^^/VjAVM^ appearing in the third line of (I2.52D , yield 



the following contribution in the coincidence limit 

^limR^\iAoVoVA' = ^r%,fc2[VoV,V'vV] + O , (2.64) 

which is non-zero. Such terms are effectively of order 1 / M^, they contain four deriva- 
tives, they have sufficient number of gamma matrices ([7*^, 7'] [ja, 7b] in the particular 
example we are discussing here, following from (|2.52[) ), and, most importantly, they 
contain an even number of powers of ko (Rq in the present example) so that the inte- 
gral over Rq is non-zero. One has to extract carefully all such terms from the power 
series expansion of the exponential, although they are more rare compared to the other 
terms for which the mismatch in the scaling is irrelevant. Finally, the fourth criterion 
is also used to eliminate all cross terms of A'^Dq/M^ with the 0{1/M'^) terms arising 
in the expansions (|2.49|) or (|2.59|) and alike, because they contain odd powers of ko in 



the coincidence limit. This is the reason why such higher order terms were omitted 



from the expansions (|2.49|) and (|2.59|l as irrelevant. 

The various terms that meet all four criteria are still too many to display them one 
by one. Here, we will only sketch the steps we have taken for their determination by 
grouping them together in different classes. The intermediate details are straightfor- 
ward but very cumbersome and they are left to the interested reader to complete. 
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In the course of the calculation we also need the following integrals over the space 
of three-momenta (whereas the integration over Rq is very simple). 



_ (2n-l)!! 



In ^ I d'k e-^"k^- = n''\ (2.65) 



as well as the curved space analogue of equation (|2.32[) , which now reads 



dhe-^'k^''kjkj=^-^5jj. (2.66) 
3 

The proof is easily done by introducing tangent space indices I, } (to be distinguished 
from the spatial indices i, j on the slices E3 of R x E3), so that the unit vectors fc, that 
are defined by ki = kkj are decomposed with respect to the dreibeins associated to 
the metric G/y = 5jje\e^j as follows, ki = e^ikj (see also Appendix C.2 for the corre- 
sponding integral in the space of four-momenta). The rest proceeds as in flat space, 
choosing, in particular, the components of kj as ki = sin^sin^, ^2 = sin^cos^ and 
^3 = cos9 in terms of spherical coordinates in momentum space. 

Likewise, we evaluate the following integrals that will all be needed in the course 
of the calculation (in fact, up to n = 3), 



and 



J e-''\'" UMl = ^ i^iM + ^iK^jL + SilSjk) (2.67) 

/J 
dhe-^k^"" kjkjkKkikMkN = ^ij {^kl^mn + ^km^ln + ^kn^lm) 

+ ^IK i^JL^MN + ^JM^LN + ^JN^Lm) 
+ ^IL i^JK^MN + ^JM^KN + ^JN^Km) 

+ Sim {Sjk^ln + ^jl^kn + ^jn^kl) 

+ Sin {SjkSlm + Sjl^km + SjmSkl) ■ (2.68) 

After these explanations, we are now in position to describe the structure of the 
various terms that contribute to the gravitational anomaly of Lifshitz fermions in the 
limit M — ?► 00 and x ^ x' . We will also provide some identities that are necessary to 
cast them in familiar form. It is most convenient to split the terms in two groups as 

A{t,x) = Ai{t,x) + A2{t,x) , (2.69) 

where Ai{t,x) denotes collectively the terms that contain no derivatives of and 
A2{t, x) denotes all other terms that contain two derivatives of and which can be 
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evaluated using the fourth-order Synge-DeWitt tensor. As noted before, terms that 
contain only one derivative of A^; vanish in the coincidence limit, whereas higher order 
derivative terms do not enter in the calculation. 

(i). Ai terms: These terms are the easiest to describe since the coincidence limit can 
be taken from the beginning. There are two different type of such terms, so we write 

(1) (2) 

Ai{t,x) = A\ '{t,x) + A\ '{t,x), which can be read directly from the power series 
expansion of exp(X + Y) (note in this respect that [X, , Y] and all other higher order 
commutators of X and Y can only produce additional terms that contain derivatives 
of prior to the coincidence limit). First, we display 

A^i\t,x) = ^^^Tr J dhe-^'js [7°, 7'] [la, lb] x 

+ 3k^ (^R^^i + 2R"\ikik^^ (^D^ + AkjkoiD^^ 
- 18 kjkk^^ I^R'^, + lR'\ikik^^ D^D^ 

-R''\,D^-R'\,[d\d,] + ], (2.70) 



which contains the cross product terms between (|2.49|) and (|2.52[) excluding all deriva- 
tives of . In particular, the first line in (|2.70|) includes the product of terms of order 



1/M with terms of order 1/MP , the second line the product of terms of order 
with 1 / M^, the third line the product of terms of order 1 / M-squared with terms of 
order 1 / M^, and, finally, the fourth line includes the terms of order 1 / appearing 
in (|2.52|) . In writing (|2.70[) , we introduced unit vectors ki and integrated over /cq, thus 
picking up a factor of tt^^^. Then, we can easily perform the integration over three- 
momenta based on the relations (|2.65|) - (|2.67|) , which are used here up to n = 2, and 



find that the resulting terms are either of the form J^^^^D^ or R^^QiiDi, D']+. These 
terms cancel against each other separately before taking the trace of the gamma ma- 
trices, and, therefore, we obtain 

Af^(i,x)=0. (2.71) 

Next, we display the remaining terms of this group that contain the cross product 
terms of order 1/M* and 1/M^ found in the expansions (|2.52[) and (|2.59|) . They yield 



Af^(i,x) = Y^T^Y^^Tr j d\e-''k' 75 [l^ Y] [la, It] \lK 7^] \lc, Id] x 

jflfe T)cd I A-nah T)cd t.tm i ouflfc lycd t tl i Qr?ab -ncd tTJt.tm 



R'orR'jk + ^R'^rR'^jJ^iT" + lR'\iR''jkKk^ + SR'\iR'^ ^Ji,k%V^ (2.72) 
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after introducing the unit vectors fc, and integrating over kQ. These terms are anal- 
ogous to the group of terms (|2.30|) that contribute to the axial anomaly of a Lifshitz 



fermion coupled to gauge fields, but they are more complicated now because they 
involve the product of 75 with eight gamma matrices. As we will see shortly, they 
provide a non-vanishing contribution to the gravitational form of the anomaly, which 
has no analogue in the relativistic case. Integrating over the three-momenta, based on 
relations (|2.65|l - (|2.67|) , which are used again up to n = 2, we arrive at the expression 



^i'^^^'^^ = 45 ■ 512712 Tr (75 l'']\[la, Y]\bc, 7d]) X 

(53R"\,R''j, + SR'\^R'' j, + 80,^0""^%,^'^%) (2.73) 

which is written here using the inverse vierbeins E^/^' that trade the space-time gamma- 
matrices 7?' {}i = 0, k) with their tangent space-time counterparts 7"'. 

At this point we employ the trace identities found in Appendix A in order to evalu- 
ate A\ {t, x). Using the trace of 75 with the product of eight gamma-matrices, which 
is given by equation (|A.10|) in terms of flat space-time indices and which should be 
appropriately adapted to current notation, we obtain the resul1§ 

Af ) it,x) = {RomR'^jk - Rjm^'^rn) ■ (2.74) 

A useful identity that can be used here and later to combine terms of the form Rg/Oi^'^'/A: 
and Rj^QiR^^in (these are the only terms that can arise in the computation of the gravi- 
tational anomaly in non-relativistic theories) is provided by 

le'^'^^RamR'' = e''^' {RomR'^jk - RjmR'\n) ■ (2.75) 

This look rather odd at first sight, but it is in fact true in general for all four-dimensional 
geometries. The proof of this identity follows by brute force, writing down all terms 
that arise from summation over repeated indices. Thus, overall, we obtain 

Ai{t,x) = ^^-^e'^j^RMrR'^jk , (2.76) 
which has the right form but comes with positive sign. 



^In the course of this calculation we also find terms of the form e'^'^*^Ro/;l^'""mn/ but they vanish 
because Koiji: + ^o/fci + ^Oici) = 0. Such terms also arise in the computation of A2. We also employ the 
identity 



in 



2 

to simplify the intermediate expressions and express everything in terms of i^o/O/^'^'y/c srid Rji(QiR^"i 
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(ii). Aj_ terms: This group contains all other terms that meet the four criteria out- 
lined before and they have derivatives of = kySJ ^^SJ'^ cr{x,x'^. Such terms will 
arise from the power series expansion of the exponential, which now apart from 
exp(X + y) it also receives contributions from the commutators [X, Y], [X, [X, Y]], 
[Y, [X, Y]] and [Y, [X, [X, Y]]]. The relevant terms are too many to write them 

down explicitly We only mention here that they can also be divided in two sub- 

f 1) 

groups: A^^ it, x) that encompasses the terms for which the scaling of Aq and A; is 
same as fco ^rid k[ and Jvr^ {t,x) that encompasses the terms for which the scaling re- 
lation is anomalous (their naive order with respect to the scaling of A^ is 1 / and 
respectively, as explained earlier). In the first subgroup we may simply replace 
V^,Vt;Ao(x,x') by fco[V^,Vi/Vo VV] and V^,Vt/A,(x, x') by fcy[V^Vi,ViVV] in the co- 
incidence limit X x', whereas in the second subgroup the anomalous scaling term 
M^ko[VpVvViV°cr] shown in (|2.47|) should be used instead to replace V^VvA;(x,x'). 



In either case, we employ expression (|B.18|) for the fourth order Synge-DeWitt tensor. 



[V^ Vi^VkVa^t] = ^ (K^,Av + R^Xkv) , (2.77) 
and obtain after a very long calculation that also involves the term-by-term integra- 



tion in momentum space (using all relations (|2.65|) - (|2.68|) up to n = 3) the following 
result 

19 



1 9 

= — ^e°'^■*^R«M^'^•fc • (2.78) 

2880712 "^"^ ^ ' 



Note that the contribution of A2{t, x) comes with the opposite sign compared to Ai{t, x) 

Combining these expressions, it turns out that all terms that contribute to the sum 
A{t,x) = Ai{t,x) +A2{t,x) yield a net result 

Ait,x) = . (2.79) 

Thus, in conclusion, the axial anomaly of a Lifshitz fermion in the background of a 
metric field, which equals 2A(t, x), turns out to be identical to that of a Dirac fermion, 
i.e., 

V,/f = - ^e''i'R''\K, ^ - Y^Tr (R A R) (2.80) 
with the same overall numerical coefficient, as advertised before. 

Note that if we were computing the anomaly of a Weyl-Lifshitz fermion, as op- 
posed to the Dirac-Lifshitz fermion that we have been considering so far, the overall 
coefficient would be half of (|2.80)) . 
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Finally, we end this section with some general remarks regarding alternative ap- 
proaches to the axial anomaly which can also be used in non-relativistic field theories. 

Although the Fujikawa method for computing the anomaly is straightforward to 
implement in principle, it turns out to be very cumbersome for higher order opera- 
tors, as we have just seen. Later, in section 4, we will provide an alternative method 
for determining the overall coefficient of the topological density Tr(R A R) based on 
the computation of the ^/-invariant of the three dimensional Lifshitz operator on E3 
and obtain the corresponding index theorem associated to the integrated form of the 
gravitational anomaly for the (3 + 1) -dimensional Dirac-Lifshitz operator. In effect, 
this will provide a faster though more mathematically oriented way to arrive at the 
same result. 

There is also an alternative physical method based on supersymmetric quantum 
mechanics to compute efficiently and reliably the axial anomaly, as described in the 
classic work IITOl for the cases of the Dirac and the Rarita-Schwinger fermion opera- 
tors. In that context, one has to find the appropriate supersymmetric two-dimensional 
non-linear sigma and reduce it to + 1 dimensions so that the index of the fermion 
operator can be obtained by path integral methods of the underlying supersymmet- 
ric quantum mechanics model. The general philosophy of this procedure is to find a 
one-dimensional quantum mechanical system defined on the Riemannian space-time 
manifold M4 such that its Hamiltonian is the square of the fermion operator. Then, 
the expression for A{t, x) is interpreted as the partition function for an ensemble with 
density matrix |0 = 75 exp[— /3(z7?^D^)^] at temperature j6~^ (/3~^ = for relativistic 
Dirac fermions) and its evaluation is equivalent to the high temperature expansion 
of the system under periodic boundary conditions for both bosons and fermions. The 
details of the calculation are usually much simpler in that approach than in Fujikawa's 
treatment of the same problem as they circumvent the very many terms that can arise 
otherwise. This approach also allows to relate the result to the computation of the 
fermion propagator in a constant uniform magnetic field [i46il (but see also pTi], in 
particular p. 100), thus providing a more physical interpretation to the origin of the 
topological density that enters into the anomaly (in effect, Schwinger's computation 
provides a physical derivation of the so called A-roof genus of the fermion operator on 
M4). We have not considered this possibility for Lifshitz fermions so far, but we think 
that the right framework should be provided by sigma models with non-relativistic 
supersymmetry that generalize the construction reported recently in [48] to theories 
with multi-component scalar fields associated to a z = 3 Lifshitz non-linear sigma 
model in 1 + 1 dimensions. If this step is successfully implemented in the future, the 
treatment of anomalies in quantum field theories with non-relativistic fermions will 
be on par with the relativistic theories. Work in this direction is in progress and will 
be reported elsewhere. 
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3 Instantons of Hofava-Lif shitz gravity 



In this section we present an overview of the instanton solutions of z = 3 Horava- 
Lifshitz gravity following our earlier work on the subject ||38| . After some general 
considerations that set up the notation and framework for the definition of instantons 
in terms of geometric flows, we focus on the explicit construction and classification 
of solutions with SU{2) isometry. These will be used later as metric backgrounds 
to evaluate the integrated form of the axial anomaly and provide examples of chiral 
symmetry breaking in Lifshitz theories induced by gravity It will turn out that such 
phenomena become possible only in a unimodular phase of Horava-Lifshitz gravity, 
which is thought to govern the deep ultra-violet regime of the theory. Our treatment 
of instantons is rather complementary to the results presented in f38l giving a more 
qualitative picture that is also appropriate for comparison with similar solutions that 
arise in Euclidean Einstein gravity. The presentation is self-contained as several of 
these results are not widely known. We also reveal the curvature characteristics of the 
solutions that make them relevant for applications to chiral symmetry breaking. 

3.1 General considerations 

Horava-Lifshitz gravity is a non-relativistic theory that is defined by assuming that 
space-time is of the form IR x E3. It is naturally formulated using the ADM decompo- 
sition of the space-time metric G^v, 

ds^ = -N^dt^ + gij (^dx' + N'dt^ {dx^ + N^di^ (3.1) 

with lapse and shift functions N and Nj, respectively. The spatial slices E3, which 
are the leaves of a foliation, have induced metric gij and extrinsic curvature (second 
fundamental form) Kij, 

% = ^ i^^S^j - V.N^- - VjN,) , (3.2) 

which provides the conjugate momentum to the metric gij in the canonical formalism. 
The momenta are simply expressed using the superspace metric as 

K'i =\^d^g'j^^Ki,i. (3.3) 

Recall that superspace is defined as the infinite dimensional space of all Rieman- 
nian metrics on E3, and, as such, it provides the arena for geometrodynamics. It is 
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naturally endowed with a metric 

Q'''' = l{gV + g'g'')-^g''g'' (3.4) 



that generalizes the standard DeWitt metric 11491 using an arbitrary parameter A. There 
is also the inverse metric 

1 A 

Gijkl = 2 {gikgjl + gilgjk) - ^j^^gijgkl (3-5) 

which satisfies the relation 

G'^^'Gklmn = \ (/Jn + • (3-6) 

In general relativity A = 1, but in Hofava-Lifshitz gravity A can take more general 
values. It is important to note that the superspace metric is positive definite when 
A < 1/3 and it is indefinite when A > 1/3; in the latter case, which includes the 
value A = 1, the indefiniteness of the DeWitt metric accounts for the conformal factor 
problem in all gravitational theories. The value A = 1/3 is special since the DeWitt 
metric becomes degenerate. It is closely related to the limit A — > ±oo for which the 
inverse DeWitt metric becomes degenerate. In these special cases, the conformal factor 
of the three-dimensional geometries decouples and a unimodular phase of the theory 
emerges; we will say more about this later as it provides a very important class of 
models for the purposes of the present work. 

With these definitions in mind, the action of Hofava-Lifshitz gravity takes the fol- 
lowing form written as a sum of kinetic and potential terms, p0ll3T|| , 

S = -| y did^x^di^ N K^jO'l'^'K,,, - y y dtd^x,/Mg N WQiju^^^ , (3.7) 

where = 32nG is the four-dimensional gravitational coupling and is the gradi- 
ent of a suitably chosen local functional W[g], 

E'> = _ 1 . ,3.8) 

Here and in the following we are only considering the form of the theory which is 
said to satisfy the detailed balance condition, meaning, in particular, that the potential 
term is derived for a superpotential functional W[g]. Other forms of the theory have 
also been studied in the literature, but they will not be in focus in the present work. 
Ordinary general relativity resembles the form of the action (|3.7|) but the potential is 
given by minus the Ricci scalar curvature R of the metric on E3. In that case, one 
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may still derive the potential from a superpotential functional, which is known as the 
Hamilton-Jacobi functional of general relativity, but it is a non-local functional^ of the 
metric g in general. 

The resulting theory exhibits scaling anisotropy in space and time depending on 
the choice of W[g]. We are going to consider models with anisotropic scaling expo- 
nent z = 3 by choosing as superpotential functional the action of three-dimensional 
topologically massive gravity on E3 which is defined as follows, Il50l , 

^[g] = \ I d^x./M^ (R - 2Aw) + -Wcs . (3.9) 

The first term is the usual Einstein-Hilbert action in three dimensions with gravita- 
tional coupling Kw and cosmological constant A^, whereas the second term is the 
gravitational Chem-Simons action [511, which is written in terms of the Christoffel 
symbols of the metric g, 

= Tr (^cv Adcv +'^cv Acv Acv^ (3.10) 

using the fully antisymmetric symbol in three dimensions with e^'^^ = 1 or equiva- 
lently using the spin connection one-forms co (they should not be confused with the 
coupling constant co of the Chern-Simons term). Obviously, this theory is not invari- 
ant under parity, since orientation reversing transformations on E3 flip the sign of the 
Chern-Simons coupling co. This affects the vacuum structure of topologically massive 
gravity, since it introduces chirality in the geometry, and it also has implications for 
the instanton solutions of Horava-Lifshitz gravity. 

The classical equations of motion that follow by varying W[g] take the form 

^^7 - l^S^i + ^^Stj - —Qj = , (3.11) 

Z CO 

where C/y is the Cotton tensor of the metric g, which is defined as follows. 

The Cotton tensor is a covariantly conserved and traceless tensor, which is third order 
with respect to the space derivatives of the metric, and it accounts for the anisotropic 

^Non-local functionals can also be used formally to recast Hof ava-Lifshitz theory without detailed 
balance into detailed balance form, but these generalizations will not be pursued further. 
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scaling exponent z = 3 of the associated Hof ava-Lif shitz gravity. A special case arises 
when Kw becomes infinite and the Einstein-Hilbert term drops from the action W; 
the resulting three-dimensional theory is conformal gravity whose vacua are confor- 
mally flat metrics satisfying the field equations Qj = (see, for instance, Il52l for more 
details). Another special case arises when co becomes infinite and the gravitational 
Chern-Simons term drops out from the action W; the resulting three-dimensional 
theory is ordinary Einstein gravity and in this case the anisotropic scaling of the as- 
sociated Horava-Lif shitz gravity reduces to z = 2. In all cases, the vacua of three- 
dimensional gravity provide static solutions of the 3 + 1 non-relativistic theory (|3.7)) . 

Hofava-Lifshitz gravity, in its original formulation that we are using here, does 
not exhibit general coordinate invariance but only a restricted symmetry associated 
to foliation preserving diffeomorphisms of space-time. Then, for general values of A, 
the physical viability of the theory becomes questionable because it contains a spin-0 
graviton mode that is unwanted in the infrared limit and it should decouple. For cer- 
tain values of A this problem does not arise; we will say more about this later, since 
most of our present work will concentrate to those special values of A. For now it 
suffices to say that we are adopting the projectable version of the theory, and, further- 
more, in view of the applications, we choose the lapse and shift functions 

N(i) = l, Ni{t,x)=0 (3.13) 

exactly as in the case of Lif shitz fermions, unless stated otherwise. Setting N{t) = 1 
amounts to using proper time. 

Gravitational instanton solutions can be obtained by considering the analytic con- 
tinuation of Horava-Lif shitz theory in the time coordinate ||38l. The four-dimensional 
Euclidean action that follows from (|3.7|) by inverting the potential term is bounded 
from below provided that the superspace metric is positive definite. Thus, assuming 
that A < 1/3, we have the following relations 

SeucI. = ^/ dtd^xy^K,jG'i^^Ku + '^ J dtd^xy^ E^igij,,,E^^ 

+2 J dtd^Xy^MgK.jE'i . (3.14) 

We also implicitly assume that E3 is a compact manifold with no boundary in order 
to avoid unnecessary complications that may arise from additional boundary terms 
in the Euclidean action (persistent boundary effects in gravitational theories require 
more careful treatment as in the monopole sector of gauge theories). It is then clear 
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that the action is bounded from below as 

/f \ f t^W 

dtd\^/d^K,jE'j = TJ dtd^xy/d^ E'jdtgij = ±2 y • ^^-^^^ 

The lower bound is saturated for special configurations satisfying the first order equa- 
tions in time 

which describe gradient flow equations for the metric on Z3. They provide extrema of 
the four-dimensional action, and, as such, they also satisfy the classical equations of 
motion which are second order in time. 

The instantons of Horava-Lif shitz gravity are defined to be eternal solutions of the 
geometric flow equations (|3.16|) that exist for all time — 00 < ^ < +00. The ± sign 



refers to instantons and anti-instantons, which are mutually related hy t ^ —t. These 
flow lines interpolate smoothly between any two fixed points of (|3.16|) , i.e., between 
different extrema of W[g], and they have finite Euclidean action 



_ 1 

'instanton — 2 



|AW| (3.17) 



given by the difference of W at the two fixed points reached as i — ?■ =Foo. It can be easily 
seen that W changes monotonically along the flow lines, ±dW / dt > 0, provided that 
A < 1/3 so that the superspace metric is positive definite. Therefore, the instanton 



action is always positive definite, justifying the use of the absolute value in (|3.17|) . 
This definition is analogous to the instanton solutions of point particle systems that 
interpolate smoothly between different degenerate vacua (minima of the potential). 
In the present context, the absence of singularities along the eternal solutions of the 
flow equation (|3.16|) ensures that the corresponding space-time metric is complete 



and regular. Other flow lines that may only exist for a finite time interval before they 
become extinct are not physically acceptable as they are inflicted with singularities. 

Specializing W[g] to the action of three-dimensional topologically massive gravity, 
we obtain from (|3.16|) the following third order non-linear equation for g, 

3.x, = [r, - ^(3X31) - 31^*.) ± TJC, , (3.18) 

which is called Ricci-Cotton flow [38J. Its eternal solutions are the instantons of z = 3 
Hof ava-Lifshitz gravity that will be studied later in detail focusing, in particular, to a 
mini-superspace model that allows us to derive explicit results and obtain complete 
classification of all such special configurations. 

Although the flow lines of (|3.18|) depend explicitly on A, the fixed points are in fact 
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independent of it (they are vacua of topologically massive gravity irrespective of A). 
As such, they satisfy equation (|3.11|) , which is equivalently written as 



^tj-lR8tj-—Cti = and R = 6A„ . (3.19) 



The first relation is the traceless part of (|3.11|) and the other is the trace telling us. 



in particular, that the two end-points of the instantons are metrics on Z3 with Ricci 
scalar curvature R fixed by Aw The curvature of E3 changes along the flow lines. The 
flow equation becomes independent of Aw only when A is infinite, in which case the 
end-points of the instantons can have arbitrary R (in this case, the fixed points satisfy 
only the traceless part of (|3.11|) ). We also note for completeness that the Ricci-Cotton 



flow (|3.18|) becomes the celebrated Ricci flow when oj — )■ ±00. Actually, in that case, 
we obtain Hamilton's Ricci flow for A = 1/2 (see, for instance, the collected papers 
on the subject ||53l and the textbooks ||54ll55| ). Finally, when Kw — > 00, equation (|3.18|) 



specializes to the Cotton flow that was introduced in the literature in more recent years 

nsa. 



In all cases, the instantons of Horava-Lifshitz gravity are regular space-time con- 
figurations of the form R x Z3 that extend for all time — 00 < i < +00 and they are 
suspended from two fixed points. They can be visualized as infinitely long cylinders 
which provide the envelope of the geometric flow on Z3 as shown in Fig.l. The slices 
(which are depicted here by suppressing two of the three spatial dimensions) are the 
portraits of the geometry shown at different instances of proper time. Similar pictures 
can also be drawn for the instantons of Einstein gravity in proper time, but the cylin- 
der is semi-infinite in that case, since the solutions are suspended from removable 
singularities (nuts or bolts) at one end. We will say more about these differences in 
due course. 



Figure 1: Schematic picture of an interpolating tnstanton solution IR x E3 
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3.2 Unimodular phase at special values of A 



Next, we examine the properties of Horava-Lif shitz gravity at special values of A that 
will play important role in the applications. First, we consider the case A = 1/3 which 
amounts to 

G'^^'gM = , n^^ = ^ Vd^ G'^"K,i = I [k^^ - Ig'^K^ ■ (3.20) 

The first relation states that the superspace metric is degenerate as it develops a null 
eigen-vector provided by g, whereas the second relation shows that the conjugate 
momentum becomes traceless. More generally, projects any tensor to its traceless 
part. At the same time, the inverse metric in superspace is not well defined (it becomes 
infinite), and appropriate projection is required to make sense of the theory at A = 1/3. 
Thus, at A = 1/3, we are led to define 

= \ ( A^'' + ^'g^^) - \gH' (3.21) 

and 

1 1 

Qijkl = 2 {gikgjl + gilgik) - ^gijgkl / (3.22) 

so that Qijj^i also projects any tensor to its traceless part. In this case we have 

G'^^'Gklmn = \ {5[A + - \g'^gmn , (3.23) 

which differs from the corresponding relation (|3.6|) by subtracting its trace part. This 
modification is necessary for consistency of the projection. 

The projection to traceless tensors is clearly related to the role that the conformal 
factor (volume of space) has in the theory. The conformal factor of the metric g is 
a propagating degree of freedom with positive kinetic energy in superspace when 
A < 1/3 and with negative kinetic energy when A > 1/3. At the boarder value 
A = 1 /3 its kinetic energy vanishes, and, therefore, the conformal factor has no dy- 
namics at all. This does not necessarily imply that the theory is invariant under rescal- 
ing, but it means that the volume of space can only appear as spectator. Then, the 
action of Horava-Lif shitz gravity (|3.7[) at A = 1/3 takes the form (here we reinstate 
momentarily N and N, in order to make some general comments) 



X 



1 1„ \ 1 



w 



1 1„ ,A 1 



( ^^7 - s^S^j - -Qj :^{R'^- -M^ ] - -P ■ (3.24) 



K^«7 V 3 



CO 
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The underlying three-dimensional theory of topologically massive gravity becomes 
unimodular in this case, since the volume of space E3 is held fixed. Its field equations 
are solely described by the traceless condition 

^'7 - l^Srj - ^Qi = (3.25) 

without imposing any restriction on R. Then, Aw can assume arbitrary values as it 
has the interpretation of an integration constant arising from the contracted Bianchi 
identity (this is the standard viewpoint for the cosmological constant in all unimod- 
ular theories of gravity). Similar reasoning applies to the limiting case X ^ ±00 that 
also leads to the action (|3.24[) . 



The kinetic term of the action (|3.24[) is invariant under the following transformation 

gij n2(i, x)gjj , N n.^{t, x)N , Ni n^{t, x)Ni (3.26) 

that changes the measure of integration as ^/detg N a/ detg N. If W[g] were 

the action of three-dimensional conformal gravity, by letting Kw — > oo, the potential 
term of Horava-Lifshitz gravity would also remain invariant under this transforma- 
tion, since the Cotton tensor transforms as C/y Ci~^Cij. In that case, the theory 
would exhibit anisotropic Weyl invariance, as it was first noted in IISOllSTI , thus turn- 
ing the rigid anisotropic scaling of space and time coordinates {t,x) — > {a^t,ax) into 
a larger local symmetr}]^ However, the deformation that arises when is held fi- 
nite breaks anisotropic Weyl invariance because the traceless part of the Ricci tensor is 
invariant under conformal transformations and the potential term of the action (|3.24|) 



depends explicitly upon O. Thus, at A = 1/3 (and likewise at A = ztcxj) the action of 
Horava-Lifshitz gravity is not Weyl invariant, in general. It is only invariant under 
the restricted group of foliation preserving diffeomorphisms that keep the volume of 
E3 unchanged. They account for the additional Hamiltonian constraint = 0, which 
now follows from the unimodular projection of topologically massive gravity on Z3. 

It is also interesting to note that the scalar (spin-0) graviton mode which haunts 
Horava-Lifshitz gravity for generic values of A decouples in this case and only the 
usual tensor graviton mode remains in the physical spectrum. This was first noted 
in the literature by considering linear perturbations around the Minkowski vacuum 
of the anisotropic Weyl invariant theory of gravity Il30i l3ll , but the proof was later 
generalized to the more general theory defined by (|3.24|) ||57| . Thus, the phase of the 



theory arising at A = 1 /3 is not plagued with the inconsistencies that otherwise haunt 



''Only in that case it is appropriate to consider the non-projectable version of the theory, since the 
lapse function necessarily depends on all space-time coordinates; even if we were choosing N to de- 
pend only on t, the transformation | |3.26|| would make it to depend on both t and x. 
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the non-relativistic theories of gravity. We will stick with it in the applications that 
will be considered later for reason that will become transparent in due course. Two 
important questions will remain unanswered though, but they are both lying beyond 
the scope of the present work: one is the extrapolation to large scale physics and the 
other is the fate of Weyl invariance and unimodular symmetry upon quantization. 

The instantons of the gravitational theory (|3.24|| arising at these special values of A 
are eternal solutions of the flow equation 



3.S., = T^(R,;-3%;j±-C,, (3.27) 



which can be obtained formally from (|3.18|) by taking the limit A — > ±00 (we opt the 



limit A — > —00 as it fits in the range A < 1/3 that was considered earlier). We call 



equation (|3.27[) normalized Ricci-Cotton flow because the driving terms of the metric 
deformations are traceless tensors that preserve the volume of space, as consequence 
of unimodularity in three dimensions. The fixed points of the flow lines satisfy equa- 
tion (|3.25|) and they depend on the volume of space, which, however, is appearing now 
as a spectator (as noted before, the traceless Ricci and Cotton tensors scale differently 
under conformal transformations of the metric). As a result, the end-points of these 
instantons can have different Ricci scalar curvature. This should be contrasted to the 
curvature of the fixed points of the unnormalized Ricci-Cotton flow (|3.18|) , which is 
the same. 

Next, we obtain explicit solutions by considering instantons with sufficiently large 
group of isometries so that the equations become manageable. Instantons with less or 
no isometries may also exist but we are lacking the mathematical tools to investigate 
them properly. A few general remarks about this will only be made at the end of 
section 3.4. 



3.3 Instantons with S 11(2) x U{1) isometry 

Instantons with 51/(2) isometry arise by choosing E3 ~ and introducing homoge- 
neous metrics of Bianchi IX type 

ds^ = ^r{t){a'f + ii{t){a^f + 73(i)(^')' (3-28) 

in terms of the left-invariant one-forms of the group Sll{2) that satisfy the defining 
relations 

da^ + ^e^]K<^^ A (7^ = . (3.29) 
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Remarkably, this ansatz leads to consistent reduction of the Ricci-Cotton flow to a 
system of ordinary differential equations for the metric coefficients 7; which are taken 
to depend only on time. The same applies to all other homogeneous model geome- 
tries on E3, but we are only considering here the mini-superspace model of the flow 
equations with SU{2) symmetry. The resulting equations are highly non-linear and 
unfortunately they can not be integrated when the geometry of is totally anisotropic 
with 7i 7^ 72 7^ 73- It is possible to investigate these equations numerically and study 
some of their qualitative features, which turn out to be sufficient for establishing the 
existence of SU{2) gravitational instantons, in general, following 1381 . We will not re- 
peat the same line of presentation here, but rather give a complementary account of 
the solutions and then expand on their properties that are most relevant to the present 
work. 

In view of the applications that will be considered in the next section, it suffices to 
restrict attention to axially symmetric configurations (also known as Berger spheres) 
by letting 

71 = 72 = ^ , 73 = ^ (3-30) 

and present all gravitational instantons with enhanced isometry SU{2) x U{1). Here, 
L is a characteristic length scale so that the volume of space equals In^L^. We also 
confine ourselves to the normalized Ricci-Cotton flow, in which case the volume of 
space is preserved by the evolution (a few remarks about the unnormalized flow as- 
sociated to more general values of A will be made later). Then, in this case, the system 



(|3.27[) reduces consistently to a single differential equation for x{t) in proper time f. 



dx , 4k^ / I .\ / 1 coL 



whose solutions will provide instantons in the unimodular phase of Hofava-Lifshitz 
gravity. The variable x is positive definite and it parametrizes the shape modulus of 
space. The value x = 1 corresponds to a round sphere, where x > 1 (respectively 
X < 1) describes a geometrically deformed sphere which is squashed (respectively 
stretched) along the axis of symmetry. Of course, the choice of axis is arbitrary because 
one can permute the three principal directions of S^, and, as a result, our construction 
of instantons is unique up to a Z3 permutation symmetry. 

The instanton solutions are by definition eternal solutions of the flow equation 
that exist for all time —00 < t < +00 and interpolate smoothly between any two 
fixed points, provided that there are more than one degenerate vacua in the system. 
It is clear, by inspecting equation (|3.31|) , that the number of available fixed points 
with SU{2) X U{1) isometry depends crucially on the sign of the gravitational Chern- 
Simons coupling co. For a; > 0, there is only one fixed point located at x = 1, whereas 
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for a; < there is an additional fixed point with 



L2 fl2L2 

which is axially symmetric with shape modulus specified in terms of the couplings by 



71 = 72 = ^, 73 = ^, (3.32) 



coL 

(3.33) 



3k2 



Thus, for CO > 0, there is no instanton, whereas for a; < there is a single instanton 
solution with Euclidean action that turns out to be 



_ 47r2 

'instanton — "i T 

\co\ 



a -if (^4fl4 + 8a^ + 12a^ + 2fl + . (3.34) 



Note that the two fixed points coalesce for special values of the couplings associated 
to the choice a = 1, in which case the instanton collapses to a point and its action 
(|3.34[) vanishes, as required on general grounds. For a < 1 (respectively a > 1) the 
anisotropic fixed point is located to the right (respective left) of the isotropic fixed 
point on the x-axis and the corresponding instanton interpolates between a squashed 
(respectively stretched) sphere and the round one. The reduced equation (|3.31|) can be 
easily integrated to express t{x) in closed form, as in [38], but the particular expression 
will not be needed in the present work. 

This construction provides the prime example of an instanton in Hof ava-Lifshitz 
gravity, which has the important property of being chiral for general couplings of the 
theory. The reason is that the Cotton tensor is odd under parity, as it involves the 
totally antisymmetric tensor eijj^ in its definition, and, therefore, it flips sign under 
orientation reversing transformations on S^. Instanton and, likewise, anti-instanton 
configurations only exist for one orientation (sign of a;) but not for the other. This 
occurrence can be intuitively explained by comparing the effect of the driving curva- 
ture terms of the Ricci-Cotton flow. For a; > all curvature terms work in the same 
direction to form the totally isotropic fixed point, but when a; < the Cotton tensor 
competes against the other curvature terms so that more than one fixed points can 
arise. Then, interpolating flow lines between the two become possible. Chiral instan- 
tons are not commonly used in physic^ but we have to work with them now as they 
arise naturally in z = 3 Horava-Lifshitz gravity. 



^We note in passing that analogous configurations are known to exist in some simpler non- 
relativistic integrable field theory models in 1 + 1 dimensions, as in a derivative variant of the non- 
linear Schrodinger equation, which possesses chiral solitons (for a brief account see, for instance, ||58l 
and references therein). The chiral nature of such solitons is also attributed to the presence of Chern- 
Simons terms in a (2 + 1) -dimensional gauge theory that undergoes dimensional reduction and gives 
rise to these integrable systems (see also 159] for more technical details on this subject). 
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It is instructive at this point to provide an effective point particle description of 
the instanton solutions associated to the normalized Ricci-Cotton flow in the mini- 
superspace sector of axially symmetric Bianchi IX metrics. In this case, it can be easily 
verified that the action (|3.24[) of the unimodular Horava-Lifshitz gravity consistently 
reduces to the action 



'eff. - 



dt 



16k4 [ ^ .-4/3 _ ^-3/3 



, H e ^ — e ^ e 

dt ) co^L^ \ 3k^ 3kI^ 



(3.35) 



for a single degree of freedom /3 = logx that ranges in the entire real line. This variable 
is commonly used in the mixmaster model of the universe (see, for instance, ||60ll6ll ) 
and it is required to bring the kinetic energy in canonical form. The effective potential 
follows from a superpotential W which is the action of topologically massive gravity 
evaluated for the class of Bianchi IX metrics. Upon analytic continuation in time, the 
instantons of the point particle model provide our gravitational instanton solutions 
satisfying equation (|3.31|) , which now reads in terms of the new variable /3 

Furthermore, the gravitational instanton action is equal to the Euclidean point particle 
action evaluated for the corresponding instanton solution, thus providing a simple 
derivation of (|3.34|) . 

Next, we plot the effective point particle potential V{fi) appearing in the action 



(|3.35|) for different values of the coupling co (positive or negative) to get a picture of 
the resulting vacuum structure in mini-superspace. The results are shown in Fig.2 
where the origin of the axis j6 = represents the round sphere. When a; — > ±oo, the 
anisotropic scaling exponent of Horava-Lifshitz gravity reduces to z = 2 since the 
Chern-Simons term drops out from the action. In all other cases the theory exhibits 
anisotropic scaling z = 3. The value a; = is also special in that the action is given by 
the gravitational Chern-Simons term alone, whereas the Einstein-Hilbert term drops 
out. This limit is best described by letting Kw — ^ oo so that the effective potential 
remains finite. 

Fig.2a depicts the potential V{fi) for infinitely large |a;|, in which case the flow 
equation we are considering becomes the normalized Ricci flow that was studied long 
time ago in all generality (we also refer the reader to the work [621 for its reduction to 
homogeneous model geometries). In our case, the effective potential has a hill peaked 
at = (log4) /3 where 7i = 72 = 473. Fig.2b depicts the potential for a; = (or better 
to say Kw infinite), in which case the flow equation becomes the Cotton flow |j56l. Here, 
the plot of the effective potential resembles Fig.2a with the difference that the hill is 
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now peaked at j6 = (log2)/3 where 7i = 72 = 273. For all other positive values of co 
the potential looks alike with the hill now being peaked somewhere between the two 
extreme values (log2)/3 and (log4)/3; for this reason V{^) is not drawn separately. 

Next, we consider the case of negative co (with respect to a given orientation of 
S^), which is different because the potential exhibits two vacua located at j6 = 1 and 
j6 = — logfl, using the parameter a given by equation (|3.33|) . Fig.2c represents the 
typical form of the potential for a < 1, whereas Fig.2d depicts the potential for a > 1. 
As the parameter a approaches 1, the two valleys (and, likewise, the two hills) are 
beginning to merge before they cross each other. For a = 1 the potential has a single 
hill and its plot resembles Fig.2a. Finally, when a becomes larger and larger the left hill 
is pushed further and further away to the left until it is completely taken over by the 
infinitely steep wall of the potential arising as j6 — ^ — 00. The fact that the anisotropic 
ground state is not accessible anymore as a 00 is also apparent from the instanton 
action (|3.34|) which becomes infinite; in that case one returns back to the Ricci flow 
potential represented by Fig.2a, as noted before. 
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In Fig.2, the red bullets stand for the fixed points of the flow equation (|3.31|) , 
whereas the green star represents the additional fixed point that has been added at in- 
finity by formulating the same equation in terms of the variable j6 given by (|3.36|) . We 
note in this respect that when x — > oo the speed dx/dt approaches a limiting constant 
value ^4k^/3k^L^ as seen by equation (|3.31|) , whereas in terms of the variable /3 we 
have d^ / dt = at infinity In either case, the potential looks flat at infinity Its value is 
simply lowered to in the j6-parametrization, thus putting the asymptotic configura- 
tion on par with the other degenerate vacua. Only when becomes infinitely large 
there is no height difference between the two parametrizations as illustrated by the 
red bullet placed at infinity in Fig.2b for the effective potential of pure Cotton flow. 
The configuration associated to the point jS = +00 is a completely squashed sphere, 
which is nevertheless regular as it does not exhibit curvature singularities. It arises 
in the correlated limit keeping the volume of S'^ fixed. The components of its Ricci 
curvature tensor are Rn = 1 = 1^22/ = = R^^ and R33 = = R^^, whereas 
all components of the Cotton tensor vanish. This degenerate configuration also plays 
important role in general relativity as it describes a bolt in the zero volume limit. 

It is customary to think of instantons as describing the motion of a point particle in 
the inverted potential — V(fi) that follows by analytic continuation of time. Att = —00 
the particle is located at the top of a hill (location of red bullet in the inverted potential) 
and then starts rolling down until it comes to a complete stop at the top of a nearby 
hill at t = +00. If the potential develops asymptotically a flat direction, the point 
particle will also roll down the hill and reach the plateau after infinitely long time. 
The limiting speed of the particle depends of the height difference between the hill 
and the plateau, which it is zero for V{fi). The height difference is not zero when the 
potential is written in terms of the variable x, in which case the particle will reach 
infinity with finite speed. This solution can also be regarded as a bona fide instanton 
solution with finite action. In all cases, an anti-instanton corresponds to a particle that 
evolves reversely in time and has the right initial velocity to reach the top of a nearby 
hill without overshooting or falling short of speed after sufficiently long time. Based 
on this picture, it is then natural to count as instantons not only the flow lines of the 
Ricci-Cotton flow that interpolate between vacua marked with red bullets but also 
between configurations marked with a red bullet and a green star, as shown in Fig.2. 

In our previous work [38] we have not included the latter solutions in the list of 
instantons with SU{2) isometries, because we were only considering flow lines that 
interpolate between two different fixed points. But we may also add to list those 
configurations that extend from the fully squashed sphere (green star) to a nearby 
fixed point (first red bullet appearing to the left), since they also describe complete 
and regular metrics on R x with — cxj < ^ < +00 (although the extrinsic curvature 
of their slices does not vanish at both ends). The construction of instanton solutions 
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in the unimodular phase of Hof ava-Lif shitz gravity is now complete and can be used 
further for general purposes. Note, however, that these additional instantons will 
not play substantial role in the present work, because the Ricci scalar curvature of 
their spherical sections is everywhere non-negative. But they have some common 
elements to the instantons of ordinary Einstein gravity, which also admit an effective 
point particle description and reach asymptotically a flat direction. These remarks will 
become more transparent in section 5. 



3.4 Sectional Ricci curvature of instantons 



For later use, let us consider in detail the variation of the Ricci scalar curvature R of 
the three-dimensional spatial slices Z3 of the gravitational instanton solutions. First, 
restricting attention to axially symmetric spheres with coefficients (|3.30|) , we find that 



R 



(4x' 



(3.37) 



Fixing the volume of S^, which is determined by L, we observe that R as function 
of the shape modulus x attains its maximum value in the fully isotropic case x = 1. 
Squashing of the sphere corresponds to x > 1, in which case R{x) decreases mono- 
tonically and it becomes asymptotically zero in the degenerate limit of a completely 
flat sphere. Stretching of the sphere corresponds to x < 1, in which case R{x) also de- 
creases monotonically but without a lower bound. There is a critical value x^ = 1/4 
where the curvature vanishes and then turns negative by stretching the sphere more 
and more. The infinitely stretched sphere that arises as x ^ has infinitely nega- 
tive curvature and in this sense it is not a regular configuration, unlike the infinitely 
squashed sphere that is regular. The range of allowed values of R(x) is depicted in 
Fig.3. 



R(x) 




Figure 3: Dependence of curvature of Berger spheres upon squashing and stretching 
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Clearly, R changes monotonically along the flow lines (|3.31|) , since the instantons 
have either x > 1 or x < 1. For this it suffices to recall that x = 1 is a flxed point and, 
as such, it can never arise as intermediate configuration. Also, x{t) is a monotonic 
function of time (the same also applies to /3(f)) as can be easily seen by following the 
point particle motion in the effective potential. The same conclusion is also reached 
by computing dR/ dt for the special trajectories (|3.31|) . 

It is important to realize at this point that the chiral instantons of z = 3 Horava- 
Lifshitz gravity, which arise for a; < 0, will develop spherical sections with negative 
curvature after some time, if the bound > 4 is satisfied. On the other hand, chiral 
instantons with relatively small coupling co satisfying < 4 have positive sectional 
Ricci scalar curvature for all time t. Then, by Lichnerowicz's theorem fTTj] . there can 
be no harmonic spinors on the spherical slices unless negative curvature is allowed 
to occur by tuning the parameters of the theory. As will be seen later, this suffices to 
show that the index of the four-dimensional fermion operator vanishes on all back- 
grounds with positive sectional curvature. Horava-Lifshitz gravity also admits in- 
stantons with non-vanishing fermion index because there can be solutions which are 
suspended from very stretched spheres with sufficiently negative Ricci scalar curva- 
ture; in fact, it will turn out, that this can only happen provided that > 16. Such 
extreme geometric deformations of space are not allowed to occur in the instanton 
slices of Einstein gravity, which always have positive sectional curvature. 

Horava-Lifshitz gravity with more general values of the superspace parameter A 
also admit instantons. These solutions, however, do not permit extreme geometric 
deformations of their spatial slices. In that case, the fixed points of the unnormalized 
Ricci-Cotton flow have Ricci scalar curvature R = 6Aw which is determined by the 
size of the three-dimensional cosmological constant Aw Then, all SU{2) x U{1) in- 
stantons that arise for A^ > (assuming that A < 1 /3 so that the superspace metric is 
positive definite) have everywhere non-negative sectional Ricci scalar curvature. This 
can also be inferred from Fig.3 (but see also [38J which explains why we are not con- 
sidering such instanton solutions when A^ < 0). The result should be contrasted to 
the normalized Ricci-Cotton flow whose fixed points are three-geometries with dif- 
ferent Ricci scalar curvature, in general, and, as a result, the interpolating instantons 



^When < 4, all SU{2) instantons of the theory have enhanced SU{2) x !J(1) isometry. When 

> 4, there are more general SU{2) instantons that interpolate between the two fixed points without 
necessarily having an extra 11(1) symmetry. These instantons have a modulus which can be adjusted 
to yield the axially symmetric instanton solution we are considering here. We also note that when 

> i there are two additional totally anisotropic fixed points with R = which coexist with the 
axially symmetric point l|3.32t : then, there can also be SU{2) instantons that interpolate between any of 
these totally anisotropic fixed points and the axially symmetric or the totally isotropic ones, but none 
of these has enhanced symmetry. More details on this and related matters can be found in 1 38 1 . None 
of these more exotic possibilities will be further discussed here as we are only considering instantons 
with SU{2) X U{1) isometrics. 
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have slices that may extend from positive to negative values of R. This also explains 
why we are only considering the case of normalized flow while searching for models 
of chiral symmetry breaking induced by gravity. 

The above considerations were based on the large isometry group of the specific 
instanton solutions. In order to make more general statements for the variation of the 
Ricci scalar curvature on L3, it is necessary to compute the time derivative of R under 
the flow equation 

dtgrj = X,j, (3.38) 

where X/y is the driving term of the deformations. For example, for = —2Rij we 
have the Ricci flow, for X/y = Qj we have the pure Cotton flow, and for an appropri- 
ate combination of Ricci and Cotton curvature terms we have the more general flow 
equation (|3.18)) . Since dtg^^ = —X^i, we find that the Christoffel symbols of the metric 



evolve as 

^t^tj = \f (V,X^7 + V^X,/ - V/X,y) , (3.39) 
which in turn yield the time evolution of the Ricci scalar curvature, 

■d^R = VrVjX'j - AX - RijX'j , (3.40) 

where X = g'^X/y is the trace of the driving curvature terms and A = VfV is the 
Laplace operator with respect to g. 

Ideally, we would like to use (|3.40|) to put bounds on R and examine its behavior 
against the possible formation of singularities. For the Ricci flow equation we have 
Xjj = —2Rij, and, therefore, {dt — A)R = IRjjR^^ which is manifestly positive definite. 
Then, in this case, if R is positive at some initial time, it will stay positive throughout 
the evolution. Also, by appealing to the maximum principle, we have more generally 
that ii R > c at t = (for any c G R), then the curvature will stay bounded for all 
subsequent times, as R > c/[l — {ct/3)]. Also, for c > 0, it turns out that the Ricci 
flow becomes extinct since the three-geometry develops a curvature singularity at fi- 
nite time T < 3/c (for this and related issues see, for instance, the textbooks Il54| and 



[|55| and references therein). The normalized Ricci flow is straightforward to study 
too. Similar estimates can not be presently made for higher order flow equations (Cot- 
ton or Ricci-Cotton flows) because the maximum principle does not apply to them 
directly. Thus, without imposing any isometries, it is not at all clear what are the gen- 
eral criteria for having eternal solutions of the Ricci-Cotton flow and how R varies 
along the flow lines without hitting a singularity. 

These last remarks show our current mathematical limitations to study the Ricci- 
Cotton flow equation on general grounds. They also make us appreciate the existence 
of simple instanton backgrounds for the propagation of fermions coupled to gravity. 
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4 The index of Lif shitz f ermion operator 



We are now in position to combine the results of the previous sections and consider the 
index of the Dirac-Lifshitz operator associated to the integrated form on the anomaly 
on gravitational instanton backgrounds of z = 3 Horava-Lifshitz theory. The compu- 
tation can be easily done based on the spectral flow in Euclidean time t and the result 
turns out to be identical to the index of the Dirac operator on these backgrounds, as 
suggested by the local form of the anomaly. We make some general remarks first and 
then specialize the discussion to fermions propagating on backgrounds R x with 
axially symmetric homogeneous geometries on S"^. This restriction allows us to il- 
lustrate some special features of the non-relativistic fermion-gravity models, based 
on explicit mathematical results for harmonic spinors fVT] on Berger spheres ||39l |40l 
that are neatly summarized in Appendix D. It also allows to compare them directly 
to known results for relativistic field theories. The implications will be discussed to- 
wards the end of this section together with some more general remarks about non- 
relativistic theories of gravity in the deep ultra-violet regime. 

We consider the coupling of massless Lifshitz fermions (|2.1|) to gravity given by 



the following non-relativistic action, which is written in ADM form. 



k2 SW ^ SW 



, (4.1) 



using the action of three-dimensional topologically massive gravity on E3 with metric 
gij as superpotential functional W and the extrinsic curvature K/y of the spatial slices 
E3 in space-time. As such, it exhibits anisotropic scaling z = 3 in both fermion and 
boson sectors. For the extended theory (|4.1|) , a systematic coupling constant expansion 
may be given in which the Fermi fields enter only at the first correction to the lowest 
approximation. Thus, according to the standard lore (see, for instance, ||63B , as well as 
|[64| section ILG), the backreaction of the Fermi fields can be consistently ignored in 
the Bose field equations and solve the Fermi equation in the background of an external 
Bose configuration. This is the approximation that will be adopted through out this 
section while studying fermions in gravitational instanton backgrounds. 

The chiral nature of instantons in Hof ava-Lifshitz gravity should be contrasted to 
the behavior of gravitational instanton solutions in Euclidean Einstein theory under 
parity. In ordinary gravity, which preserves parity, a change of orientation in space 
simply exchanges instantons with anti-instantons, and, as a result, the background 
geometry shows no chirality. Based on this difference, it is then natural to inquire 
whether there is violation of chiral symmetry for Lifshitz fermions in the background 
of such chiral gravitational instanton solutions in sharp contrast with the negative re- 
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suit obtained for Dirac fermions in the background of Einstein gravitational instantons 
(with boundaries) [|18li24ll25ll29| . This question is addressed in detail in this section 
and the answer depends crucially on the geometry of the three-dimensional spherical 
slices that the instanton transverses as t varies from — oo to +00. For Berger spheres, 
in particular, for which explicit results can be obtained relatively easy, the index of the 
four-dimensional fermion operator is non-zero provided that the instanton slices can 
become sufficiently stretched. In fact, as it turns out, there is a critical value in the 
couplings of the theory beyond which the background geometry becomes sufficiently 
chiral to allow level crossing in the spectrum and lead to chiral symmetry breaking by 
the axial anomaly 

4.1 Atiyah-Patodi-Singer index theorem 

Let us consider the fermion operator D = i^^T>^ for either Dirac fields = D^) 
or Lifshitz fields {V^ = E)^), which are treated here together in a unified way. They 
act on four-component sptnors Y(f, x) which we split into two subspaces of chirality 
±1 with respect to the operator 75, i.e., 75 Y± = ±Y±. Introducing the pair of cor- 
responding spin bundles A± over the space-time M4 and assuming that the first and 
second Stiefel-Whitney classes vanish so that the manifold is orientable and it admits 
a spin structure, we define the spin complex by 

V: C~(A+) ^C~(A_), 

C~(A_) ^ C~(A+) . (4.2) 

Then, the index of the spin complex (also called index of the operator V) is 

Ind(P) = dim Ker V - dim Ker P+ = n+ - W- , (4.3) 

where n± is the number of normalizable solutions of the equation i^^T>^^ = with 
chirality ±1, respectively. These definitions are analogous to the familiar spin com- 
plex associated to the ordinary Dirac operator. In either case, the operator V^I) is 
elliptic on spaces with Euclidean signature which will be considered throughout this 
section (though the order is different for the Dirac and Lifshitz operators), and, thus, 
the general mathematical framework for computing the index by analytic means is 
applicable. 

The Atiyah-Singer index theorem on compact manifolds without boundary | [T2| 
states that n+ — n_ equals the integral of the A-roof genus on M4. In physical terms, 
this is provided by the local form of the axial anomaly, which is the same for both 
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Dirac and Lifshitz fermions, and the index theorem states 



lnd{V) = -^l Tr(f A f ) - j Tr(K A R) (4.4) 

when both the gauge and metric field contributions are combined together (it is called 
the twisted spin complex); here we are only considering the simplest case of a Lf(l) 
gauge field, but, apparently, the result easily generalizes to non-Abelian gauge fields. 
Thus, when 9M4 is empty, the integrated form of the anomaly yields the relation 
Ind(D) = Ind(D). This is an interesting result in its own right, but it is not directly 
applicable to the cases we have in mind. First of all, the Lifshitz models of gauge theo- 
ries that are available in the literatur exhibit anisotropic scaling z = 2 and not z = 3. 
But even if this mismatch is not a major concern in the coupled spinor-vector Lifshitz 
theory, the problem is that the space-time has boundaries as the theory is naturally de- 
fined on M4 ~ R X L3. Furthermore, we have no instanton solutions in our disposal 
for Lifshitz gauge theorie^ in order to investigate the relation between the index of 
the fermion operator and the instanton number of the background configurations, as 
in relativistic Yang-Mills theories. For all these reasons we will not consider the gauge 
field sector in the remaining part of this paper. As for the gravitational sector, we need 
to work with the generalization of the Atiyah-Stnger index theorem to manifolds with 
boundaries in order to account correctly for the end-points of our instanton solutions. 

The index of the fermion operator T> in the background of a four-dimensional 
metric on a manifold M4 with boundary 9M4 takes the following form, according to 



^*^The Lifshitz gauge field theory (Abelian or non-Abelian) is naturally defined by the following 
action, which is written here in flat space-time in 3 + 1 dimensions, |37|, 



dt d^x [Tr(E,E') - ^Tr (^{DiF'^){DjF!i,] 



where E, = 3f A, and F/y = 3jAy — djAi — i[Ai, Aj] provide the electric and the magnetic fields, respec- 
tively, in the axial gauge Ao(f, x) = 0. S satisfies the detailed balance condition using as superpotential 
functional W the action of three-dimensional Yang-Mills theory with coupling constant g. Augment- 
ing W by the Chern-Simons term for the gauge fields A,, so that the three-dimensional action is that 
of topologically massive gauge theory |50l , simply shifts D/f ''^ by a term proportional to e'^^Fjj in the 
potential term of S. In either case, the anisotropic scaling parameter can not be higher than 2 contrary 
to the gravitational case for which the addition of a Chern-Simons term increases z from 2 to 3. 

^^The instantons equations for Lifshitz gauge theories can be derived by completing the square of 
the Euclidean action S, as we did for the gravitational case. They correspond to eternal solutions of the 
first order equations in time 

dtAi{t,x)=±(DjFii + pieijkFi'^] 



taking also into account the presence of a Chern-Simons term in W with coupling constant pi. This 
equation is a variant of the gauge field flow in three dimensions, which, however, has not been studied 
in detail so far. We hope to return to it elsewhere with explicit results. 
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the Atiyah-Patodi-Singer index theorem ||T9ll (but see also lfT3ll20| ), 



Ind(D) = --^ [ Tr(R A R) + [ Tr(0 A R) - ^Vvm,) . (4.5) 

19271^ JMi 19271^ JdMi 2 

In writing this formula it is implicitly assume that we are computing the L^-index 
of the fermion operator acting on square integrable spinors on the manifold M4 with 
respect to the inner product 

(Y, Y') = / rf^xVdetG Y+(x)Y'(x) . (4.6) 

For manifolds of the form R x E3, in particular, with compact slices without bound- 
aries, this means that we are considering spinors with asymptotic behavior 

Y - e~^+^ as ^ ^ +00 with fc+ > , (4.7) 

Y - e~''-^ as ^ ^ -00 with < , (4.8) 

which is characteristic of bound states. These are often called Atiyah-Patodi-Singer 
boundary conditions, which are assumed from now on. When the slices E3 are com- 
pactified due to the physical boundary conditions imposed on the fields at infinity, 
the index theorem remains valid as it stands. However, for non-compact spaces that 
account for permanent boundary effects, one has to device a variant of the index the- 
orem that is applicable to open spaces |65l , as in the case of monopole configurations. 
None of these generalizations will be further pursued in this paper as we restrict at- 
tention to Z3 ~ S"^ exclusively. 

The bulk contribution to the index (|4.5|) is given by the integral over M4 of the 
characteristic class Tr(R A R), which is appropriately normalized, as before. In the 
presence of boundaries one adds two surface terms. The first surface term is the inte- 
gral over 9M4 of the secondary characteristic class Tr{9 A R) that is computable from 
the second fundamental form 9 determined by the choice of the normal to the bound- 
ary, and, as such, it accounts for the difference of the metric from a cross-product form 
at the boundary. Introducing a cross-product metric Gq on M4, which agrees with the 
original metric G on dM^ and serves as reference frame in the calculations, we then 
define the second fundamental form as the difference of the corresponding connec- 
tion one-forms 6 = co — co^. By definition o;*^ has only tangential components to dM^. 
Furthermore, we have Tr(R A R) = dQ{iv,iv^), where 

Q(a;,a;°) = Tr (lO AR+'^9 AO A9 - 29 Acv A9 - 9 Ad9^ (4.9) 



is the Chern-Simons form [|5T| and Tr(Ro A Rq) = for the cross-product reference 
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metric. A simple calculation yields 



Q(a;, w") [boundary = Tr(e A R). 



(4.10) 



Then, the surface correction term has to be subtracted from the bulk term with the 
same normalization in order to have a well-defined quantity when the metric is not 
of a cross-product form. Finally, the last surface term is the so called //-invariant 
determined by the eigen-values of the tangential part of the operator restricted to the 
boundary. It is a non-local term that accounts for the asymmetry between the positive 
and negative modes on 9M4 and it is commonly defined using the spectral Riemann 
zeta function. It is also corrected by the dimension of the zero modes if they are present 
at the boundary. In all examples that will be considered later the //-invariant will be 
computed explicitly by algebraic methods. 

When the index is computed on a given metric background, the bulk term and 
the first boundary correction are identical for the Dirac and Lifshitz operators on M4 
because the normalization constants are the same. We claim that a much stronger 
statement is actually true, namely that the indices of the two operators are equal 



on spaces with boundaries. It means that the //-invariants of the tangential part of the 
Dirac and Lifshitz operators are the same even though their eigen-values differ. We 
will prove this relation by direct computation only for space-times whose slices have 
homogeneous and axially symmetric metrics. This is sufficient to demonstrate that 
the index of the Dirac operator can still be used to determine whether there is chiral 
symmetry breaking on a given background with SU{2) x U{1) isometry in the non- 
relativistic theory of Horava-Lif shitz gravity. Note, however, that this is also sufficient 
to establish the general validity of the relation //d = //jd, beyond mini-superspace 
models, without putting extra effort. According to the general theory (see, in partic- 
ular, fl9]), the //-invariant should differ by a constant from the Chern-Simons invari- 
ant (appropriately normalized). Since the local form of the anomaly yields the same 
Chern-Simons invariant in both cases, it suffices to show that the constants are also the 
same. This can be easily established using the standard metric on S'^, as will be seen 
later, and, therefore, the //-invariants of the Dirac and Lifshitz operators ought to be 
equal, as they are, for all metric^ll. Verifying this general result for all Berger spheres 
seems superfluous, but we will go through it anyway. It serves as an independent 
consistency check based on spectral methods that the local form of the gravitational 

^^Care should be taken when the operators at the boundary admit zero modes, in which case one has 
to add the necessary harmonic corrections to rj. 




(4.11) 
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anomaly is the same for Dirac and Lifshitz fermionj^ 

The computation of the index can be alternatively described as spectral flow of the 
three-dimensional Dirac operator on the slices E3. This approach is very important 
for understanding intuitively the meaning of the index. It also puts in better context 
the results that will be given later. Let us assume first, for simplicity, that the metric 
on M4 has the cross-product form 

ds^ = N{tQ)dt^ + gij{to, x)dx'dxj (4.12) 

around t = tQ (we may also set N = 1 without loss of generality). Then, the fermion 
operator takes the block-diagonal form (see Appendix A for the fermion operator in 
Euclidean space, which can be either Dirac or Lifshitz) 



j d/dt - iaiEj% \ 

(4.13) 

\ d/dt + iaiEfVi / 



where cj are the Pauli matrices and Ef are the inverse dreibeins associated to the 
metric gij (summation over the space indices i and the tangent space indices I are 
implicitly assumed). The 2x2 blocks d/dt ± iaiEiVi are operators mapping the two- 
component Weyl spinors Y± to Yrp and they are mutually related by conjugation as 
{d/dt + iaiEi'ViY = -(d/dt - iajEj'V,). Equation ij^V^Y{t,x) = acting on four- 
component spinors reduces to the following system of Weyl equations on M4, 

^ ± mEi'V,^ Y±(f,x) = . (4.14) 

The cross-product form of the metric on M4 allows to introduce the separation of 
variables Y±(f,x) = exp{—Et)W±{x) and obtain 

± iaiEi'ViW± (x) = EY± (x) (4.15) 

in the vicinity of ^ = tQ. Thus, the solutions of the four-dimensional fermion equa- 
tion i^^T>yW{t,x) = reduce to the eigen-value problem of the three-dimensional 
fermion operators ±icrjEj'T?i, which is analogous to the familiar reduction of the time- 
dependent Schrodinger equation to the time independent one (though, here, the time 



^^We intend to revisit this problem elsewhere by different methods based on the integral representa- 
tion of the //-invariant of an operator D, 



\s=0 
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is Euclidean). 

In reality, however, M4 may not have a cross-product metric for all t and the three- 
dimensional operators Vj may also depend upon t. To remedy the situation we make 
the simplifying assumption that the time evolution is very slow so that in the adiabatic 
approximation we may write 

Y±(i,x) = exp dt'E{t')^ Y±(x) (4.16) 

with 

± iaiEi%{t)'f±{x) = E{t)Y±{x) (4.17) 

replacing the naive relations given before. In this more realistic (but not yet ultimate) 
picture, the solutions of the four-dimensional fermion equation still reduce to the 
eigenvalue problem of the corresponding three-dimensional fermion operators but 
with time dependent eigenvalues E(t) that are slowly varying with time. Normaliz- 
able solutions x) on R x S'^ correspond to those cases that E{t) is positive for 

i — ^ +00 and negative for i — ^ — 00 (these are the Atiyah-Patodi-Singer boundary 
conditions (|4.7|) and (|4.8|) with k± = E(t = ±00)). This, in turn, implies that the three- 



dimensional fermion operators ±icrjEj^T>i should admit zero modes for some interme- 
diate values of i in order for the spectral flow E{t) to be able to produce the necessary 
level crossing from negative to positive eigen-values. Otherwise, there will be no 
normalizable solutions of the four-dimensional fermion operator ij^^V^ and its index 
will be zero. In either case, the eigen-value problem zcr/Ej'DjYi (x) = ±E{t)W± (x) is 
central to the subject and level crossing from negative (resp. positive) eigen-values to 
positive (resp. negative) eigen-values is the dominant effect as t varies continuously 
from —00 to +00. Every normalizable solution contributes one unit to the index (ei- 
ther to n+ or n_ depending on the chirality), as required on general grounds. Letting 
t — ^ —t in the applications, so that instantons and anti-instantons are interchanged, 
simply amounts to replacing n± by tizp. 

Finally, it can be shown that the index does not depend on the adiabatic approxi- 
mation used to establish its relation to the difference of the normalizable zero modes 
of the the operators A = d/dt + iaiEiV^ and = —d/dt-\- iaiEj^Vj. It is obvious 
that the zero modes of A are also zero modes of A^A, and, likewise, the zero modes 
of A^ are also zero modes of AA^ . These two sets are in general different, whereas 
the non-zero modes of A^A and AA^ coincide (for this note that if Y is an eigen-state 
of AA^ with eigen-value E 7^ 0, the state E~^A^Y will be eigen-state of A^A with 
the same eigen-value E). Then, the difference of the number of zero modes can be 
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equivalently expressed by 

n+-n- = dim Ker A - dim Ker = dim Ker (A+A) - dim Ker (AA^) . (4.18) 

The operators A^A and AA^ are better behaved than A and A^, since they elliptic for 
either Dirac or Lif shitz cases, and they are more suited for the argument we are about 
to make. Smooth variations of the background metric that keep the boundaries invari- 
ant deform the entire spectrum of the operators used in the adiabatic approximation, 
but they can not affect the index. Indeed, if a non-zero eigen-value of A^A deforms 
to zero, the same eigen-value of AA^ will also be forced to become zero. Then, both 
n± will increase by 1 and their difference with stay the same. By the same token, any 
non-zero mode of A^A that crosses from negative (resp. positive) to positive (resp. 
negative) values will be accompanied by the same level crossing of the operator AA^. 
Hence, the result of the adiabatic approximation extends without modification to the 
index of the fermion operator D in the background of instanton metrics, which inter- 
polate continuously between the two end-points of time, irrespective of the interme- 
diate details. 

For a more detailed account of the index theorem and related subjects we refer the 
reader to the original papers | IT2l[T9| and also to the mathematics textbooks | [T3l[T4ll20| 
and the physics report ||29|. An excellent account of the physical applications of the 
index theorem to fermions in topologically non-trivial background fields is provided 
in the modem textbook ||66| (see, in particular, part III) with many references to the 
original works. It is beyond our purpose to list separately all important contributions 
made in the literature in this area of research over the years. 



4.2 Lifshitz operator on Berger spheres 

We will determine the spectrum of the three-dimensional Lifshitz operator and com- 
pute its ^/-invariant for the class of axially symmetric Bianchi IX metrics. Our pre- 
sentation is parallel to that for the Dirac operator. Thus, we consider the third order 
operator 

iyiD, = 1 [(zVD0(-D2) + (-D2)(zVD,)] (4.19) 

with = DiD\ which is the tangential part of i'y^'D^ restricted to any given slice of 
space-time (eventually it will be taken at the boundary), and examine the eigen-value 
problem 

zYD,Y(x) = z(rjEj'DjY(x) = Z Y(x) (4.20) 

for two-component spinors Y(x). This equation admits zero mode solutions (the ana- 
logue of harmonic spinors for the three-dimensional Dirac operator) when the geom- 
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etry has sufficient negative curvature. These modes are very important in our dis- 
cussion as they provide the critical points for level crossing under geometric flows, 
which, here, will have the interpretation of spectral flow. Note in this respect that the 
number of zero modes is not a topological invariant as their existence solely depends 
on the geometry of space. 



Using the three-dimensional analogue of Lichnerowicz's formula (|B.9|) for the square 
of the Dirac operator, we may substitute 

- d2 = {iYD.f - 1r (4.21) 
into the Lifshitz operator (|4.19)) , which now takes the following form 

zVD, = {ij'D.f - iHiYD,) - hyiV.R) . (4.22) 

This equation is valid in general for all metrics on a given manifold E3 and forms 
the basis of our construction. Note that the Lifshitz operator does not necessarily 
commute with the Dirac operator since there is a geometric obstruction 

[iYD,, ij^lDj] = ^(V'R)D, + ^V^R (4.23) 

that involves the derivatives of the Ricci scalar curvature R on E3. Here, is the 
three-dimensional Laplacian acting on scalars as opposed to that acts on two- 
component sptnors. 

When the obstruction term (|4.23|) vanishes the eigen-value problem of the Lifshitz 



operator reduces to that of the Dirac operator. A class of metrics that realize this 
possibility is provided by homogeneous model geometries because R becomes purely 
algebraic: in these cases the curvature depends on the metric coefficients, since the 
geometry is not isotropic in general, but it is independent of the coordinates on E3 
(it is the same at all points) as simple consequence of homogeneity. Then, the three- 
dimensional Lifshitz operator simplifies to 

iyiD, = in%f - ^Rih'D,) (4.24) 

and the problem of its diagonalization reduces to that of the three-dimensional Dirac 
operator. In fact, any class of metrics arising in Bianchi's classification of homoge- 
neous three-geometries will serve this purpose, but we will not delve into the general 
discussion. It is also appropriate to note at this point that the use of homogeneous 
geometries also removes the factor ambiguity in the definition of the Lifshitz fermion 
operator, which is now uniquely defined (recall that different factor orderings differ 
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by total derivative terms of the curvature tensor which vanish for all homogeneous 
metrics). 

Here, we restrict attention to Bianchi IX metrics on with isometry group SU{2). 
Still the spectrum of the fermion operator is not easy to determine in closed form 
unless there is an additional symmetry that makes the problem tractable. We consider 
axially symmetric metrics with isometry SU{2) x 11(1), as in section 3, which we write 
in the form 

ds^ = 7 Ua^f + ((r2)2 + ^2(^3)2j (425) 

setting 7i = 72 = 7 and 73 = S^y. The Dirac and Lifshitz operators scale uniformly 
by y/j and (^7)'^, respectively, under conformal transformations of the metric. Thus, 
we may set without loss of generality 7 = 1. The zero modes, in particular, are inde- 
pendent of 7 and they only depend on S. The entire spectrum can also be determined 
as function of 3. The parameter S varies from to 00, depending on the degree of 
spatial anisotropy, and this induces a spectral flow that will also be studied in the 
following. 

Based on these observations, we immediately see that both the Dirac and Lifshitz 
operators have the same eigen-states Y on Berger spheres, whereas the eigen-values 
Z are simply expressed as 

' " ^ ^ ' (4.26) 



in terms of the eigen-values ^ of the Dirac operator. Recall at this point that the spec- 
trum of the Dirac operator on Berger spheres is given by |39| (but see also Appendix 
D for the notation and essential details) 

1 



^^ = 4^25 V ^^'P^ +iP-'^V' (p. ^) e N . (4.27) 

The multiplicity of these eigen-values isp + q for each pair {p,cj). The positive integers 
p and q are not ordered which means that the conjugate pair {q, p) also yield the same 
eigen-values ^± with the same multiplicity. Of course, p and q can also be equal to 
each other. Also, there are additional eigen-values arising for q = 0, 

with multiplicity 2p. Equally well we could have set p = and let q G N, but there 
is no solution when both integers are equal to zero (^ = ^ /4 is not allowed to oc- 
cur). These formulae provide the complete spectrum of the Lifshitz operator on Berger 
spheres together with the multiplicities. We will denote the eigen-values by Z± and 
Zq depending on the choice ^± and ^0 raade, respectively, in (|4.26|) . 
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The general results simplify in three special cases that are worth commenting as 
for the Dirac case. First, we consider ^ = 1 that corresponds to the homogeneous 
and isotropic constant curvature metric on S^. The spectrum of the Lifshitz operator 
becomes 

^ = 1: Z = ± — (2n + 1) (4n^ + 4n - 5) , n G N with multiplicity n(n + l) (4.29) 

and it is equally distributed to positive and negative values. Next, we consider ^ = 
that corresponds to a fully squashed three-sphere along its third principal axis. In 
this case, the eigen-values Zq become infinite. The eigen-values Z± also tend to ±oo, 
respectively, as long as p ^ q, whereas for p = q = n finite eigen-values arise 

^ = : Z = ±-n{2n^ - 1) , n G N with multiplicity 2n (4.30) 

and they coincide with the spectrum of the Lifshitz operator on the round S^. Finally, 
we consider the limiting case ^ — > oo that corresponds to a fully stretched three-sphere 
along its third principal axis. All eigen-values become infinite and exhibit the follow- 
ing universal behavior, 

^ — ^ 00 : Z ~ with infinite multiplicity . (4.31) 

The spectra in these cases should be compared with the corresponding expressions 
(|D.23|) , (|D.24|I and (|D.25|) found in Appendix D for the Dirac operator. 



Zero modes of the Lifshitz operator arise for those special values of S that Z = 
and, as such, they are invariant under rescaling of the metric. Note at this end that 
the eigen-values Zq and Z+ are positive definite for all 5, whereas Z_ are negative for 
(J < 4. This follows from the fact that + (^2 _ 4^ /g 

in (|4.26|) stays positive definitJ^ 



for all values ^0 arid ^± irrespective of 5. Then, the sign of the eigen-values Z is 
completely determined by the sign of ^. Zero modes arise when ^_ = 0, i.e., when 



^*The validity of this statement, which is very important in the following, can be easily shown case 
by case. First, we consider arid find 



^§ + ^ - 4) = ^ (3^^ + 4(p - 2)3^ + 4p2) 



16^2 

which is manifestly positive definite for all p > 2. For p = 1 the terms in the parenthesis take the form 
33'^ — 4:3-^ + 4 = 23'^ + {3"^ — 2)^ which is also positive definite for all 3. Next, we consider ^± and note 
that ^2 -(- ((^2 — 4)/8 > is equivalent to the inequality 

93^ + 16i2pq - 3)3^ + 8 ((p - + 8{2pq - 1)^) 3^ + 64(p - qfi2pq - 1)3^ + 16(p - q)^ > . 

This is manifestly true for all integer values of p and q apart from p = q = 1 that should be examined 
separately In the latter case, ^± = 3/4:±l and, thus, + {3^ - 4)/8 = (3^^ ±S3 + 8)/16. Since the 
discriminant is negative, the sign is positive for all values of 3, which completes the proof. 
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there are positive integer solutions {p,cj) to the equation 



S^ = 2 y4^^^27^^r^ ^ (4.32) 

as for the Dirac operator (see equation (|D.26|) ). Thus, the necessary condition for the 
existence of Lifshitz zero modes is provided by the bound (5 > 4 on the anisotropy 
parameter of Berger spheres. The first zero modes arise at ^ = 4 with multiplicity 2. 
When the Berger sphere is stretched more and more, zero modes occur at other special 
values of S with the same multiplicities as for the zero modes of the Dirac operator 
described in Appendix D.2 (where we refer the reader for details and examples). 

The presence of zero modes implies level crossing under spectral flow when S 
varies. As soon as an eigen-value Z_ crosses zero, it will stay positive for all higher 
values of S and never cross back to negative values. Fig.4 provides the spectral flow 
of Z_ as functions of 3 for the lowest values of p and q. We plot the results for 
(p^q) = (1,1), (1,2), (1,3), (2,2), (1,4) and (2,3) for which level crossing occurs 
at ^ = 4, 5.67, 6.95, 8, 8.03 and 9.80, respectively. The eigen-values Z- with p = q pass 
from the values —1/2, —7, • • • that correspond to the negative modes of the Lifshitz 
operator on S^, whereas all other ones tend to — oo as ^ 0. Direct comparison can 
be made with the corresponding eigen-values ^_ of the Dirac operator shown in Fig.6 
in Appendix D.2. Of course, the dependence of Z_ upon S is more complicated now. 
All eigen-values behave asymptotically as Z = 33^/64: when S ^ oo (compare to the 
asymptotic linear slope lines ^ = S/Ain Fig.6). 



Z_(6) 




Figure 4: Level crossing of the modes Z- (S) from negative to positive values 
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In anticipation of the general results for the j/-invariant of the Lifshitz operator on 
Berger spheres, we examine the spectral asymmetry between positive and negative 
eigen-values in the special cases ^ = 1 and S = 0. This computation can be done by 
elementary means and can be used for comparison with the more general formulae 
that will be derived next. Since there is no spectral asymmetry for S = 1, //o vanishes 
in this case, as //d also does. For ^ = the finite eigen-values (|4.30|) are equally 
distributed to positive and negative values and the same thing applies to all other 
eigen-values Z± that tend to ±oo. In this case, however, there is an excess of positive 
modes provided by Zq. They have multiplicity 2p and depend upon p as up to an 
overall constant 1/ (8^'^). Although this normalization constant is infinite at ^ = 0, 
the j/-invariant does not depend upon it, since j/o is inert to uniform rescaling of the 
spectrum. We immediately find //d(s) = 2^(3s — 1), and, therefore, //]d = //d(s = 
0) = 2^(— 1) = —1/6. The final result is identical to the //-invariant of the Dirac 
operator for ^ = although in that case the spectrum depends linearly upon p and, 
thus, 77d(s) = 2^(s — 1) (see the relevant footnote in Appendix D.3); setting s = 
yields the same answer 2^(— 1) for both operators. 

These preliminary results make us suspect that //o = J/d holds for all values of S 
and not only for and 1. We will confirm this relation by direct calculation and also 
explain its occurrence in the context of Atiyah-Patodi-Singer theory that was outlined 
before. 

4.3 Computation of the ;y-invariant 

In analogy with the Dirac operator discussed in Appendix D.3, we employ the spectral 
Riemann zeta-function and define the regularized difference between positive and 
negative eigen-values of the Lifshitz operator by the general formula 

^d(s)= E (signZ)|Z|-« (4.33) 

eigenvalues 

excluding zero modes. Then, the ?/ -invariant of the three-dimensional Lifshitz op- 
erator is naturally defined by appropriated analytic continuation at s = 0, as //d = 
//d(s = 0), accounting for its spectral asymmetry. 

First, we consider Berger spheres with ^ < 4 so that Zq and Z+ are all positive 
definite and Z_ are all negative. Then, up to uniform rescaling of the spectrum that 
does not affect the final result for tj^j (recall that when the spectrum scales uniformly 
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the j/-invariant remains the same because it is evaluated at s = 0), we have 



M>0 




+ x 



+ x 



+ 




y + p) -4) 



(4.34) 



setting X = y^AS'^pq + (p — qY for notational convenience. The first line refers to the 
contribution of Z_|-, the second to Z_ and the third to Zq. We will compute the indi- 
vidual sums by expanding all fractions in powers of 5 (up to the appropriate order) 
and then set s = 0. The steps we follow are analogous to those outlined in Appendix 
D.3 for j/d, but, of course, the intermediate expressions are much more involved now. 

The contribution of the modes Zq{5) is the easiest to evaluate. Using the power 
series expansion 




+ p 



(4.35) 



we find that the last term in (|4.34[) assumes the following expansion in terms of Rie- 
mann zeta-functions. 



7o(s) = 2^(3s - 1) + s5^ (4^(3s + 1) - 3^(3s)) + 



s(9s - 1)^^ 



i(3s + 1) + 



(4.36) 



The terms that are omitted vanish at s = since they contain the factor s^(3s + n) with 
integer n >2 and ^(s) is absolutely convergent for Res > 1; they include a term of the 
form S^/p^ in (|4.35|) as well as all terms of order and higher which come multiplied 
with l/p"+i with n > 2. Then, taking into account that ^(-1) = -1/12,^(0) = -1/2 
and that ^(s) has a simple pole at s = 1 with residue 1 (and, thus, 3s^(3s + 1) equals 1 
at s = 0), we obtain 

1 4^2 ^4 



^o(O) 



+ 



12 



(4.37) 



The contribution of the modes Z±{S) to the j/-invariant is more difficult to extract 
because of the double sums that are involved. Again, we expand the fractions in 
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power series of S 



:4- 



~8~ 



^2 

±^ + x 



4) 



9s 



X2 



1 8(3s + 2) 



X3 



X3s 

3s(s-l)(3s + 2)c^^ 
16X3 



+ • • • (4.38) 



omitting all terms of the form l/X""'"^ with n > 2 that do not contribute to the final 
result when s = 0, as will be seen shortly. Then, the first two sums in the general 



expression (|4.34|) for //d(s), which we denote respectively by I±{s), are expanded as 
follows. 



7+(s)-7_(s: 



-3sS^f 



2,/3s + l\ s(3s + 2)^*/ 2A-/3S + 3 



2 J 8 

when they are combined together, setting for convenience 

f(s) = y l±l 



16 + 3(s-l)^^ / 



+ 



(4.39) 



(4.40) 



The function /(s) also arises in the computation of the ^/-invariant of the Dirac oper- 
ator. It is absolutely convergent on the complex s-plane for Res > 3/2, which justifies 
the suppression of the higher order terms in the power series expansion (all such terms 
are multiplied with s and vanish at s = 0). We also recall from Appendix D.3 that /(s) 
has simple poles at s = 1/2 and s = 3/2 with residues given by equation (|D.38|) . These 
are the values of (3s/2)/((3s + l)/2) and (3s/2)/((3s + 3)/2) at s = and they are 
equal to (^^ _ / 5 ^^d 1 / (2^^^ respectively. Taking these into account, we obtain 



I+(0)-I_(0) 



12 



(4.41) 



The final result for the //-invariant of the Lifshitz operator on Berger spheres with 
^ < 4 is identical to the result for the Dirac operator, i.e.. 



is' - If 



(4.42) 



This is not surprising in retrospect because the axial anomaly is identical in both cases 
including the overall normalization constant. Then, by the Atiyah-Patodi-Singer 
index theorem, the gravitational Chem-Simons action should be related to the ?/- 
invariant of the Lifshitz operator exactly as the Chern-Simons action is related to 
//-invariant of the Dirac operator by equation (|D.45|) . From this point of view, the com- 
putation of ?/]D provides an alternative derivation/ verification of the axial anomaly for 
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Lifshitz fermions without path integral methods. Hence our interest in pursuing this 
calculation to the end. It can also be seen here that there are more terms contributing 
to //d(s)/ compared to tj-ois), but they are weighted with different factors. The two ex- 
pressions are clearly different when s ^ (for comparison see the corresponding steps 
taken in Appendix D.3 for the Dirac operator). Nevertheless, the individual terms add 
up to produce the same final result at s = 0, as in the path integral computation of the 
local form of the axial anomaly. 

Next, we consider Berger spheres with S > 4. In analogy with the Dirac operator 
discussed in Appendix D, let us denote by 

C{Sc) = {ip,q) e N2; Si = 2^4pqS^ + {p - q)^} (4.43) 

the set of all positive integers {p,cj) associated to zero modes of the Lifshitz operator 
at each one of the special values < ^- Since the multiplicity of the zero modes at Sq 
isp + q for each {p, q), the quantity 

S(^)=Ef E (P + ^)] (4-44) 

Sc<S {p,q)eC{3c) 

provides the total number of modes that cross from negative to positive values as S 
varies from ^ < 4 to any given value ^ > 4. Taking into account the multiplicity of 
eigen-values, it follows that S{5) is an even integer number. Thus, for S > 4, the rj- 
invariant of the Dirac operator for Berger spheres is shifted by twice the number of 
these modes and equals 

//D = -g(^'-l)' + 2S(cJ). (4.45) 

S(S) is the same number for the Dirac and Lifshitz operators, since level crossing oc- 
curs at the same values of S and the mode multiplicities are also the same. S(S) simply 
vanishes when S <4. Thus, we have 

t]c = t]u (4.46) 

as advertised above. This formula is actually valid for all < J < oo, including the 
special values 3c, where zero modes appear, and it will be used in the following. 



4.4 Chiral symmetry breaking by instantons 

We apply the index theorem to the Dirac-Lif shitz operator defined on Hof ava-Lifshitz 
instanton backgrounds with SU{2) x U{1) isometry Our treatment is quite general, 
with no reference to the details of the particular gravitational solutions, in order to 
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identify the sector of the theory that leads to chiral symmetry breaking. As will be 
seen in the following only the gross qualitative features of the geometry matter for the 
analytic derivation of the index formula. 

Using the Euclidean four-dimensional metric on M4 ~ R x S'^ written in terms of 
the left-invariant one-form cr^ of SU{2) with coefficients 71 = 72 = 7, 

ds^ = dt^ + j{t){(r'f + ^{t){a^f + ^3{t){^^f , (4.47) 

we introduce the orthonormal coframe of vierbeins given by the one-forms 

= dt , = (ji , e^ = , e' = ^ cr' (4.48) 

and obtain the corresponding connection one-forms computed by de^ + co'^i, A = 0, 

, ,0 _ 1 ^7 1 ^ n _ 1 ^7 7 ^ n _ 1 ^73 ^3 

1 — ;~ ^ / CO 2 — ;~ ^ / CO 3 — —- ; — 6 , 

^ 2j dt ^ 2j dt ^ 273 dt 

2 _ \/73 1 3 _ x/73 2 ,,1 _ 27 - 73 3 

CO 3 — -— — e , CO I — -— — e , co 2 — --^^ — 1= e • (4.4y) 
27 27 27^ 

Then, the curvature two-forms computed by R^i, = dco"i, + co^c A co'^i, turn out to be 



R\ = A - -r^:j- ^ A e^ (4.50) 

dt^ j ^y/l^dt V7 ' 



,0 _ 1 / d^V^\ .0 . .2 1 d f 73\ 3 1 



^ 2 = ^ e'' A - {^]e^Ae', (4.51) 



,0 _ 1 / ^^^\/73 \ .C3 , 1 ^ /73\ 1 , 7 



^3 = =^V^ e^Ae^ + -^-^ e^AeS (4.52) 

^\ dt^ IJ^dt V7 



and after some calculation the topological density takes the form of a total derivative 



^^The left-invariant one-forms of SJJ(2) should be distinguished from the Pauli matrices ctj used 
earlier to express the three-dimensional fermion operator as i^'T>i = iaiEi'Vj in terms of tangent space 
indices. 
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term, without using any equations of motion. 



(4.56) 



We also introduce a variant of the four-dimensional metric that has the cross- 
product form 



(4.57) 



using a fixed parameter fo instead of t in the metric components. It serves as reference 
frame for the computations in the vicinity of i = Iq. Since the connection one-forms 
of the cross-product metric have components co'^i = and co^j are the same as above 
with t simply replaced by Iq, we infer that the second fundamental form 9 = co — coq 
has the following non-vanishing components at i = ig. 



9^ 



\t=to 



9'. 



\t = to 



dt ' - - 

Then, Tr{9 AR) at any given slice i = io is given by 



9^ 



3 

\t=to cr 



\t = to 



dt 



(4.58) 



(4.59) 



Applying the Atiyah-Patodi-Singer index theorem (|4.5|) to the fermion operator V 
(it can be either D or D) we obtain 



Ind(P) 



1 

12 



i.3(|ios:^l+l--i 



1 

24 



73 

7 

f = + 00 
t= — 00 



t=+co 

+ 

t= — oo 

t=+co 



Vv 



t= — oo 



(4.60) 



after performing the integral in the time variable t and the angular coordinates of S'^, 
using the volume element A A = sm9 d9 A dcp A dtp, which yields a factor 
of 167r^. This result is quite general, as it stands, because it does not depend on the 
specific form of the metric coefficients 7(f) and 73(^) for intermediate times, but it 
relies on their existence as t varies from —00 to +00. 

Furthermore, we observe that the time derivative terms cancel against each other 
without reference to any field equations, and, thus, equation (|4.60|) simplifies to 



Ind(P) = -- 





2 


t=+oo I 




t= + 00 


) ] 


t=-oo 2 


Vv 










t= — oo 



(4.61) 
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The only information needed for computing the index are the boundary values of the 
metric coefficients, whereas the boundary values of their time derivatives are irrele- 
vant. Thus, this formula is equally valid for the instantons of Horava-Lifshitz gravity 
that interpolate between any two fixed points of the Ricci-Cotton flow with ^ = 1 and 
< ^ 7^ 1 as well as for the instantons that extend all the way to the fully squashed 
sphere with ^ = 0. Setting 73/7 = d^, as we did before, and using the general expres- 
sion for the //-invariant on Berger spheres, y]x> = —{p^ — 1)^/6 + 1S{S), we arrive at 
the final expression 



given by the number of level crossings (including their multiplicities) that have oc- 
curred on the spatial slices, via deformation of 5, in the entire history of time (for the 
definition of S(^) see equation (|4.44|) or (|D.47|) ). This is precisely what was expected by 
the spectral flow argument. Instantons and anti-instantons will have opposite signs 
of AS, if it is not zero, since they are mutually related by i — ^ — i. It is also obvious 
now that the index vanishes on all instanton backgrounds whose end-points have 
non-negative sectional Ricci scalar curvature. 

The only instantons that can give rise to level crossing, and, thus, to non-vanishing 
index occur in the unimodular version of Horava-Lifshitz gravity with superspace 
parameter A = 1/3. They are eternal solutions of the normalized Ricci-Cotton flow 
whose end-points are not constrained to have the same Ricci scalar curvature. When, 
in particular, the Chern-Simons parameter 00 is negative (for a given choice of orien- 
tation on S^), there are flow lines that interpolate between the totally isotropic config- 
uration and a Berger metric with anisotropy parameter 



Recall from section 3.3 (to which we refer for the notation) that such chiral instantons 
are associated to the tunneling of an effective point-particle between the red bullet 
vacua of the potential well shown in Fig.2c and Fig.2d; as such, they correspond to 
continuous deformations of the round sphere to configurations with d < \ and ^ > 1, 
respectively. 

The necessary and sufficient condition to have non-vanishing index is provided by 
^ > 4, which singles out Fig.2d with the left red bullet pushed sufficiently far away 
from the origin. This, in turn, amounts to the following inequality among the various 
parameters of the theory. 



Ind(P) 



AS 



(4.62) 




(4.63) 




(4.64) 
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Here, we have expressed the characteristic length scale (size) L of space in terms of the 
volume of S'^ which equals to 2n^L^ and it remains constant in time. Then, the index 
is given by the total number of modes that changed sign as S was varying from 1 to 
fl"^/^ and equals 

Ind(P) = ±S(fl3/2) (4.65) 

with plus or minus signs referring to the way that the cylinder R x is transversed 
(backward or forward in Euclidean time). This is in agreement with the spectral flow 
interpretation of the index of the spin complex. 

In more physical terms, a non-vanishing index implies that the axial charge asso- 
ciated to the (anomalous) conservation law of the axial current , 

Q5= [ dhy^jl{t,x), (4.66) 

is not necessarily conserved but it rather changes as 

AQ5 = Ind(I?) (4.67) 

in all processes mediated by the corresponding gravitational instanton backgrounds 
and it leads to lepton and baryon number non-conservation. The index is an integer 
numbei0 (this is consistent with the fact that M4 admits spin structure), which can 
become arbitrarily large when the parameter a is sufficiently large. In terms of the 
inequality (|4.64|) , it implies that Hof ava-Lifshitz gravity can exhibit chiral symmetry 
breaking when the relative coupling of the Cotton term (^ k^/co) is sufficiently small. 
Said differently, chiral symmetry breaking becomes possible when there is a lower 
bound on the volume of space for fixed couplings and oj, 

-^j . (4.68) 

A Lifshitz universe that exhibits a bounce ||67l |68l |69l may help to realize this novel 
possibility in practice. We intend to study the implications of this scenario elsewhere 
in more detail. 



^^S(i5) is an even integer and thus Ind(X') that counts the difference of positive and negative chirality 
zero modes of the four-dimensional fermion operator acting on four-component spinors is even. This 
means, in particular, that Weyl f ermions can be consistently coupled to Hof ava-Lifshitz gravity without 
suffering from global gravitational anomalies. Global anomalies can only arise when the number of 



level crossings is odd, in which case the sign of y det{ijf''Dpi) becomes ambiguous fTU\. Apparently, 
there are no large diffeomorphisms of (embedded in the foliation preserving diffeomorphisms of R x 



S^) that can flip the sign of ^ dei{i^T^V^) and trigger an inconsistency. Global gravitational anomalies 

may nevertheless be present in higher dimensional generalizations of the theory, as in ordinary gravity 
in 8fc or 8A: + 1 dimensions, where disconnected diffeomorphisms can arise and play important role 1391 . 
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Another interesting question is the behavior of inequality (|4.68|) for chiral sym- 
metry breaking under the renormalization group equations of topologically massive 
gravity. The one-loop beta functions for the three dimensional Newton coupling 
and the cosmological constant A have been computed explicitly in the recent paper 
llTOl , where it was also found that the Chern-Simons coupling co does not change with 
the energy scale. One may try to recast these results into a form that is directly appli- 
cable to the unimodular variant of topologically massive gravity, in which A assumes 
the role of an integration constant and the volume of three-space appears as a running 
coupling, and inquire whether both sides of the inequality (|4.68|) are renormalization 
group invariant (at least to one-loop level). It seems that the answer is negativjl^ 
A weaker (and perhaps more appropriate) question that can be posed in this context 
is whether the inequality (|4.68|) is not overturned by the renormalization group flow 
once it is satisfied at a given energy scale; and if so, what will be the deeper mean- 
ing of this? We are primarily interested in studying this question in the vicinity of a 
non-trivial fixed point, which is UV attractive in all directions and has negative A, but 
there are still some issues left open in the analysis given in |l70l that prevent us from 
drawing definite conclusions at this moment. Intertwining our results with renor- 
malization group ideas and exploring their meaning in the four-dimensional theory 
of Horava-Lifshitz gravity are some important problems to which we also hope to 
return elsewhere. 



4.5 Topological invariants of the instantons 

The index of the Dirac-Lifshitz operator on our gravitational instanton backgrounds 
is not a topological invariant of space-time. This is mere reflection of the fact that the 
dimension of the space of harmonic spinors on S'^ is not bounded by the topology. It 
only depends on the geometry and it can grow without bound for metrics with very 
negative Ricci scalar curvature. Thus, to complete our discussion, we will compute 
explicitly the topological numbers of the instanton spaces, namely their signature and 
Euler numbers, using the Atiyah-Patodi-Singer theorem for the Hirzebruch and de 
Rham complexes, respectively. They both turn out to be zero, as expected by the 
cylindrical form of M4 (but see also the exact arguments below). 

Hirzebruch signature complex: The signature of a Riemannian four-manifold 
with boundaries is closely related to the formula for the index of the fermion oper- 



ator. It reads |[19| 



t(M4) = -^/ MRAR) + -^[ Tr(0AR)-7s(aM4). (4.69) 



^"^Other important relations in the three-dimensional world, such as the defining relation of the chiral 
point of topologically massive gravity, also fail to be renormalization group invariant, as noted in Ii70i . 
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The first two terms are identical to those appearing in (|4.5|) for Ind(X>), but they are 
now multiplied with 1 / lAn^ that differs from 1 / 192/1^ by a factor of 8 (this is also con- 
sistent with the fact that a compact four-dimensional spin manifold without bound- 
aries has signature that is an integer multiple of 8). Thus, their particular form is 
known from before (see equations (14.561) and (|4.59|) ). The last term is the associated t]- 
invariant for the signature complex^n, which has been computed by Hitchin on Berger 
sphere^ |[7T| and takes the following form, 

7s(^) = -^('52-1)2 (4.70) 

for all values of the anisotropy parameter S. Since it is also 8 times larger than fjjy/l = 
— {S'^ — 1)^/12, without counting the effect of level crossing, as it does not occur here, 
we find immediately that t(M4) = 0. 

Euler-de Rham complex: The Euler number of a Riemannian four-manifold with 
boundaries is given by 

^(^4) = 7^ I eabcd R% A R'd - / eaicd Uh A R'd - \o\ A e\ A d'^ . 

6ln^ JMi 1671^ JdMi \ 3 J 

(4.71) 

The bulk term is the integral of the Gauss-Bonnet topological density and the bound- 
ary term is the integral of the corresponding Chem-Simons secondary class which 
is required for well definiteness. In this case there is no analogue of the ^/-invariant 
boundary term. Explicit computation shows that all our gravitational instanton back- 
grounds obey the following relation, 

R\ A R'^ = rff A ^ e% A R'j - A e\. a e'A = (4.72) 

^ [ 1 f f^I y _ 4, + 3T3 V ^ f ^) 1 rf< A A A 

dtl2j^y\dtj ' 'y\dtj -y^ \dt J \ 

and, therefore, we easily obtain that ;\;(M4) = 0. This relation follows from the particu- 
lar form of the four-dimensional metric ds^ = df- + 7(i)((7^)^ + 7(i)((7^)^ + 73(i)(^'^)^ 
without using any equations of motion. 

^^To define the //s one considers the self-adjoint operator B acting on forms co on the manifold 3M4 
as B{iv) = (— l)P(e*rf — d *) co, choosing e = 1 for 2p-forms and e = — 1 for (2p — l)-forms. B 
preserves the parity of forms on dM^ and commutes with to { — 1)P *lo, so that B = B^^^" ® 6°'^'^ 
and B'^v'^" is isomorphic to B°'^'^. In particular, the (/-invariant of the operator B is twice the //-invariant 
of B*^^*^". Normally, the index theorem l|4.69t would involve one-half the //-invariant of B, but this is 
finally expressed in terms of //g associated to B*^^*^" alone, explaining the factor of 1 /2 that is missing. 

^^The computation is straightforward and it relies (once more) on the relation between the rj- 
invariant and the Chern-Simons action implied by the Atiyah-Patodi-Singer index theorem 1 19|. Note, 
however, that there is a discrepancy by a factor of 2 in the final expression for tjc, compared to the result 
derived in 171,1 (it is apparently a typo). 
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Another way to view these results is to think of M4 ~ R x S'^ as being topologically 
equivalent to a compact manifold Mc — with two points removed, say the "north" 
and the "south" poles. This is best seen by considering the canonical de Sitter metric 
on with radius L and r running from to 00, 



(1 + (r/2L)2)2 



dr2+4r2((£ri)2 + ((ri)2 + ((ri)2 



(4.73) 



Passing to the proper coordinate t given by the change of variables tan(i/2L) = r/2L, 
the metric takes the form 

rfs2 = dt^ + sin2 ( L) ( (a^f + ((7i)2 + {a^)A (4.74) 



with t running from to ttL. It resembles the form of the four-dimensional metrics 
that we have been considering all along having homogeneous and isotropic slices with 
7i = 72 = 73 = sin2(f/L). In this frame, S'* is described as fibration over S^. At 
the two poles which are located at f = and ttL, respectively, the three-dimensional 
slices shrink to a point. Then, removing the two poles yields a compact space with 
boundaries that has the topology of R x S"^. Next, using the following general relations 
t(M4) = t(Mc) and = Xi^^c) - 2, and knowing, in particular, that t(S^) = 

and x{^^) = 2, we obtain immediately that t(M4) = = xi^i)- 

It is now appropriate to compare these results to the topological invariants of in- 
stanton solutions of Einstein gravity, focusing, in particular, to the Taub-NUT and 
Eguchi-Hanson spaces. These metrics take exactly the same form in proper time as 
the instantons of Hof ava-Lifshitz gravity, and, therefore, one may (erroneously) think 
that their signature and Euler numbers also vanish. Here, we focus to the topological 
aspects of these spaces (their analytic properties are discussed in the next section to 
which we refer for details) and explain in brief how the presence of removable singu- 
larities affects the computation of x arid t". Such singularities are not present in the 
instanton metrics of Hofava-Lifshitz gravity, by construction, and, naturally both x 
and T vanish there. 

The removable singularities that arise at one end of the instanton solutions of Ein- 
stein gravity form the fixed point set of a Killing vector field. Then, one has to employ 
the Lef schetz fixed-point theorem, which is a special case of the more general G-index 
theorem, to account for the presence of the fixed points in the computation of the Eu- 
ler number (see, for instance, [29 J for a rather elementary account). According to this 
theorem, an isolated fixed point (nut) adds one unit to the Euler number of the mani- 
fold and a surface of fixed points that form an invariant submanifold (bolt) adds to 
it two units. As a result, the Euler number of Taub-NUT space is 1 and that of Eguchi- 
Hanson space is 2. The computation of the signature requires some additional ingre- 
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dients in the case of Eguchi-Hanson space: the j/-invariant of the signature complex 
(|4.70|) should be made equivariant with respect to a discrete Z2 symmetry that iden- 
tifies antipodal points on (see, for instance, [|7T|) and integration over yields a 
factor of Stt^ instead oilSn^. Taking these into account one finds that the signature of 
Eguchi-Hanson space is ±1 (the sign refers to the choice of instanton or anti-instanton 
metric), whereas the signature of Taub-NUT space vanishes. In either case we have 
xiM^) = |t(M4) I + 1, as required on general grounds for this class of Einstein spaces. 



5 Comparison to relativistic theories 

This section is devoted to instanton solutions of Einstein gravity and the associated 
index formulae for the Dirac operator. It will make more transparent the issues we 
addressed earlier but in a more familiar context. 

The index of the Dirac operator D in the background of a four-dimensional metric 
on a manifold M4 with boundary 9M4 takes the same form (|4.5)) as for the Lifshitz 
operator, according to Atiyah-Patodi-Singer index theorem IIT9B , where t/d is the t]- 
invariant of the three-dimensional Dirac operator at the boundary, which is discussed 
extensively in Appendix D. Ind(D) counts the difference in the number of normaliz- 
able positive and negative chirality spinors obeying the appropriate boundary con- 
ditions, as usual. Our main objective is to compute the index of the Dirac operator 
in the background of Taub-NUT | j2T| and Eguchi-Hanson | j22j l gravitational instan- 
ton metrics so that direct comparison can be made to the results derived in the non- 
relativistic case. Our exposition is somewhat extensive containing all relevant details 
and it includes a relatively new (alternative) derivation of the index in the case of 
Eguchi-Hanson space. 



5.1 Gravitational instantons with 511(2) x 11(1) isometry 

The Taub-NUT and Eguchi-Hanson metrics provide the only complete and regu- 
lar gravitational instanton solutions of Einstein gravity with SU{2) x U{1) group of 
isometries. They also admit an effective point particle description in the context of 
mixmaster dynamics ||60| that parallels the description of S 11(2) x U{1) instantons of 
Hof ava-Lifshitz gravity. 

Recall that the potential of general relativity is minus the Ricci scalar curvature of 
the spatial slices E3 in the ADM formulation of the theory on R x Z3 [|44|. Using the 
Bianchi IX ansatz, we find that the gravitational action with metric coefficients (|3.30|) 
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takes the following formcj in proper time (with N{t) = 1), 



setting for convenience = 32 ttG and 

= logx , n = logL . (5.2) 

Note that the superspace metric is not positive definite for Einstein gravity, since one 
has to take A = 1, f^], and this is reflected in the negative kinetic energy of the vol- 
ume modulus found in the effective action. Because of this difference, the instantons 
of Einstein gravity admit a slightly different description (compared to the instantons 
of Horava-Lifshitz theory), which is not directly based on Euclidean action bounds. 
Now the action can in principle be made arbitrarily negative because of the indefinite- 
ness of the DeWitt metric associated to the conformal factor. Yet, the Euclidean action 
of Einstein instantons is finite and it is given by a boundary term, since the bulk term 
vanishes. Also, the instantons of Einstein gravity have self-dual Riemann curvature 
tensor, whereas in Hof ava-Lifshitz theory self-duality is not present. 

Fig.5 is a plot of the effective potential V{fi) appearing in the action (|5.1|) as function 
of the anisotropy parameter /3 with fixed O. 

V08) 




Figure 5: Effective potential for SU{2) x U{1) mixmaster dynamics in Einstein gravity 

^'^Integration over S'^ yields a volume factor IbTZ'^ that fixes the normalization used in Sgff.. This is 
correct for Taub-Nut space but not for Eguchi-Hanson space whose slices are topologically In 
the latter case the volume factor is Stt^ and the overall normalization of Sgff should be Stt^/k^. We 
overlook this difference here in order to make the presentation uniform, but we will come back to it 
later as it affects the action of the instantons. 



69 



In this case, however, the volume of space changes with time and V{fi) scales appro- 
priately. The plot is closely related to the one shown in Fig.3 for R{x), since V = —R, 
but the potential is now expressed as function of j6 rather than x. It should also be com- 
pared to the form of the potential in Horava-Lif shitz gravity. All these differences will 
be explained next from a particular point of view. 

Following IITSi , we invert the potential V{fi) and represent an SU{2) x U{1) instan- 
ton by a point particle which rolls from the top of the hill located at /3 = towards the 
plateau that is reached asymptotically as /3 — ^ +00. But the potential also varies with 
time through its scaling factor, thus adjusting continuously the height difference in the 
course of the motion. More importantly, the two end-points of the particle trajectory 
have n = ±00. The top of the hill will correspond to a removable nut singularity in 
space-time if O = —00 and to a Euclidean (or conical infinity) if O = +00. Likewise, 
the plateau reached at the other end of the motion will correspond to a removable bolt 
singularity in space-time if Q = —00 and to a Taubian infinity if O = +00. As a matter 
of fact, the self-dual Taub-NUT metric corresponds to a trajectory which starts with 
zero velocity at the top of the hill as O = —00 and reaches the plateau at infinity as 
n = -|-oo. On the other hand, the Eguchi-Hanson instanton corresponds to a trajec- 
tory that starts from the plateau at infinity as O = —00 and reaches the top of the hill 
as n = +00. These are the only two possible instanton solutions of Einstein grav- 
ity with SU{2) X U{1) isometry that yield complete non-singular manifolds. In both 
these cases, the proper time variable t extends from a finite value (say 0) to +00 since 
there is a nut or a bolt forming at i = and space-time ends there. The anti-instantons 
follow by letting i — > — i. 

This qualitative picture can be made precise by writing down the first order equa- 
tions that /3 and Q satisfy in time. They are gradient flow equations derived from a 
superpotential which is the Hamilton-Jacobi functional of general relativity computed 
for Bianchi IX metrics. In our case, the potential takes the form 



J 4 V 9n 



(5.3) 



and the first order equations satisfied by the instantons are 



dl_ sn^J^ ^dn._ ,^dW 

dt~ a/3 ' dt ~ ^ an • ^ ^ ^ 

Note, however, that W is not uniquely determined because of the indefiniteness of the 
DeWitt metric. There are two possible choices of superpotential given by 

Win = ^ (e'^^ - 4.-/^/2 j or Weh = ^ e^"" (e'^^ + 2e^) , (5.5) 
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which correspond to the celebrated Taub-NUT and Eguchi-Hanson instantons, re- 
spectively. In the first case, exp(— 2Q + 2/3) — exp(— 2Q + /3/2) = 1/ (16m^) is a first 
integral of the system (|5.4[) and the integration constant m is called nut parameter. In 
the second case, the integral is provided by exp(4n + 2/3) — exp(4Q — /3) = with a 
being the modulus of Eguchi-Hanson metric. These integrals allow to determine the 
metric of gravitational instantons in closed form. 

It can be easily verified that the first order equations (|5.4|) solve the second order 
equations of motion following from the action (|5.1|) upon analytic continuation in time. 
The relative minus sign between (9>V/9j6)^ and (3W/3n)^ in the potential is also at- 
tributed to the indefiniteness of the DeWitt metric. The plus or minus signs appearing 
in the defining relations (|5.4[) refer to the instanton or anti-instanton configurations. 



which are related to each other by time reversal. These equations also follow from the 
self-duality relations by integrating once in time and making the appropriate choice 
of integration constants (we refer the reader to the original work 11181 for these tech- 
nical details). It so happens that the first order equations for the metric coefficients 
of the self-dual Taub-NUT space are identical to the Ricci flow equations on with 
metric (|3.28|) provided that 7; are traded with 7^. There is no analogous interpretation 
of the equations that determine the Eguchi-Hanson metric. 

The Euclidean point particle action for these trajectories is evaluated by completing 
the squares and using equation (|5.4[) . We obtain 

(5.6) 

which is solely determined by the difference of W at the end points of proper time, 
t = and ±cx). As such, it resembles the form of the action for Horava-Lifshitz in- 
stantons. However, W blows up in the asymptotic region of both instanton solutions 
(Taubian or conical infinity). This can be easily seen by employing the first integrals 
and evaluating W at f = ±00. At the other end of the instanton, which corresponds to 
a removable singularity (nut or bolt) reached at f = 0, W is finite. Thus, Sinstanton/ as 
evaluated above, turns out to be infinite and looks problematic. To remedy the situa- 
tion one has to subtract the infinite contribution of the flat space metric which serves 
as reference frame in the calculations. 

Notice in this respect that the integral of the mean curvature of any given slice E3 
is the derivative of its volume, which, in turn yields the following relation using the 
first order equations (|5.4|) , 

i/^Vd^SK = = |vol(L3) = 2.^^.30 ^ ^3^.^ . 

(5.7) 
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The mean curvature is the trace of the extrinsic curvature, K = g^^Kij, with 2Kij = 
dgij/ dt since we have made the choice of lapse and shift functions N = 1 and N; = 0. 
The volume of space is Vol(E3) = In'^L^ with L = expO, and, finally, we have 
dW/dCi = 2yV due to the particular dependence of W upon Q for both Taub-NUT 
and Eguchi-Hanson spaces. Based on this observation, we recast the Euclidean in- 
stanton action (|5.6|) in the form 

Sinstanton = -4 / d^X^/detgK. (5.8) 

To make the action finite one simply has to subtract the contribution of a reference 
round sphere with volume ISn^t^ at both ends of space-time which we write as 

Sinstanton = "4 / rf^X^di^ {K - Kq) . (5.9) 

The integral oiK — Kq over dM/^ is precisely the so called Gibbons-Hawking bound- 
ary term and it is rewarding to derive it here in the canonical ADM formalism by a 
different line of arguments. Taking this into account, the regulated instanton action 
of the Eguchi-Hanson metric turns out to be zero and the difference in the normal- 
ization of Sgff. (being instead of 16n^) becomes irrelevant. The regulated action of 
Taub-NUT instanton is not zero but finite. 

Recall for completeness that in the conventional approach to the problem one adds 
the Gibbons-Hawking boundary term to the bulk Einstein-Hilbert action \72\ (but see 
also |t73| ). This is provided by the difference between the trace of the second funda- 
mental form of the boundary in a given metric, K, and its value in the flat space metric, 
Kq. More precisely, we have the Euclidean action (setting = 32 ttG, as before) 

Sgravity = -\ I d^xVdetG R[G\ [ d^x^/d^ {K - Kq) . (5.10) 

This does not affect the classical equations of motion but it makes the variational prob- 
lem well-posed in the presence of boundaries. Then, the Euclidean action of the in- 
stanton metrics is provided by the Gibbons-Hawking boundary term alone, since the 
bulk term vanishes, and it is finite, as required on general grounds. 

In summary, the instantons of Einstein gravity with SU{2) x U{1) isometry re- 
semble the instantons of Horava-Lifshitz gravity that interpolate between the round 
and a fully squashed sphere. The volume of their spherical slices is modulated with 
time. In the non-relativistic case the extrinsic curvature vanishes at the end-points of 
space-time, since these are the fixed points of the Ricci-Cotton flow, and the Gibbons- 
Hawking term becomes obsolete. Nevertheless, it will be interesting to investigate the 
form of the boundary terms in Horava-Lifshitz gravity in more general situations. 
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5.2 Index of Dirac operator for Taub-NUT metric 



The Taub-NUT metric is conveniently written in terms of a radial coordinate r G 
[m, oo) (m is the nut parameter) as 



ar + [r 

r — m 



m 



2\2 



4m" 



(r + m, 



[a 



3\1 



(5.11) 



The end-points of r determine the two boundaries of Taub-NUT space: the removable 
nut singularity and the asymptotic infinity Proper time can be alternatively used, 
setting 

I r -\- m r i—^ 

(5.12) 



, X + m 
m \l x-m 



m arcosh h 

m 



and it runs from to ±oo (the sign depends on whether one imposes self-duality or 
anti-self-duality condition on the Riemann tensor). As for the three Euler angles that 
are used to represent a\ they range asO<0<7r, 0<^< 2n, < xp < An, 
and, therefore, the slices that arise at fixed r have the topology of (xp is extended to 
the double covering of the rotation group SO (3)). Thus, integrating the volume form 
Acr'^ Acr^ = sinO dO Adcp A dip yields a factor 16n^. 

The computation of the index of the Dirac operator D on the Taub-NUT space 
relies on the general formula that was obtained in section 4 without reference to any 
equations of motion for the background geometry. Using equations (|4.56|) and (|4.59|) , 
we have, in particular. 
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(5.13) 



The derivation of equations (|4.56|) and (|4.59|) was done before using the proper (Eu- 
clidean) time t, but it can also be done using the radial coordinate r. The only dif- 
ference is that the end-points of the instanton are now located at r = m and r = oo 
(equivalently at i = and i = cx)) as shown above. 

Here, the slices arising at fixed r are Berger spheres with anisotropy parameter 
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(r + m)2 

and the //-invariant takes the following form, according to formula (|D.42|) , 

(r + 3m)^(r 



(5.14) 



m 



VD\r 



(5.15) 



6(r + m)4 

Note that J/D|r=m = 0, since the nut is conformally related to an endowed with the 
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constant curvature metric (i.e., S = 1). We also note that 7/D|r=oo = —1/6, since the 
Taubian infinity is conformally related to the fully squashed (i.e., S = 0). These 
special cases are discussed in detail in Appendix D and they provide the two extreme 
values of ^ in Taub-NUT space. 

We have, in particular, ^ < 1 so that the spherical slices have positive Ricci scalar 
curvature for all r > m. Thus, there is no level crossing from negative to positive 
eigen-values of the three-dimensional Dirac operator as the slices are deforming ge- 
ometrically by varying r and the index vanishes. This is also apparent from equation 
(|5.13|) . Thus, there is no chiral symmetry breaking induced by this gravitational instan- 
ton background reconfirming the results that were originally obtained in the literature 
|23ll241 long time ago. The Hirzebruch signature of Taub-NUT space is eight times the 
index of the Dirac operator - there are no caveats in this case - and it vanishes, as noted 
before. 



5.3 Index of Dirac operator for Eguchi-Hanson metric 

Another example is provided by the Eguchi-Hanson metric which is conveniently 
written using a radial coordinate r G [a, cx)) as 



dsjiTT — . , . + 



r 



2 



crY + (crY + ( 1 - ^ )(c.3)2 . (5.16) 



The end-points of r determine the two boundaries of Eguchi-Hanson space: the re- 
movable bolt singularity and the asymptotic locally Euclidean infinity. Proper time 
can also be used in this case, setting 

^ / 1 1 3 r(3/4) 

which is given in terms of the hypergeometric function 2F1 and runs from to ±00, 
as for the Taub-NUT space. However, the three Euler angles range differently here, 
as < 9 < n, < (p < 2n, < xp < 2n, and, therefore, they cover only once 
the rotation group S0(3) (S0(3) ~ SU{2) /'Z.2). Letting xp range up to 2n instead of 
4/1 is necessary to remove the bolt singularity from the space-time metric. Then, in 
this case, the slices that arise at fixed r have the topology of the real projective space 
RP-^ ~ S'^/Z2 and angular integration yields a factor 8n^ instead of 1 6/1^. 

Repeating the steps that were taken before, we find that the index of the Dirac 
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operator for the Eguchi-Hanson space takes the form 



The numerical factor of the first term is 1/24 instead of 1/12 that was encountered 
before because of the difference in the angular integration. The slices at constant r are 
Berger spheres with anisotropy parameter 

— = = 1 - 4 ' (5-19) 

which varies monotonically from to 1 as r runs from r = a to 00. 

The j/-invariant of the three-dimensional Dirac operator at any given slice is more 
tricky to find as it is now defined on RP^. One has to use the spectrum of the Dirac 
operator on Berger spheres given in Appendix D and implement the appropriate pro- 
jection by Z2 while computing //d. RP'^ has two inequivalent spin structures provided 
by the embedding of Z2 into SU{2). For the canonical spin structure, the spectrum of 
eigen-values and their multiplicities are given by equations (|D.21|I and (|D.22|) with [,± 
restricted to pairs of positive integers quantum numbers {p,q) with p + q even and 
^0 restricted to even positive integers p. This is so because the states \i,m > used in 
Appendix D.l to solve the eigen-value problem flip sign under rotations by In when 
j assumes half-integer values, i.e., when 2j + 1 is even (the canonical choice). For the 
other spin structure, the eigen-values ^± are restricted to pairs of positive integers 
{p,q) with p + q odd and ^0 are restricted to odd positive integers p so that 2j + 1 is 
odd. Overall, the spectrum of the Dirac operator on S'^ consists of the disjoint union of 
the spectra on RP-^ associated to the two inequivalent spin structures (further details 
as well as generalizations to lens spaces can be found in ||40ll ). 

Using the canonical spin structure on RP^^, we can easily compute j/d from first 
principles at the end points of the Eguchi-Hanson instanton (there is no need to know 
the exact equivariant expression for //d ^or generic values of the anisotropy parameter 
S). At one end, r = a, 'RP^ is fully collapsed to a two-sphere (it is the removable bolt 
singularity of the instanton) and ^ = 0. At the other end, r = 00, the space becomes 
asymptotically locally Euclidean and ^ = 1. The results differ from t] computed at the 
two end points of Taub-NUT space (for which ^ = and 1) because of the Z2 moding 
that affects the spectral asymmetry. 

First, we consider the case ^ = 0. Only the eigen-values ^± with quantum numbers 
p = q = n ^ 'N tend to finite values and yield the spectrum ±n of the Dirac operator 
on (see the special cases discussed in Appendix D.l). They are the limiting finite 
eigen-values on RP-^ endowed with the canonical spin structure, since p + q = 2n is 
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even. The other eigen-values ^± tend ±00. All these eigen-values are equally dis- 
tributed showing no asymmetry. The remaining eigen-values ^0 with even quantum 
numbers p G 2N also become infinite at ^ = 0. These are the relevant ones here as 
they provide the excess modes on the fully squashed KP^ endowed with the canon- 
ical spin structure. Thus, the spectral asymmetry is determined by computing their 
contribution to the //-invariant, 

VDis)\r=a = L'} = ^L 7^ - - 1) ' (5-20) 

p even P n=l ^ ) 

which is expressed in terms of the Riemann zeta-f unction ^(s). Setting s = 0, we 
obtain the //-invariant on the fully squashed RP"^ 

r]B\r=a^m-l)^-\. (5.21) 

This limiting case also shows why the second spin structure on RP^ should be ex- 
cluded from the discussion. All eigen-values ^± and ^0 associated to it tend to ±00 at 
5 — Q leaving no finite modes behind to account for the spectrum of the Dirac opera- 
tor on S^. This is not surprising in retrospect as admits only one spin structure that 
follows from the canonical spin structure on RP^; it also explains the choice of spin 
structure that was made on the slices of Eguchi-Hanson metric for all other values of 
5 that connect continuously to 0. 

Next, we consider the case 5 — 1. The spectrum turns out to be ±(2n ± l)/4 with 
multiplicities n(n ± 1) for all even positive integers n (this follows easily by projecting 
the results found in Appendix D.l to RP"^ endowed with the canonical spin structure 
and the canonical round metric). This projection creates an asymmetry in the spec- 
trum compared to the round S^. The allowed eigen-values are now restricted to the 
subset • • • , -11/4, -7/4, -3/4, 5/4, 9/4, 13/4, • • • instead of the more extended list 
of eigen-values • • • , —7/ 4, —5/4, —3/ 4, 3/ 4, 5/ 4, 7/4, •• • that arise without the pro- 
jection. Then, taking into account the multiplicities of these eigen-values, we obtain 
the corresponding spectral Riemann zeta-function 

1^(s)\r=^ = E - E = i - 2) - W^)) (5.22) 

n even \^"' ~ -^J n even \^"^ ^) ^ 

where j6(s) is the so called Dirichlet (or Catalan) beta function, which is defined by the 
following infinite sum over all integers n > 

B(s) = y ^} . (5.23) 



76 



Then, the ^/-invariant at infinity follows by setting s = and it is found to be 



VD\r=oo = \ (/3(-2) - ^(0)) = 1 (E2 - Eo) = , (5.24) 

where jS(— fc) = are given by the Euler numbers Ej^ for all integers k > 0; we 

have, in particular, Eq = 1 = — E2. 

Combining these results, we find that the index of the Dirac operator for the Eguchi- 
Hanson gravitational instanton background is zero. This follows from (|5.18|) , term by 
term, as 

Ind(D) = 



24r« 



1 / 1 1\ 1 1 

2 -4 + 3 =^-^ = °- 



Thus, chiral symmetry breaking is not allowed to occur in this case either, reconfirm- 
ing the results derived in the literature Ii281 (but see also the review [i29il ). 

An alternative derivation of the //-invariant on RP"^ endowed with the canonical 
metric is provided in the literature by employing the G-index theorem (for a physicists 
description see, for instance, 11291 and references therein). This method was originally 
used in [28] to compute the index of the Dirac operator on Eguchi-Hanson space (but 
see also [1261 for the computation of boundary terms in all A-D-E gravitational instan- 
ton backgrounds by similar methods). Our derivation here is different, as it relies on 
spectral methods, and we have included it since it is not discussed that way elsewhere 
to the best of our knowledge. 



5.4 Miscellaneous remarks 

The index of the four-dimensional Dirac operator on gravitational instanton back- 
grounds M4 with boundaries ought to be zero as simple consequence of Lichnerow- 



icz's theorem | IT7| . Using (ij^Dpi) = —D^D^, which is the specialization of formula 
(|B.9|) to four-metrics with zero Ricci scalar curvature, one immediately sees that the 
zero modes of the Dirac operator should necessarily be covariantly constant spinors. 
Such solutions exist in general, but they can not be normalizable on spaces that extend 
to infinity. Thus, the index theorem that counts the difference between the positive 
and negative chirality normalizable zero modes (with respect to the Atiyah-Patodi- 
Singer boundary conditions) vanishes, as noted in [181. On the other hand, if M4 is 
compact without boundaries, there exist covariantly constant spinors and the index 
is non-zero; in particular, one finds that there are two covariantly constant modes of 
one helicity and no modes of the opposite helicity on K3. Normalizability of the zero 
modes is not an issue in the latter case because of the compactness of space. 

The calculations presented in this section provide a complementary viewpoint 
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based on the ADM decomposition of the four-metric that is most appropriate for com- 
parison with Lifshitz theories (there is no simple analogue of Lichnerowicz's formula 
for the Dirac-Lifshitz operator that can lead to similar arguments for the index). In the 
relativistic case things work out in such a way so that the three-dimensional sectional 
Ricci scalar curvature on E3 remains positive definite as one transverses M4 from one 
end to the other (0 < (5 < 1 in both examples consider above). Then, by Lichnerow- 
icz's formula, which is now applied to the three-dimensional Dirac operator, there can 
be no zero modes (harmonic spinors) on E3. This, in turn, implies the absence of level 
crossing under the spectral flow induced by the deformation of the three-metric along 
the radial direction r and the index of the four-dimensional Dirac operator vanishes, 
as required on general grounds, but through the application of a different argument 
akin to the canonical ADM formalism. 

A novel possibility arises when the Taub-NUT space is taken to another region 
—m<r<m. In this case, one has to flip the overall sign of the metric to main- 
tain Euclidean signature while going beyond the coordinate singularity r = m. Then, 
the anisotropy parameter S = 2m/ {r + m) is allowed to become arbitrarily large as r 
comes closer to the curvature singularity located at r = —m. In fact, the Ricci scalar 
curvature of the spherical slices becomes zero at r = and then turns negative. Har- 
monic spinors first arise at r = —mil, where ^ = 4, and level crossing becomes 
possible for all —m < r < —m/2 according to the analysis given in Appendix D. This 
possibility was recognized in literature quite early | j24l . It seems interesting as it al- 
lows for non-vanishing index of the Dirac operator. The computation proceeds as for 
r > m with the difference that the ^/-invariant has to be replaced with the more gen- 



eral formula (|D.48|) . Thus, by spectral flow, the index equals to the number of modes 
that have crossed zero, S{S). There is a problem, however, in that the space can not be 
truncated to any given value of r other than r = —m for otherwise the metric will be 
incomplete. This is also seen pictorially by considering the point particle interpreta- 
tion of the gravitational instanton in the effective potential shown in Fig.5. As r varies 
from m to —m, the shape modulus /3 varies from to —00. Then, the point particle rolls 
down the infinitely deep cliff of the inverted potential — y(/3) and there is no way to 
stop it before hitting the curvature singularity. In this case, the singularity is inevitable 
for otherwise the metric will be incomplete. Thus, formally, the index of the Dirac op- 
erator in this region of Taub-NUT space is infinite. This example illustrates the fact 
that the index of the Dirac operator is always zero in Einstein gravity provided that no 
curvature singularities are present in the geometry. Certain variations of this theme 
have also been considered in the physics and mathematics literature in recent years 
(see, for instance, | l74li75| V 

Finally, we note that the results obtained for Lifshitz fermions resemble more the 
behavior of relativistic Rarita-Schwinger spin 3/2 fields coupled to gravity rather than 
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spin 1/2 fermions. In that case, the axial current also exhibits an anomaly which is 
given by Tr(R A R), since it is a total derivative term, but its coefficient differs by a 
factor of —21 compared to the spin 1/2 anomalous term l|76l [77l [781; consistency of 
the Rarita-Schwinger equations of motion also requires that the metric background 
should be Ricci flat, which is certainly true for all instanton solutions of Einstein grav- 
ity Then, the integrated form of the anomaly for a relativistic spin 3/2 field reads 

using the corresponding //-invariant at the boundary of space-time, and it provides 
the index of the Rarita-Schwinger operator on M4. Explicit computations have shown 
that 73/2 equals 2t (twice the Hirzebruch signature of M4) on all asymptotically locally 
Euclidean gravitational instanton backgrounds of Einstein gravity | j27ll28ll26l[79| (but 
see also |29J for an overview). This follows by applying the G-index theorem to these 
special spaces and not by computing rj^ ji directly by spectral methods (this compu- 
tation has not been carried out in the literature to the best of our knowledge). For 
instance, one finds J3/2 = 2 for the Eguchi-Hanson space. As a result, chiral symme- 
try breaking can be induced by gravitational instanton effects in supergravity, which 
contains spin 3/2 fields, and it may give rise to helicity changing amplitudes. These 
results should be contrasted to the index of the Dirac operator that vanishes on all 
gravitational instanton backgrounds with boundaries; only on ¥3, which is the unique 
compact gravitational instanton without boundary, one obtains a non-vanishing re- 
sult 7i/2 = t/8 = —2, whereas J3/2 = — 2IJ1/2 = 42, as required by the local form of 
the anomalies. 

It will be interesting to consider the non-relativistic analogue of Rarita-Schwinger 
fields in Hof ava-Lifshitz gravity, examine the structure of its quantum anomalies, and 
also compare them to (super)gravity. More generally, one may consider higher spin 
fermion fields and compute their anomalies. Let A± (M) be the ± chirality bundles on 
the four-dimensional space-time M and let S'^A-t (M) denote their r-fold symmetric 
products. Then, following ||76l |79| (but see also ||29| ), we introduce the general spin 
elliptic complexes D^/2, „/2 : ^mii, w/2(M) l^n/i, m/iiM), where A,„/2, n/iiM) = 
S"^ A+ (M) (g) S" A- (M) and find that their index is related to the spin 1 /2 index by 

Im/2,n/2{M) = —[n{n + 2){3n^ + 6n-U)-m(m + 2){3m^ + 6m-U)]- 

(m + l)(n + l)Ji/2(M) . (5.27) 

It is natural to expect that this general relation extends to non-relativistic theories too, 
since I1/2 = h/2, is the same for both Dirac and Lifshitz fermions. 
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6 Conclusions and discussion 



We have established the universal form of the axial anomaly for massless fermions 
coupled to background gauge and gravitational fields in four space-time dimensions. 
The result is identical for Lifshitz fermions with anisotropic scaling z = 3 as it is 
for Dirac fermions for which z = 1. In either case, the gauge field contribution to the 
anomaly is proportional to Tr(F A F) and the gravitational contribution is proportional 
to Tr{R A R) . This result is not surprising in retrospect, since the anomaly obstructs the 
axial current conservation law and, as such, it should be expressed by a topological 
density that is locally written as total divergence. The overall normalization of the 
anomaly also turned out to be the same for both Lifshitz and Dirac fermions, as it was 
verified explicitly by detailed computation based on Fujikawa's path integral method. 
This was also explained on general mathematical grounds based on the integrated 
form of the gravitational anomaly while showing that the ^/-invariant of the three- 
dimensional Dirac and Lifshitz operators, which are closely related to the value of the 
Chern-Simons action, are in fact equal to each other. Thus, the index of the Dirac- 
Lifshitz operator coincides with the index of the relativistic Dirac operator. The index 
is non-zero on space-times of the form IR x E3 provided that the three-dimensional 
slices E3 can become negatively curved to allow for level crossing. These results are 
of general value and they may be used in various applications in future work. 

Let us further remark that the computation of the anomaly can be extended to more 
general situations, where the fermionic Lifshitz operator is modified by a relevant 
term of the form M^Y^i- This additional term coincides with the spatial part of the 
conventional Dirac operator and influences the transition of the fermion system from 
the z = 3 UV fixed point to the z = 1 Dirac fermion theory (coupled to the same 
background). The anomaly is not affected by the flow, since it is attributed to IR rather 
than UV effects, which explains the origin of our results. Furthermore, we note in 
this context that we have chosen the fermions to be spinors under the spatial rotation 
group, since in the IR we want to flow to ordinary Dirac fermions. 

Among other issues that have been raised in the main text, the most pressing one 
is the derivation of the gravitational anomaly of Lifshitz fermions by methods of su- 
persymmetric quantum mechanics that are analogous to the relativistic case IITOl . This 
will provide a more efficient and systematic way to compute the anomaly and may 
also lead to new developments in the area of non-relativistic supersymmetric field 
theories. A closely related problem, which is connected to our computation of the t]- 
invariant of the three-dimensional Lifshitz operator, is the derivation of the effective 
action induced by massless fermions in three dimensions. It is well known that the 
Chern-Simons action is induced by fermions coupled to an external gauge field |80| 
and this generalizes naturally to the case of a gravitational background field in three 
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dimensions IISTI , in agreement with general index-theoretic expectations. It was also 
pointed out in [82| that the results generalize to higher odd dimensional spaces by 
showing, in particular, that the imaginary part of the effective action is identified with 
//d (up to a factor of nil). Similar arguments apply to the spin 3/2 fields that appear 
in supergravity Since we have shown that //o = >/D/ it is natural to rephrase our result 
by saying that the parity violating part of the effective action induced by z = 3 Lif shitz 
fermions in the background of external fields in three dimensions is also given by the 
Chern-Simons action (with the same overall normalization) after integration of the 
fermions. It will be interesting to provide a diagrammatic computation of the three- 
dimensional effective action in this case (as for the axial anomaly in four dimensions) 
in order to obtain an independent field theoretic verification of this statement. Simi- 
lar inquiries can be made about spin 3/2 fields in non-relativistic theories of gravity, 
closing the general circle of investigations that have been put forward in this paper. 

We have subsequently coupled Lifshitz fermions to (3 + 1) -dimensional Horava- 
Lifshitz gravity and studied the role of axial anomaly in the background of instanton 
solutions. Particular attention was paid to instantons with 511(2) isometry group, 
since there is complete classification and there are explicit forms available for their 
metrics. These instanton backgrounds provide the simplest (yet non-trivial) exam- 
ples of more general solutions and they also have analogues in Einstein gravity to 
which direct comparison can be made. The main difference with ordinary instan- 
tons is their chiral nature which is inherited from the parity non-invariance of three- 
dimensional topologically massive gravity. The integrated form of the anomaly on 
such backgrounds, including all boundary effects, provides the index of the fermion 
operator, which can be used to study chiral symmetry breaking induced by instan- 
tons in gravitational theories. In particular, it was found that the index can become 
non-zero on certain instanton backgrounds of a unimodular phase of Hof ava-Lifshitz 
gravity provided that the coupling of the Cotton tensor term exceeds a critical value. 
The result was derived by spectral flow methods and it relies heavily on the existence 
of harmonic spinors on sufficiently deformed three-spheres, which are the leaves of 
space-time foliation. Otherwise, the geometry of the background does not allow for 
bound states of the four-dimensional fermion operator, since the index is zero, ex- 
actly as in the case of instanton solutions of ordinary Einstein gravity with bound- 
aries. Then, in the unimodular phase of Horava-Lifshitz gravity, it is possible to vio- 
late chiral symmetry by non-perturbative instanton effects, provided that the volume 
of space exhibits a lower bound that it is determined by the ratio of the Ricci to Cot- 
ton tensor couplings raised to the third power. This, in turn, can lead to baryon and 
lepton number violation processes in the deep ultra-violet regime of the theory. It is a 
rather striking phenomenon that differentiates Horava-Lifshitz gravity from ordinary 
Einstein gravity once more. Certainly, further investigation is required in order to es- 
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tablish this result in general, beyond the class of instantons with 51/(2) symmetry, and 
investigate its phenomenological and cosmological consequences. 

A variant of this theme can be put forward by considering spherical slices of the 
form /Y, instead of itself, where F is a finite subgroup of SO (4) acting freely 
on S^; these include lens spaces and certain other generalizations by crystallographic 
groups and they can arise as prime factors in the decomposition of more general com- 
pact three-manifolds. First of all, it will be interesting to consider the normalized 
Ricci-Cotton flow on such spaces and establish an analogue of Hamilton's space-form 
theorem (the latter states that if g(0) is a metric of positive Ricci curvature on a three- 
manifold Z3, then the volume-normalized Ricci flow exists for all time and converges 
to the round metric on S'^/F for appropriately chosen F C S0(4) [|83|, but see also 
|84| for a more pedestrian account of this and more recent results on the subject). It 
will also be interesting to consider eternal solutions of the flow equations that will 
provide the instantons of Horava-Lifshitz gravity in this more general set up. Such 
generalized configurations, if they exist, will qualify as non-relativistic analogues of 
the general class of asymptotically locally Euclidean instantons of Einstein gravity 



that were proposed by Hitchin | [85| and constructed by Kronheimer |[86|. These met- 
rics, however, do not admit any isometries, in general, and the same thing should be 
expected to happen here. Application of the index theorem in these backgrounds will 
then require an equivariant generalization of the ^/-invariant based on the spectrum of 
the Dirac operator on lens spaces [40] (and generalizations thereof), as for the Eguchi- 
Hanson metric in Einstein gravity where Z3 ~ S"^ / Z2 (see, for instance, for an up 
to date exposition of this subject with several references to the mathematics literature). 
The task seems formidable. 

The use of spherical spaces S^/F in the canonical (3 + 1) -dimensional formulation 
of Einstein gravity is quite standard by now and it is often associated to the concept of 
geons. Thus, what we are proposing here is to consider the non-relativistic analogue of 
geons in the context of Hofava-Lifshitz gravity. Such configurations are rather exotic 
in that they admit asymptotically trivial diffeomorphisms which are not deformable 
to the identity and they can act non-trivially on the quantum state space (for example, 
they can give rise to states with half-integral angular momentum); we refer the in- 
terested reader to the papers Il88ll89l for further details. The homotopy groups of the 
diffeomorphism group of these spaces are directly related to the homotopy group of 
their isometries, noting, in particular, that 7ro[Diff (S'^/F)] ~ no[N{T) /T], where N(F) 
is the normalizer of F (the isometries of S'^/F are those elements h G SO (4) for which 
the image of an orbit of F is another orbit of F, i.e., when h is in the normalizer N(F) 
of F defined by hTh~^ = F). The disconnected components of the diffeomorphism 
group have been calculated and tabulated in [|89| for all different type of compact 
spherical spaces. The large diffeomorphisms of E3 can lift to the foliation preserving 
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diffeomorphisms of R x E3. It will be interesting to study their implications, if any, to 
the gravitational anomalies in 3 + 1 space-time dimensions and understand if Weyl 
fermions can be consistently coupled to geometry in the quantum theory We have 
not yet delved into the details of all these generalizations, but we hope to return else- 
where and report some results on the fate of chiral symmetry breaking and its variants 
on geonic backgrounds. 

Another interesting generalization that is left open to future work is the local 
form of the axial anomaly inld + 2 space-time dimensions for Lifshitz fermions with 
anisotropic scaling z = 2n + 1. In this case, the Dirac-Lifshitz operator is naturally 
defined, modulo factor ordering ambiguities, as z^^'D^, = iy^Do + Z7'(V^)"D; for all 
integer values of n (ordinary Dirac fermions correspond to n = 0). It is practically 
useful to consider only the values n = 0, 1, 2, ■ ■ ■ ,d focusing, in particular, to the 
"maximal" case n = d. Then, in this case, it is natural to expect that the axial anomaly 
will be the same as for Dirac fermions coupled to background gauge and metric fields 
(see, for instance, IITOll for the relevant results in relativistic higher dimensional field 
theories). One should be able to explicitly verify this statement by path integral meth- 
ods or even more easily by developing alternative computational methods based on 
supersymmetric quantum mechanics, which unfortunately are lacking at the moment. 
Coupling the Fermi theory to Horava-Lifshitz gravity will then allow to compute the 
index of the fermion operator on the corresponding instanton backgrounds (yet to be 
found) using the spectral theory of the Dirac operator on odd dimensional spheres. We 
refer the interested reader to the second reference in ||40| , which proves the existence 
of metrics admitting harmonic spinors in dimensions 2d + 1 = 3 mod 4, thus gen- 
eralizing Hitchin's results for three-dimensional Berger spheres |j39l; curiously, these 
more general results single out 2d + 2 = 4k space-time dimensions which also have 
prominent role in relativistic theories and for which the index can be non-vanishing. 
Higher dimensional analogues of Horava-Lifshitz gravity together with the associ- 
ated instanton solutions have not been studied so far in the literature, but an interest- 
ing case should be made in 2rf + 2 dimensions with anisotropic scaling z = 2rf + 1 (c.f. 
n = d). We hope to return to these generalizations elsewhere. 

Persistent boundary effects in gravitational theories arise when the spatial slices E 
have boundaries, as in the case of gauge theories with background monopole charge. 
Appropriate boundary terms should be added to the action, following the canonical 
treatment described in | [90ll9T| for 3 + 1 Einstein gravity. Carrying the analysis over to 
Horava-Lifshitz gravity should be rather straightforward and would make our anal- 
ysis of the Euclidean action bounds applicable to more general situations. Then, the 
computation of the index would require the open space generalization of the index 
theorem developed in |,65,|. 
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Another open question is the fate of Weyl symmetry upon quantization. It primar- 
ily concerns the anisotropic Weyl invariant phase of Hof ava-Lifshitz gravity, which is 
thought to govern the deep ultra-violate regime of the theory, but this should break as 
one flows to the infra-red by suitable running of the parameter A away from the spe- 
cial value 1/3. Although this issue remains largely unexplored in the literature, it is of 
central importance for future developments in the subject. It also relates to chiral sym- 
metry violation by gravitational instanton effects, since this possibility seems to arise 
only in the unimodular phase of the theory (this is a weaker form of Weyl invariant 
phase that also arises at A = 1/3). A simpler question concerns the Weyl anoma- 
lies of matter fields in the background of gravity, such as scalar, fermion and vector 
fields, which provide the quantum obstruction to the tracelessness of their energy- 
momentum tensor. The Weyl anomaly only arises in even dimensions, like the axial 
anomaly, and it has been exhaustively analyzed for relativistic field theories (for an 
overview and history of the subject see, for instance, [921 and references therein). The 
form of this anomaly remains largely unexplored for non-relativistic field theories, 
apart from the case of conformally coupled Lifshitz scalars that provide the simplest 
example of anisotropic Weyl invariant models at the classical level [ 931 1 (but see also 
Il94| ); even in this simple case there are some questions that still remain unanswered 
and call for further work. One should also be able to compute the Weyl anomaly for 
other fields, such as Lifshitz fermions, and compare the results to the relativistic the- 
ories. The coupling of the fields to geometry should be made conformal, as in [93J. 
In general, there is no a priori reason to expect that the coefficients of the anomalous 
terms will stay the same (in this respect the computation of the axial anomaly is much 
cleaner because its form is protected by topology). Finally, the breaking of Weyl sym- 
metry in the gravitational sector should be addressed by appropriate methods. 

At this end, it also seems appropriate to view our results as being part of a more 
general relation between non-relativistic fermions in condensed matter physics sys- 
tems and Dirac fermions in relativistic field theories, as they exhibit similar behav- 
ior. Other important examples of this kind (though seemingly unrelated at first sight 
to the problems we have considered here) include fermion number fractionization 
in polyacetylene (see, for instance, ll95B and references therein) and more recently in 
graphene Il96l , which are all accounted by the presence of fermion zero modes in non- 
trivial backgrounds. We should delve deeper into these problems and understand the 
fundamental reasons for this interplay. Definitely, Lifshitz type systems provide an 
interesting class of models to further this endeavor. 
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A Dirac matrices and their trace identities 



Recall that the Dirac gamma-matrices in the flat four-dimensional Euclidean space 
R"* with metric S^iv satisfy the anti-commutation relations 



[7^ Y]+=2S^\ 



(A.1) 



which follow from their Lorentzian counterpart by analytic continuation t — > it and 
7° ij^. Then, the chiral representation of the Dirac algebra is provided by the 4x4 
Hermitian matrices 



(A.2) 
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(A.3) 



and 75 = —jOjiji^s jg represented by the Hermitian matrix 



/ 1 



75 



(A.4) 



-1 



that anti-commutes with all jf^ and squares to 1. 

The trace of any product of odd number of gamma-matrices vanishes identically. 
Furthermore, we have 

Tr75 = = Tr(757^T^) , (A.5) 



whereas 



Tr(757^7^7''7^) = -^e^^^^ 



(A.6) 



written in terms of the fully anti-symmetric Levi-Civita symbol with e^^^ — 1. We 
also have Tr (7''7^) = 43^^^ and the useful relation 

^li^j" = S^'j" - S^^Y + ^""7^ + e^^^M^ls (A.7) 
from which we obtain the following trace identity among gamma-matrices, 

Tr(7?'7^7'^7'^) = 4{S^''S''^ - Sf'^S"^ + Sf^S""") (A.8) 
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and, consequently, the identity 



^jiv^KXpa _ ^jiK^vXpa _|_ ^vK^jiXpa _|_ 



^\p^}ivKa _ ^Xa^pivKp _|_ ^pa^jivKX 



(A.9) 



Finally, the trace of 75 with the product of eight gamma-matrices reduces to a sum 
of traces of 75 with the product of six gamma-matrices giving, in particular, the trace 
identity 

Tr (75 l''']\ba, 1,\\W, Y]\bc, 7i]) = \ea''''{si6,, - S^S,,) + 

le/^'''{^l'hd - ^t'^ad) + le/''''{^ihc - ^ihd) + lec'^'^'^'isiSad - ^i^u) + 
\e^''''''\di5ac - did,,) + \e,'''''\5X5ac - ^Ubc) + ^e/^''' {^'c ^ad - Si^ac) + 
\ed''''{5ihc - ^i^ac) + \eM''\5i5i - 5(5^) + - S^S^') + 

3^ c'd'i^a'^h' c'd'(^a'^h' a^'a^'S -U "^c- a'h'(^c'^d' xc'xd'\: 

Ij^ad \.0c -^b^c ) + 2^bc [i>d -^a^d) + 2^ad l^c ^b -^b^c ) + 

^ej'\5i5i - Sisf) + 4e,/''{K'4' " K'^a) + ea/^'sfsj' + eb/^'s^S',' + 

„ a'd'^b'^c' , ^ h'c'^a' ^d' , „ h' d' ^a' ^c' , „ a'c'^b'^d' , „ ^'c'xa'^d' , 
^ad h + ^hc t>fl '>d + ^ac h ^d + ^bd + ^ad h "c + 

^fec <>a <^d ~ ^ac (>b (>d ~ ^ac (>b ^d ~ ^hd <^a ~ ^bd <>a ~ 

ead''''5i5i - Cbc^'^'sisi - eJ^'sUi " eb/^'s^a'^^ (A.IO) 

This expression is manifestly invariant with respect to the simultaneous exchange of 
indices a -n- c, b d and a' -n- c' , b' -n- d' . It is written directly in a form that will 
be used in the computations presented in section 2.3, where the gamma-matrices are 
labeled by flat (tangent) Euclidean space-time indices, as for R*. The expression given 
by (|A.10|) has been put into final form using the relation 

^abc^^a'b'c'd ^ gaa' ^gbb' gcc' _ ^bc'^cb'^ _ ^ab' ^^ba' ^cc' _ gbc' ^ca' ^ gac' ^gba' gcb' _ gbb' gca' _ 

(A.ll) 

These are the trace identities that will be needed for the computation of the axial 
anomalies of fermions coupled to background gauge and metric fields. 
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B Some geometrical apparatus 



We review some disperse geometrical notions and formulae that are relevant for the 
coupling of fermions to background metric fields and the computation of the axial 
anomaly in the presence of gravity. It also sets the notation used in the text. 

B.l Spin connection and Lichnerowicz formula 

The coupling of fermions to a background metric field G^v is achieved using Cartan's 
formalism (see, for instance, [29]). To set up the notation, we introduce vierbeins 
as G^v = Sah^^fiS^v, where S^,}, provides the flat tangent space-time metric in the 
Euclidean domain. We also introduce the inverse vierbeins E/^ as G^^ = S^^Ej^Ei,^, 
which satisfy the relations Ea^ = 5ahG^''^e^y and Ej'^e^^ = S^. The components of the 
spin connection co"^ are 

= e% {d^E,- + Tl^E/) = -E,- (d^e% - rj,e%) (B.l) 

so that V^e^v = V^e^v + <^^b}i^^v = ^ii^^v — T^v^^A + "^"fe^^^v = 0. Cartan's equations 
(in the case of no torsion) provide the components of Riemann curvature tensor as 

d^cv^K - d,cv'\ + oo\^oo'\ - oo\yio'\ = R'\y = e\e\R^\, . (B.2) 

The Dirac operator i^y^D^i coupled to the metric is defined using the gamma-matrices 

J^^=YEa^, [Y,l']+ = 25a,, (B.3) 

which are expressed in terms of their tangent space-time counter-parts 7" and the 
covariant derivative operator acting on spinors 

which is defined through the spin connection. Then, we have 

[D,„ Dy] = i[7„ ^bW^v (B.5) 

in accordance to the fact that the metric analogue of the gauge field is provided by 
(z'/8)[7fl, jifjcoii"^ in the spinorial representation of the tangent space-time rotation 
group and its field strength is the Riemann curvature. More generally we define 

D^ = Vf, + h7a,7bW^ (B.6) 
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with the Christoffel affinity included, which also satisfies equation (|B.5|) and has the 
property that [D^, 7^] = 7" V^^E/ = 0. 

Using (|B.6|) , we also have the following relation for the square of the Dirac operator, 

ihW,)^ = -D,D^ - 1[7«, 7fc] be, Id] R"''' . (B.7) 
Contracting jalbleld with the identity R"^"^ + R""^^ + R"'^^" = 0, it follows that 

[7a, 7b] be, 7d] R"'"' = ^7a7b7e7d R'""" = -^7a7b R'" = SR , (B.8) 

which is solely expressed in terms of the scalar Ricci curvature. Then, one arrives at 
the celebrated Lichnerowicz formula for the square of the Dirac operator written in 
terms of the Laplacian operator acting on spinors, |[T7|, 

(Z7^'D^)2 = -D,,Df' + ^R . (B.9) 

is often called Bochner Laplacian to distinguish it from the Laplace-Beltrami oper- 
ator acting on scalars. 



B.2 Geodesic interval and Synge-DeWitt tensors 

Another item that is needed for the computations (in Appendix C and in the main 
text) is the geodesic interval [41], which is defined as one half the square of the distance 
along the geodesic between any two points x and x', 




(B.IO) 



This notion also appeared in the literature under the name "world function" (in Synge's 
textbook 1421) and plays important role in the theory of Green's functions on Rieman- 
nian manifolds with metric G^y, including manifolds of indefinite metric 143] . 

The geodesic interval is a symmetric function a{x, x') of x and x' that transforms 
as biscalar, i.e., as a scalar separately at x and x'. It satisfies the differential equation 

cr(x,xO = = ^(V;^)(V''^(r) (B.ll) 

with boundary conditions 

a(x,x) = 0, lim VuCix,x') = = lim V'^cr(x,x') (B.12) 

x^x' x^x' ^ 
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and it obeys the important relation 



lim VuVv(r{x,x') = - lim VuVXx,xO = Guv ■ (B.13) 

x^x' x^x' 

As such, it provides a natural generalization of cr{x,x') = {x — x')^ /2 from flat to 
curved space. In a general Riemannian manifold cr(x, x') is not single-valued, except 
when X and x' are sufficiently close to each other. 

The geodesic interval a is a function that can be defined through the first order 
non-linear equation (|B.11|) . This has the interpretation of the Hamilton-Jacobi equa- 
tion of a mechanical system whose evolution is described by the geodesic equation on 
a given manifold. It allows to express the invariant delta-function on a Riemannian 
manifold (used in the computation of Green's functions) as 

-J{x - x') = [ ^,'*^,<vv(xy) ^ -^4) 



VditG Vd^J (27r) 

which, in turn, shows that exp{ikiiV^(r{x, x')) provide the analogue of plane waves in 
curved space. In flat Euclidean space, we have x') = x^ — x'^ and one recovers 

the ordinary plane waves, which depend linearly on the Cartesian coordinates in the 
exponent. In general, cr can not be determined in closed form, but in the physical 
applications it is sufficient to know its expansion as x and x' come close to each other. 

The coincidence limit x — > x' of succussive covariant derivatives of the geodesic in- 
terval, viz. Vv • • • VkCT, are called Synge-DeWitt tensors and they are of the utmost 
importance in computing asymptotic expansions of Green functions of partial differ- 
ential operators. They are also very important for the computation of the gravitational 
contribution to the axial anomaly. Adopting the short-hand notation 

[V„ Vv • • • VkCt] = lim V. V, • • • VKcr(x, x') , (B.15) 

x^x' 

which is commonly used in the literature, we immediately have following results for 
the lowest rank Synge-DeWitt tensors, restating (|B.12|) and (|B.13)) , 

[cr] = , [V^cr] = , [V^Vya] = G^y . (B.16) 

For the purposes of the present work we also need the third and fourth rank tensors, 
which turn out to be (in the absence of torsion) 

[V„Vi,VK£r] = (B.17) 
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and 

[V^Vi^VkVa^t] = - (J^^kAv + RfiMv) (B.18) 

written in terms of the Riemann curvature tensor. All higher rank Synge-DeWitt ten- 
sors involve products and derivatives of the curvature tensor but their expressions are 
quite cumbersome to present here (the sixth rank tensor involves about forty different 
terms and the situation becomes increasingly more complex for tensors of even higher 
rank); fortunately, we do not need to go beyond rank four. 

The computations are based on equation (|B.11|I by taking successive derivatives 
and applying the Ricci and Bianchi identities as many times as necessary, depending 
on the rank of the tensor. In this fashion one obtains recursive relations for the Synge- 



DeWitt tensors that can be iterated starting from the basic ones (|B.16|) 
For example, differentiating (|B.11|I three times yields 



V^VvVkCt = (V;,Vt,VV)(VKVAC7) + (Vt,VV)(V^VKVAcr) + 

(V;,VV)(V.V,VAcr) + (VV)(V^V.VKVAcr) (B.19) 

that results into the following equation (term-by-term) in the coincidence limit, 

[V^^VvVkCt] = [V^VyVM + [V^VkV^^ct] + [VvVkV^ct] . (B.20) 

Simple rearrangement yields the equivalent expression 

- 2[V,,V^VkC7] = ([V^VkVvCT] - [Vf,V^.VKC7]) + ([V^.VKV,,(7] - [VyV^.V^a]) + 

([VvV^VkCt] - [V^VyV^cr]) . (B.21) 

The first pair of terms on the right-hand side gives zero because tr is a scalar func- 
tion of X and [Vk, Vv]ct = 0. The second pair of terms also gives zero because 
[Vk, V^i]cr = 0. Finally, the third pair of terms equal [[Vy, V^,](Vk(^)] = [R^kiiv{^ \^)] 
and vanishes because [VaC^] = 0. This proves the relation (|B.17|) . The proof of relation 
(|B.18|) proceeds along similar lines, but it is more lengthy to present it here in detail 
(note, however, that our sign convention for the Riemann curvature tensor is opposite 
to the one used in |41|1). 

We also refer the interested reader to the textbook ||97| (see, in particular, chapter 
8) for a comprehensive account of the recursive calculations of Synge-DeWitt tensors 
and their use in the asymptotic expansion of Green's functions on curved space-time. 
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C Axial anomalies for relativistic f ermions 



In this appendix we summarize the computation of the axial anomaly for a mass- 
less Dirac fermion coupled to a background gauge and metric field, using the elegant 
method of Fujikawa |[5H6l in Euclidean space. Thus, we will consider the regularized 
quantity 



A{t,x) = lim 



lim 



n 



(C.l) 



defined in terms of the eigen-values and orthonormal eigen-functions of the interact- 
ing Dirac operator, (z7^'D^)^„ = A„^„, and evaluate it. The axial anomaly simply 
reads V^Jg = 2A. An alternative (closely related) method of computing the anomaly 
is provided by the point-split method of the fermion bilinear in the axial current (see, 
for instance, 1*91 and references therein). 

Our purpose is to set up the notation and present the main formulae and tricks 
needed for the computation of the gauge and gravitational contributions to the axial 
anomaly. It also simplifies the presentation in the main text of the paper without 
making detours. Then, the intermediate steps and results of the relativistic theory can 
be easily compared to the corresponding contributions received by the divergence of 
the axial current in the non-relativistic theory of Lifshitz fermions with anisotropic 
scaling z = 3 in four space-time dimensions. 



C.l Gauge field contribution to the anomaly 

Minimal coupling to an Abelian or non-Abelian gauge field in flat space-time 
amounts to considering the covariant derivative = ^^ — iA^ so that [D^, Dy] = 
—iF^y is expressed in terms of the corresponding field strength. Then, we have the 
following relation, 

{iYD,f = -D,D^^ + ^[7^ Y] F,y . (C.2) 

To evaluate the anomaly it is convenient to use the eigen-values and functions of the 
free Dirac operator i^^ii' which are plane waves. This will enable to extract the gauge 
field dependence of A{t, x) using a basis that has no gauge field dependence by itself. 
Then, in terms of this basis, the primitive form of the local anomaly takes the form 

AU,x) = lim Tr / (C.3) 

^ ^ M^co J (27r)4 ^ ' 
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taking the trace on anything available on the right-hand side. Combining this with 
relation (|C.2|) yields immediately 



r d^lc r 1 / 7 



In the last step has been rescaled to M/c^. 

It is now clear, by expanding the exponential in power series of 1 / M, that only 
terms up to order will survive in the limit M — ^ oo. Furthermore, by the 

trace identities of the Dirac gamma-matrices with 75, a non-vanishing contribution 
requires having at least four such gamma-matrices and, therefore, the only term that 
survives at the end is 

Mt>^) = ^ / ^"'^''"Tr (75 |[7^ 7l |[7^ 1%.F,A^ ■ (C.4) 

Recall at this point the identity Tr(757^'7^7'*'7'^) = —/^el^^^^ for Euclidean gamma- 
matrices, which takes care of their trace. Also, the integral over the Euclidean space 
momenta k is easily performed and yields immediately 

e ^f^' = —— / dk k^e " = (C.5) 



(27r)4 (27r)4 }q 16n^ 

using spherical coordinates and vol(S^) = 27r^ for the unit sphere. 

Putting all together, the final result for the gauge field contribution to the axial 
anomaly, which equals 2A, is 



^^^'^ = -T^^''''Tr(f,,f,,) ^ -^Tr(f A f ) . (C.6) 



C.2 Metric field contribution to the anomaly 

Next, we summarize the computation of the gravitational contribution to the axial 
anomaly for Dirac fermions, following (SUll. The primitive form of the anomaly is 

A{t,x) = lim X](pJ(i,x)75e-('^"°"/^)>n(i,x) 

= lim lim Tr / 75 ,-('7^'D„/M)2g*„V/'.(xy) (^.7) 

M^oox^x' J (27r)4 ^ ' 
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by passing to the plane wave analogue basis that naturally involves the geodesic in- 
terval^. Here, the derivatives of the Dirac operator are taken with respect to x (and 
not x'), whereas the limit x — > x' should be taken after computing the action of the 
operators on exp(z7c^ W(x, x')). Then, taking into account Lichnerowicz's formula 
(|B.9|I , we obtain the following intermediate result for A{t,x), 

The first line follows by acting with the operator (27^ D^)^, in which case is simply 
replaced by + ils.^ with A^i(x, x') = kyV (j{x, x') (the curved space analogue of 
+ ik^), whereas the second line follows by taking the limit x x' , in which case 
Aj,(x, x') is replaced by k^ thanks to relation (|B.13|) for the rank two Synge-DeWitt 



tensor, and by rescaltng k^ to Mk^. Although the rank three Synge-DeWitt tensor 
(|B.17|) vanishes, so that [V^A^] = 0, all higher rank tensors do not vanish and can in 
principle give non-zero contribution in the coincidence limit x — ?■ x'. Careful treat- 
ment of the problem show that none of these additional terms have sufficient number 
of gamma matrices to survive under the trace, and, therefore, we make no error by 
taking the coincidence limit right from the beginning, as shown above. 

The rest of the computation proceeds as in the gauge field case by expanding the 
exponential in power series of 1 /M and taking the limit M — > 00. The only difference 
is that the curvature term in the exponential contains no gamma-matrices and so the 
non-vanishing contribution to the trace comes from a different combination of terms 
of order l/M^. Dropping all lower order terms, which have zero trace, and all higher 
order terms, which vanish identically as M — > 00, we obtain the following contribution 
to the axial anomaly, 

A{t,x) = TrJ^e-^i'^"j5(^^{D')' + lk^kyk,k;,DW^D''D^ 

-hc^ky {p^D^'D'' + D^D^D" + D^'D'^D^) j , (C.9) 

setting here for convenience = D^D?'. 



may choose to work in the Unearized approximation setting Gjtv — S^v + h^v and assuming 
that hjiv and its derivatives are sufficiently small. Then, exp(zA:^V''£7'(x, x')) ~ exp{ikii{xt' — x''')) and 
equation | |C.3l l follows again by letting x ^ x' after passing exp( — //c^x''') to the left of the operator 
(which only depends on x and not x'), as in the gauge field case. However, the resulting gravitational 
contribution to the anomaly will be only valid at the linear level, as in the original work jT), and its 
non-linear generalization will be left open. Making use of the geodesic interval remedies the situation, 
as it provides the means to compute the anomaly at the full non-linear level. 
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The integrals over the Euclidean space momenta k are easily performed by intro- 
ducing a unit vector with components kfi, so that fc^ = /c/c^ with kpk^ = k . Alter- 
natively, we may choose to work with tangent space-time indices fc^ = e^^ka, since 
kyM^ = kak", and note the identities 

f d^k _fc2,2ff 1. 271^ /■~„,5_fc2 1 , 



and 



/ 



^ e k^ kahkckd = TT^ i^ab^cd + ^achd + ^adhc) / {C.ll) 



(27r)4 '--(^'--a 54^2 

which are manifestly symmetric in all permutations of the indices. These integrals 
are evaluated using spherical coordinates in the space of four-momenta, choosing 
ki = sin^sin^sini/?, ^2 = sin^sin^cosi/?, kj, = sin^cos^ and ^4 = cos0 in terms of 
the corresponding angular variables that range asO<0<7r, 0<(/)<7r and 
< i/7 < 27r. The integration measure over the Euler angles in momentum space 
is swi-dsm(pddd(pd^ so that the volume of the unit sphere is vol(S^) = Iti^. Taking 
these results together with the integral (|C.5|) , one finds that equation (|C.9|) yields in 
steps 

= T^Tr(75[D„D,][D^D^]) 

= ^^T^Tr (75 ba, It] be, ld])R'\vR''^' 

]-^eatcdR'\vR''^' , (C.12) 



using ^"(^ici^lalhlcld) = ~^^ahcd foi" the tangent space-time indices. We may pass to 
the fully covariant totally anti-symmetric symbol defined as c^vkX = s^^n^^ve^^xe^ \£ahcd- 

Thus, the final result for the metric field contribution to the axial anomaly, which 
equals 2A, is 

= -38i^'''''^'"'^^ - -I^Tr(R A R) (C.13) 
written in terms of the corresponding Ricci curvature two-form. 



Here and in the main text e^^^A -^^ defined as e^^^^ = 1 / VdetG and £0123 = V detG, 
unless stated otherwise. 
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D Dirac operator on Berger spheres 



In this appendix we consider the Dirac operator defined on three-dimensional spheres 
that are endowed with SU{2) homogeneous metrics. Following ||39l (but see also 
mni), we solve the eigen-value problem of 27' D, (i = 1, 2, 3) acting on two-component 
spinors. 



(D.l) 



zYD,Y(x) =^Y(x); Y(x) = 

in the special case of Berger spheres that exhibit an additional U{1) isometry. This 
provides a tractable example where all computations can be carried out to the end 
and the results are relevant for evaluating the integrated form of the axial anomaly on 
gravitational tnstanton backgrounds with SU{2) x U{1) isometry group. Along the 
way, we point out the limitations to compute the spectrum in closed form for totally 
anisotropic geometries on S^. 

It turns out that there can be zero modes (often called harmonic spinors 113) whose 
number depends on the metric on S"^ and can grow without bound. This should be 
contrasted to results for the space of harmonic spinors on a two-dimensional compact 
Riemann surface of genus g. In that case, the number of zero modes can not grow very 
large as it is bounded by the topological invariant g + 1; also, on a Riemann surface 
there are several inequivalent spin structures, whereas in the case we are considering 
here there is only one spin structure. The presence of harmonic spinors allows for 
level crossing under the spectral flow and it is reflected in the j/-invariant of the Dirac 
operator. All these problems will be examined in detail in the following. 



D.l Spectrum of Dirac operator 

Let us consider the Bianchi IX class of three-dimensional metrics, as in section 3.3, 

ds^ = ^r{a^f + ^2{(r^f + ^^{cr'f , (D.2) 
using the left-invariant one-forms of the group Sll{2) that satisfy the relations 

d(T^ + -e^]K<^^ A (7^ = (D.3) 

and coefficients 71, 72, 73 that are, in general, unequal. There is a dual basis of vector 
fields fi associated to with < >= 5\ that satisfy the commutation relations 

Lf;, /;] = -eij^fK (D.4) 
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and represent the Killing vector fields of this particular class of geometries. Both 
and fi can be written in terms of three angular coordinates = {6,cp,xp), but the 
explicit expressions will not be needed here. 

We also introduce the orthonormal coframe = y/ji cr^ associated to the dreibeins 
e^i with = e^idx^ and the orthonormal frame Ej = fil a/tI associated to the inverse 
dreibeins Ej' with E/ = E/'B/, which allow us to define the Dirac operator on as 

^VD, = i^'U [b, + 1[7;, 7x]a;^''.) • (D.5) 

The formulae are identical with those found in Appendix Bl. The only difference is 
that the space-time index \i is now replaced by the space index i and the Dirac matrices 
7" are now replaced 7^; the latter satisfy the anti-commutation relations [7^, 7^] = 
7.5^^ (since I and / are tangent space indices) and they should be identified with the 
Pauli matrices (|A.3|) as 7^ = Cx, 7^ = —cfi and 7"^ = 03 (the Pauli matrices should 
not be confused with the one-forms of SLZ(2) used above). Then, ej/xi?^/ 7^]/8 
provide the spinorial representation of the Lie algebra generators // in (|D.4|) . 

Explicit calculation shows that the connection one-forms of Bianchi IX metrics are 

= r, I (-71 -72 + 73)g^/ (D.6) 
2V717273 

= J- (7l-72 + 73)g^ (D.7) 

2^717273 

cv^3 = r. . ^ ^ (71 - 72 - 73)g^ • (D.8) 
2V717273 

Using the fact that the square of the Pauli matrices is the identity matrix, it follows 
that the spin connection term of the Dirac equation (|D.5|) equals to 



Then, the Dirac operator for all homogeneous metrics on S'^ is the following first order 
differential operatoj^ written in terms of the three vector fields 



f 2/3/x/73 ifl/y/ll- fl/y/ll ^ 

V ^7i/\/7r + /2/x/72 -if3/^/l3 ) 



1 

+ —==(71 + 72 + 73) • 
4V717273 

(D.IO) 



^^Comparison with ||39l (see, in particular, p. 28) can be readily made by setting f\ = —£2/2, fi = 
—£3/2, /3 = — ei/2 and 71 = A2/4, 72 = A3/4, 73 = Ai/4 in Hitchin's basis of the Lie algebra 
generators. 



97 



The eigen-values of this operator can only be found explicitly when there is an 
additional U{1) isometry associated to axially symmetric geometries. Letting any two 
coefficients of the metric become equal, say 71 = 72 = 7 (all other choices are mu- 
tually related by a Z3 permutation symmetry of the three principal axes of S'^), we 
obtain 



7 



(D.ll) 



where = 73/7 measures the degree of anisotropy of the associated Berger sphere. 
Since the Dirac operator scales uniformly as y^, we may drop the conformal factor 7 
and only consider the Dirac operator on with line element {c^)^ + (c^^)^ + S^{(r^)^, 



( ifs/S ih-h\ 
\ ifi+fi -ifs/^ J 



(D.12) 



The following simple but very important observation allows to determine the spec- 
trum of the Dirac operator (|D.12[) on Berger spheres. Writing 27' Df = Q + {S^ + 2)/4:S 
and setting f± = ifi =p /2, we note that the operator + Q/S takes the diagonal form 



o 1 





I f+f- 



ifi-ifs)/^' 







(D.13) 



f-f+-ifi + if3)/^' J 



after making essential use of the Lie algebra commutation relations (|D.4|) . The eigen- 
value problem of the resulting second order operators 



J'±=f±fT-^ifi^if3] 



(D.14) 



can be easily solved by noting that they commute with the Laplace operator A act- 
ing on functions on the Berger sphere (actually, 2 A = + J^-). As such, they all 
have a common system of eigen-states provided by the appropriate representation of 
SU{2) associated to If we were using unitary irreducible representations of SU{2) 
with definite spin ; = 0, 1 /2, 1, 3 /2, • • • , the action of the generators on the base states 
>; — j < m < ;} would be /±/=p|;,m >= (; ± m)(; =f m + m > and 
fz\i,m>= im\i,m >, giving 



T±\i,m>= ±m)(; =F m + 1) + ^m(m =F 1) 



I;, m > 



(D.15) 



However, one has to be slightly more careful in applying group theoretical methods 
here as the appropriate representation is not irreducible. 
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The problem at hand is closely related to the determination of the energy levels 
and eigen-states of a three-dimensional rigid rotor with Hamiltonian 

in the special cas^£§ that the moments of inertia are Ii = I2 ^ I3. We simply have 
to set // = iL I to turn (|D.4|I into the commutation relations of the quantum angular 
momenta operators (with Planck's constant normalized to 1). But there is an impor- 
tant difference between the two problems related to the representations of the SU{2) 
symmetry group. The angular momenta Lj are differential operators of two angular 
variables {6, cp) that parametrize in a spherical coordinate system in R'^. As a re- 
sult, the eigen-states of the rigid rotor are spherical harmonics provided by \j,m > 
with integer quantum number / = 0, 1, 2, • • • . Setting Ii = I2 = I, we find the energy 
eigen-values [/(/ + !) — ni^]/2I + m^/ll^ which appear with multiplicity 1 if m = 
and multiplicity 2 if m 7^ 0. On the other hand, the Killing vector fields // are differ- 
ential operators of three angular variables {6, cp, xp) that parametrize S^, and, as such, 
they are the generators of SU{2) in the regular representatiorl^ Fourier analysis on 
SU{2) decomposes the regular representation into direct sum of unitary irreducible 
representations with all / = 0, 1/2, 1,3/2, •• • including half-integer ones. Further- 
more, each irreducible components appears with multiplicity equal to its dimension 
2j + 1 (see, for instance, | |98| and in particular p. 57-60). This, in turn, determines the 
multiplicity of eigen-values of the Laplace operator A, and, hence, of the operators T± 
which commute with it on Berger spheres. The resulting spectrum and multiplicities 
differ from those of rigid rotor despite of the similarities. 

Returning back to expression (|D.13|) , we find that the eigen-values z of the operator 
Q must satisfy the quadratic equation 

z 1 
+ - = (;■ ± m)(; ^ m + 1) + ^m(m ^ 1) (D.17) 

for all y = 0, 1 /2, 1, 3 /2, • • • . Since the two sign options are mutually related by send- 
ing mto —m, we may choose the lower sign (without loss of generality) and obtain the 



■^^The more general case h 7^ h 7^ h corresponds to a totally anisotropic Bianchi IX metric, which 
can not be solved exactly in closed form. 

^^The regular representation of a group G is defined as U{g)f{x) = f{gx) on the space of all square 
integrable functions on the group manifold endowed with the Haar measure. It satisfies U{gi)U{g2) = 
U{g\g2) and it is unitary. For SU{2), the local coordinates on the group manifold are x' = (9, q), ip) and 
gx is the left action of the group generated by the corresponding left-invariant vector fields fj. Similar 
considerations apply to the regular representation U{g)f{x) = f{xg^^) acting on the other side, which 
also satisfies U{gi)U{g2) = U{gig2) and it is unitary. 
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eigen-values of Q 



^ (^-1 ± ^J4:Sm - m)(; + m + 1) + (2m + 1)2^ (D.18) 



that appear with multiplicity 2j + 1. Both sign options are admissible in (|D.18|) as long 
as j 7^ m, in which case both components of the eigen-spinor Y(x) are non-zero. The 
choice j = m is special and requires separate investigation. In this case, the eigen- 
value problem for Q is solved directly as 



Q 



( 



\ 



\]'-] 







V 



\]'-] 







i>\ 



(D.19) 



Comparison with (|D.18|) shows that only the upper sign yields the correct eigen-value 
z = j/S. Likewise, one finds 



Q 
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> J 




V 





(D.20) 



showing that the eigen-value z = j/5 arises with multiplicity 2(2/ + 1). 

Setting p = j + m + l and q = j — m, we find that the spectrum of the Dirac operator 



(|D.12|) on Berger spheres is given by 



(D.21) 



as long as q ^ 0. The multiplicity of these eigen-values is p + q = 2j + 1 for each pair 
{p,q)- The additional eigen-values arising for q = are 



Co 



4: 2S 



p G N 



(D.22) 



and their multiplicity is 2p = 2(2; + 1). Thus, the complete spectrum of the Dirac 
operatoii§ is provided by ^± and forming a two-dimensional state lattice labeled 
by the quantum numbers {p,q)- The eigen-values ^+ and are positive definite, 
whereas ^- can also become zero or negative by tuning the parameters. Finally, we 
note that the eigen-values ^± simplify for p = q as they depend linearly on S, i.e., 
C± =^/4±p. 

Special cases: Three special cases arise from the general discussion that are worth 



^^Compared to |39|, there is a discrepancy by an overall factor of 2 in the allowed values of ^ that 
results from the different normalization of the SU{2) one-forms a^, but this is irrelevant. 
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noting. They reduce the general problem to simpler ones that are better known in the 
literature and can be used for comparison and consistency checks. These special cases 
are also important for understanding the behavior of the eigen-values under spectral 
flow as 5 varies from to oo. 

First, we consider 5 — 1 that corresponds to the homogeneous and isotropic con- 
stant curvature metric on S^. It can be easily seen by little combinatorics that the 
spectrum of the Dirac operator becomes 

1 

5^1: ^ = ±-(2n + 1) , n e N with multiplicity n(n + 1) . (D.23) 

In this case the eigen-values are equally distributed to positive and negative values 
and there is no spectral asymmetry in the problem. 

Next, we consider 5 — Q that corresponds to a fully squashed three-sphere along 
its third principal axis. The resulting configuration is non-singular as it exhibits no 
curvature singularity. It can be related to a two-dimensional sphere endowed with 
constant curvature metric. Then, the eigen-values become infinite and the same is 
true for which tend to ±oo, respectively, as long asp ^ c\. Finite spectrum arises 
only when p — q = n and it is given by 

3 = 0: ^ = ±n, wgN with multiplicity In . (D.24) 

This is the spectrum of the Dirac operator for the round which is equally distributed 
to positive and negative values and shows no spectral asymmetry However, there is 
spectral asymmetry in the three-dimensional problem, because provide the excess 
of positive eigen-values that become infinite as collapses to (it explains later why 
the //-invariant of the three-dimensional Dirac operator does not vanish for 5 — 0)). 

Finally, we consider the limiting case 5 ^ oo that yields a singular metric on as 
the sphere becomes infinitely stretched along its third principal axis. Then, all eigen- 
values become infinite exhibiting a universal behavior, 

(5 — )■ 00 : ^ ~ 4 ^ith infinite multiplicity . (D.25) 
D.2 Harmonic spinors and level crossing 

Harmonic spinors are normalizable zero modes of the Dirac operator. They appear 
when ^_ = 0, i.e., when there are positive integer solutions {p, q) to the equation 

5^^1 ^J^^pqS^ + ip- qf ■ (D.26) 
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Thus, a necessary condition for the existence of harmonic spinors is the following 
bound for the anisotropy parameter 5 of Berger spheres, 

^ > 4 , (D.27) 

in which case the metric becomes sufficiently stretched to allow for negative curvature. 
This is in accord to Lichnerowich's formula (|B.9|) for the square of the Dirac operator, 
which now takes the form (27' D,)^ = — D/D' + R/4. It tells us that if the Ricci scalar 
curvature is positive definite the Dirac operator will not admit zero modes. Zero 
curvature is also excluded in this case, as it corresponds to the special value 5 = 2 that 
falls well below the bound (|D.27|) for the existence of harmonic spinors. 

Note that the harmonic spinors and their multiplicities are invariant under rescal- 
ing of the metric on S^. They have clear geometrical interpretation but no topological 
origin. Thus, for each pair of positive integers that solve equation (|D.26|) there is 
an associated zero mode with multiplicity p + cj. The simplest ones arise for p = q in 
which case 3 = 4p. For p = q = 1, in particular, we have ^ = 4 and the space of har- 
monic spinors has dimension 2. More generally, for ^ = 4p the dimension of the space 
of harmonic spinors is at least 2p. An example that exhibits two different solutions 
arises for S = 260 in which case p = q = 65 and p = 528, q = 8 are both admissible 
choices and the multiplicity of zero modes turns out to be 2 • 65 + 2 • (528 + 8) = 1202 
in total (noting that p = 8, q = 528 has also to be counted). Conversely, the values of S 
that allow for harmonic spinors are provided by 8pq + 2y^l6p^q^ + (p — q)^ for each 
pair of integers {p,q), without restriction, and they are generally irrational numbers 
in the interval [4, 00) . 

The possibility to have zero modes implies level crossing under the spectral flow, 
which is very important for the integrated form of the axial anomaly. Let us examine 
how the spectrum of the Dirac operator changes as S varies. Note that the eigen- 
values ^+ and ^0 given by (|D.21|) and (|D.22|) , respectively, are always positive irre- 
spective of S. On the other hand, ^_ are all negative for ^ < 4, but for ^ > 4 some 
become positive or zero. The first level crossing occurs when the anisotropy parame- 
ter passes the critical value ^ = 4. More level crossings occur when S becomes bigger 
and bigger. As soon as an eigen-value crosses zero, it will stay positive for all higher 
values of S and never cross back to negative values, according to the general formulae. 
Eventually, as ^ — > 00, all eigen-values turn positive leaving no negative ones behind. 

Fig. 6 provides the spectral flow of the eigen-values ^_ as functions of S for the 
lowest values of the positive integers p and q. Obviously, all eigen-values depend 
linearly upon S whenever p = q; these lines pass from the values —1, —2, . . . when 
S = (they are the negative eigen-values of the Dirac operator on S^), whereas all 
other spectral lines are pushed to infinite as collapses to S^. Here, we plot the 
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results for {p,q) = (1,1), (1,2), (1,3), (2,2), (1,4) and (2,3) for which level crossing 
occurs at ^ = 4, 5.67, 6.95, 8, 8.03 and 9.80, respectively, where harmonic spinors make 
their appearance. Similar plots arise for all higher values of (p, q). Obviously, the lines 
(p, q) and {q, p) are identical. Note that all curves reach asymptotically the slope line 
^ = 3/4: thanks to the universal behavior (|D.25|) as ^ — > oo. We have omitted from the 
figure the spectral flow of the eigen-values ^+ and as they never cross zero, and, 
hence, they are uninteresting for our purposes. 



US) 




Figure 6: Level crossing of the modes ^_ (S) from negative to positive values 



D.3 The ^-invariant 

The results for the spectrum can been used to compute the ^/-invariant of the Dirac 
operator on Berger spheres. This invariance measures the asymmetry between the 
number of positive and negative modes in the spectrum, and, as such, it vanishes 
when (5 = 1. This is seen by computing the //-invariant in closed form. The result is 
very important for evaluating the index of the Dirac operator on gravitational instan- 
ton backgrounds by the Atiyah-Patodi-Singer index theorem ||T9l (see also l|T3ll20| ). 

First, we examine the case 3 < 4: which is the simplest as it does not involve level 
crossing. There are infinitely many positive modes ^+ and and infinitely many 
negative modes ^- . Considering the difference of the spectral Riemann zeta-function 
between the positive and negative eigen-values, excluding zero modes, 

f]^is)= (sign 1^1 (D.28) 

eigenvalues 
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allows to define the ;y-invariant of the Dirac operator by appropriate analytic contin- 
uation at s = as 

VD = 7d(0) . (D.29) 

The //-invariant is ill-defined without the zeta-function regularization as it is the dif- 
ference of two infinite numbers. Note that under rescaling of the metric the spectrum 
of the Dirac operator scales uniformly but j/d remains invariant since it is evaluated at 
s = 0. 

Then, taking into account the multiplicities of the various eigen-values, we have 
explicitly that 

Vd{s)= E^(p + ^)[(f + x) '-(-^ + x) '] + E2p(f + p) ' (D.30) 



p,q>0 

setting 



X = JU^pq + {p- qf (D.31) 



for notational convenience. The first term refers to the contribution of the second 
to ^_ and the third to ^o- The first two terms are a bit tricky to evaluate, whereas the 
third one is more straightforward. We will present the essential tricks and details since 
they are also needed in the main text to evaluate more complicated sums. 

The computations are most easily done by expanding all terms in power series of 
5. For the last one we have 

1 f. sS^ s(s + l)c54 \ ^^^^^ 



that yields the following expansion in terms of Riemann zeta-functions, 

p>0 



7o(s) = E 2P (f + p) = 2^(s - 1) - ^'s^(s) + l^Hs + l)as + !) + •••. (D.33) 



The higher order terms are irrelevant because they vanish at s = and they are omit- 
ted (they contain the factor s^(s + n) for all integers n > 2). Since 1) = —1/12, 
^(0) = -1/2 and s^(s + 1) equalfl 1 at s = 0, we obtain immediately 

W) = -l + j- (D-34) 



^^Recall that ^(s + 1) has a simple pole at s = with residue 1 since the Riemann zeta-function 
satisfies the functional relation ^(s + 1) = 2^+^7T''sin(7T(s + l)/2)r(— s)^(— s). The pole structure at 
s = follows from the identity (— s)r(— s) = r(l — s) and r(l) = 0. 
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Next, we consider the power series expansion 

2^ J X« V 2X 8X2 + 48X3 ^ ■ ■ 
that yields 

7+(s)-7_(s)= E^(p + ^)[(| + x) '-(-^ + x) 

= -^'sf - ^S's{s + l)(s + 2)/ + . . . , (D.36) 



where 



fu)= y E±l= y 111 (D37) 



The function /(s) converges absolutely for Res > 3/2 on the complex s-plane, and, 
therefore, the omitted terms in (|D.36|) are irrelevant as they vanish at s = 0. Fur- 
thermore, the meromorphic continuation of /(s) to the entire complex plane has only 
simple poles whose residues at s = 1/2 and s = 3/2 turn out to be 

^2-1 1 
res/(s)|s=i/2 = , res/(s)|s=3/2 = ^ . (D.38) 

They provide the values of (s/2)/((s + l)/2) and (s/2)/((s + 3)/2) at s = 0. All 
other terms appearing in the expression (|D.36|) vanish at s = 0, and, thus, we obtain 

7+(0)-J_(0) = --^. (D.39) 



The mathematica 
mation formulaj^^ 



properties of /(s) can be studied using the Euler-Maclaurin sum- 
that converts infinite sums into integrals - in this case of two vari- 



ables (see, for instance, ||99| ). 

^"^U a function f{x) has continuous derivatives on the interval [0, N], then the sum of f{n) with 
« e N is converted into integral as follows, 

^/(n)= / dxf{x) + Uf{N)+m)+j:^(f^^'^--'\N)-f(^'-'\0))+n. 

n=0 ■'^ ^ k=l 

Here, Bj/t are the Bernoulli numbers (B2 = 1/6, B4 = -1/30, Bg = l/42,etc) and/(2'^-^)(x) denotes the 
{2k — l)-th derivative of the function f{x). The remainder (error) term TZ is normally small for suitable 
value of p and it is bounded as 



l^l< (I^X (D.40) 
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There is no point in presenting the technical details of the residue formulae since 
they are somewhat involved. We refer the interested reader to the comprehensive 
report IllOOl for the relevant computations as well as certain generalizations to higher 
dimensional Berger spheres. We only note here that when ^ = 1 the computation of 
the residues can be done by elementary methods. In this case /(s) simplifies and it is 
expressed in terms of Riemann's zeta-function as follows. 

Then, it can be easily seen that /(1/2) = 1) — C^(O) = 5/12 is finite, and, hence, 
there is no pole, i.e., the residue at s = 1/2 vanishes. Evaluating/ (3/2) = ^(1) — ^(2) 
shows that there a simple pole at s = 3/2 given by the pole of ^(2s + 1) at s = whose 
residue is 1/2 (it is convenient to compare /(s + 3/2) with ^(2s + 1) at s = 0). These 
results agree with the general formulae (|D.38|I for ^ = 1. 

Finally, putting all these together, we add up the individual terms /o(0) + /+(0) — 
1- (0) and arrive at the following result for the //-invariant of the Dirac operator on 
Berger spheres with anisotropy parameter ^ < 4, as in ||39B , 

ni, = -\{5^-l?. (D.42) 

It can be readily seen that j/d vanishes for the isotropic metric (5 = 1, as expected. Also, 
in the fully squashed limit (J = 0, //d=— 1/6 counting the asymmetry of positive and 
negative modes that become infinitj^ (the finite eigen-values are equally distributed 
to positive and negative values). 

At this point, it is instructive to compare the //-invariant with the gravitational 
Chem-Simons action [|51| 

Wcs = Tr ^dco + '^co ^ CO ^o?j (D.43) 

evaluated for the class of homogeneous metric on S^. Using the Bianchi IX form of the 
metric (|D.2[) and the expressions for the connection one-forms oo^ j derived earlier, we 



find (up to a sign convention that depends on the orientation of S^) 

(D.44) 

The result is scale invariant, as required on general grounds. Then, choosing 71 = 



Wcs = 1671^ 



1 

1 + 7. (71 + 72 - 73)(7i - 72 + 73)(7i - 72 - 73) 

2717273 



■^^For ^ = this asymmetry is accounted by the excess positive modes ^q- Their zeta-function regu- 
larization yields the simple expression ?/d(s) = 2f (s — 1). Then, setting s = 0, we find by elementary 
means that ;/d = 2^( — 1) = —1/6 as required on general grounds. 
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72 = 1 and 73 = d'^, we obtain immediately 

^Wcs(Berger) = \{p^ - 2^' + 2) = -^d + ^ (D.45) 

that relates W^s and //d foi" Berger spheres. This is a special case of a more general re- 
lation among the Chem-Simons action and the //-invariant found in the context of the 
Atiyah-Patodi-Singer theory IIT9I1 (for a compact three-manifold without boundaries 
it suffices to consider a spin four-manifold bounded by it, extend the product metric 
near the boundary and integrate the Pontrjagin form relative to this metric). 

Next, we examine the case (5 > 4 that involves level crossing. Every time an eigen- 
value I,- crosses from negative to positive values, //d jumps by 2 to //d + 2, since the 
spectral asymmetry, as calculated for ^ < 4, changes by 2. Let us denote by 



C{5,) = {ip,q) G N2; Si = 2^^vq5l + (p - qf) (D.46) 

the set of positive integers {p,q) that account for harmonic spinors at each one of the 
special values < ^- This set has at least two elements associated to the two harmonic 
spinors at 5c = 4. Since the zero modes at 5c have multiplicity p + q, the quantity 



Sc<S 



ip,q)eCiSc 



provides the total number of modes that crossed from negative to positive values as 5 
was varying from ^ < 4 to any given value ^ > 4. Thus, for ^ > 4, the //-invariant of 
the Dirac operator for Berger spheres is shifted by twice the number of these modes 
and equals 

//D = -^(^'-l)' + 2S(^) (D.48) 



provided that 5 does not assume the special values (|D.26|) . 

On the other hand, if 5 is tuned to any one of the values (|D.26|) , including 5 = 4:, the 
corresponding zero modes have to be removed from the spectrum while computing 
the sum 7/d(s). We also have to add 1 for each one of these modes (i.e., shift t/d by 
the dimension of the space of harmonic spinors arising at those values of 5), and, 
consequently, the net effect is zero. Thus, there is no change of //d in this case. Changes 
can only occur when there is level crossing. 

Summarizing, the formula (|D.48|) for the //-invariant is the most general one and 
it is valid for all values of the anisotropy parameter 5 of Berger spheres. The result is 
finite for all finite values of 5 and it only becomes infinite in the extreme limit ^ — ^ 00 
where the geometry becomes singular. In the latter case there are no negative modes 
left behind and the spectral asymmetry is infinite, as expected. 
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